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ABSTRACT 


The  statistics  of  acoustic  signals  propagated  to  long  ranges  in  the  ocean 
are  investigated  in  detail  in  this  thesis.  The  phase  random  model  of  multipath 
propagation  is  extended  to  include  finite  bandwidth  and/or  modulated  sources  as 
well  as  multiple  source  configurations.  The  theoretical  analyses  include  the 
derivation  of  many  new  probability  density  functions  for  these  new  cases  as  well 
as  for  the  single  narrowband  source. 

y  ; 

The  probability  density  function  for  A,  the  time  rate  of  change  of  the  level 
in  decibels  for  the  single  narrowband  source  is  derived.  P- (A )  depends  only  upon 
v2.  the  single  path  mean  square  phase  rate,  which  can  be  related  to  certain  ocean 
, dynamical  processes.  The  analysis  of  finite  bandwidth  and/or  modulated  sources 
reveals  that  the  amplitude  and  amplitude  rate  densities  (including  P: (A ) )  are 
independent  of  the  finite  bandwidth  and  modulation  effects,  but  the  density  of 
the  time  rate-of-change  of  the  multipath  phase  |is  sensitive  to  these  effects. 
Thus,  fitting  P-(A)  to  histograms  from  data  to  find  is  the  preferred  method 
for  determiningAthis  important  parameter.  Bandwidth  effects  in  $can  be  neglected, 
however,  if  B<<2v  where  B  is  the  signal  bandwidth.  The  analysis  also  reveals  a 
potentially  powerful  technique  for  determining  parameters  of  the  modulation  or 
bandwidth  of  a  source  from  the  received  multipath  signal. 

The  analysis  of  multiple  sources,  applicable  to  noise  problems,  includes 
important  approximations  to  densities  which  are  intractable  analytically,  and 
would  involve  significant  computer  time  to  solve  exactly.  In  addition  to  studies 
of  the  amplitude  densities,  significant  progress  has  been  made  in  solving  for 
the  amplitude  rate  densities  and  the  joint  densities  of  amplitude  and  amplitude 
rate.  * 

In  addition  to  providing  valuable  confirmation  of  much  of  the  theoretical 
analysis,  a  computer  simulation  of  phase  random  multipath  propagation  also  con¬ 
firms  that  for  N  ^  4  paths  phase  random  multipath  conditions  begin  to  closely 
approach  the  asymptotic  conditions  for  N  +  «. 

Data  at  220Hz  and  406Hz  received  by  drifting  sonobuoys  in  the  Atlantic  at 
approximately  300  km  in  range  were  analyzed.  Values  of  v2  obtained  support  an 
internal  wave  model  for  the  relevant  dynamical  process.  The  modulation  theory 
uncovered  a  heretofore  unrecognized  modulation  in  the  data  due  to  an  error  of 
the  Doppler  tracking  system.  Predictions  of  crossing  rates  including  this  modu¬ 
lation  effect  are  in  good  agreement  with  the  data. 

Other  data  at  15Hz  and  33Hz  propagated  to  ranges  between  250  km  to  450  km 
in  the  Pacific  in  which  deliberate  modulation  was  introduced,  once  again  provide 
excellent  confirmation  of  the  theory.  Measured  values  of  v2  vary  significantly 
from  run  to  run  and  are  not  consistent  with  an  internal  wave  model,  indicating 
some  other  mechanism  (i.e.,  tidal,  rough  scattering)  must  account  for  the  fully 
saturated  phase  random  nature  of  the  data.  The  technique  for  determining  modula¬ 
tion  parameters  was  used,  and  for  the  76  modulated  runs  analyzed,  the  average 
error  in  determining  the  actual  bandwidth  of  the  modulation  of  the  source  from 
the  received  multipath  signal  was  8%. 
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ABSTRACT 


The  statistics  of  acoustic  signals  propagated  to  long 
ranges  in  the  ocean  are  investigated  in  detail  in  this 
thesis.  The  phase  random  model  of  multipath  propagation 
is  extended  to  include  finite  bandwidth  and/or  modulated 
sources  as  well  as  multiple  source  configurations.  The 
theoretical  analyses  include  the  derivation  of  many  new 
probability  density  functions  for  these  new  cases  as 
well  as  for  the  single  narrowband  source. 

• 

The  probability  density  function  for  A,  the  time 
rate  of  change  of  the  level  in  decibels  for  the  single 
narrowband  source  is  derived.  Pi  (A)  depends  only  upon  v2, 
the  single  path  mean  square  phase  rate,  which  can  be 
related  to  certain  ocean  dynamical  processes.  The  analysis 
of  finite  bandwidth  and/or  modulated  sources  reveals  that 
the. amplitude  and' amplitude  rate  densities  (including 
Pl(A))  are  independent  of  the  finite  bandwidth  and 
modulation  effects,  but  the  density  of  the  time 
rate-of-change  of  the  multipath  phase  $  is  sensitive  to 
these  effects.  Thus,  fitting  P:(A)  to  histograms  from 
data  to  find  v2  is  the  preferred  method  for  determining 
this  important  parameter.  Bandwidth  effects  in  $  can  be 
neglected,  however,  if  3  <<  2v  where  B  is  the  signal 
bandwidth.  The  analysis  also  reveals  a  potentially 
powerful  technique  for  determining  parameters  of  the 


modulation  or  bandwidth  of  a  source  from  the  received 
multipath  signal. 

The  analysis  of  multiple  sources,  applicable  to  noise 
problems,  includes  important  approximations  to  densities 
which  are  intractable  analytically,  and  would  involve 
significant  computer  time  to  solve  exactly.  In  addition 
to  studies  of  the  amplitude  densities,  significant  progress 
has  been  made  in  solving  for  the  amplitude  rate  densities 
and  the  joint  densities  of  amplitude  and  amplitude  rate. 

In  addition  to  providing  valuable  confirmation  of 
much  of  the  theoretical  analysis,  a  computer  simulation 
of  phase  random  multipath  propagation  also  confirms  that 
for  N  ^  4  paths  phase  random  multipath  conditions  begin  to 
closely  approach  the  asymptotic  conditions  for  W  ■+  ». 

Data  at  220Hz  and  406Hz  received  by  drifting 
sonobuoys  in  the  Atlantic  at  approximately  300  km  in  range 
were  analyzed.  Values  of  v2  obtained  support  an  internal 
wave  model  for  the  relevant  dynamical  process.  The 
modulation  theory  uncovered  a  heretofore  unrecognized 
modulation  in  the  data  due  to  an  error  of  the  Doppler 
tracking  system.  Predictions  of  crossing  rates  including 
this  modulation  effect  are  in  good  agreement  with  the  data. 

Other  data  at  15 Hz  and  3 3 Hz  propagated  to  ranges 
between  250  km  to  450  km  in  the  Pacific  in  which  deliberate 
modulation  was  introduced,  once  again  provide  excellent 
confirmation  of  the  theory.  Measured  values  of  v2  vary 
significantly  from  run  to  run  and  are  not  consistent  with 
an  internal  wave  model,  indicating  some  other  mechanism 
(i.e.,  tidal,  rough  scattering)  must  account  for  the  fully 
saturated  phase  random  nature  of  the  data.  The  technique 
for  determining  modulation  parameters  was  used,  and  for 
the  76  modulated  runs  analyzed,  the  average  error  in 
determining  the  actual  bandwidth  of  the  modulation  of  the 
source  from  the  received  multipath  signal  was  8%. 
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v*  .007  rad/sec.  A*  .007  rad/sec, 
and  Equation  (1.176)  . 

Histogram  of  computer  generated 
samples  of  $  with  uniform  frequency 
modulation  of  the  source, 
v=  .007  rad/sec,  A=  .021  rad/3ec, 
and  Equation  (1.176)  . 

Histogram  of  computer  generated 
samples  of  with  uniform  frequency 
modulation  of  the  source, 

\>».007  rad/sec,  A».070  rad/sec, 
and  Equation  (1.176). 

Histogram  from  Record  447  for 
X,  A,  and  $,  and  the  least 
Chi-square  fit  of  Equations  (1.4), 

(1.7),  and  (1.176),  respectively. 
Results  for  220  Hz  are  on  the  left, 
and  406  Hz  on  the  right. 

Histograms. from  Record  448  for 
X,  A,  and  <£,  and  the  least 
Chi-square  fit  of  Equations  (1.4), 

(1.7),  and  (1.176),  respectively. 
Results  for  220  Hz  are  on  the  left, 
and  406  Hz  on  the  right. 

Histograms . from  Record  449  for 
X,  A,  and  $,  and  the  least 
Chi-square  fit  of  Equations  (1.4), 

(1.7),  and  (1.176),  respectively. 


Results  for  220  Hz  are  on  the  left, 

and  406  Hz  on  the  right.  223 

Histograms . from  Record  424  for 
X,  A,  and  and  the  least 
Chi-square  fit  of  Equations  (1.4), 

(1.7),  and  (1.176),  respectively. 

Results  for  220  Hz  are  on  the  left. 


and  406  Hz  on  the  right.  224 

The  histograms  of  x  and  p  for  two 
different  sources  from  Record  447, 
plotted  with  Equations  (3.1)  and 
(3.2)  respectively.  228 

The  histogram  of  %  for  two  different 

sources  from  Record  449,  plotted 

with  Equation  (1.95)  with  L*2.  229 

The  histogram  of  p  for  two  different 

sources  from  Record  447,  plotted 

with  Equation  (1.114).  230 

The  histogram  of  A  for  two  different 

sources  from  Record  447,  plotted 

with  Equation  (1.123).  231 


Ratio  of  phase  period  to  amplitude 
period,  versus  amplitude,  for 
Records  447,  448,  442,  and  424, 
respectively.  Equation  (3.6)  is  the 
theoretical  ratio  for  a  non- 
modulated  narrowband  source.  The 
correction  factors  v./v*  (Table  XV) , 
and  Equation  (1.209)A  ? 

(Gaussian  modulation)  have  been 

applied  to  obtain  the  corrected 

curves  for  220  Hz,  and  406  Hz 

accounting  for  the  modulation.  237-240 

The  source  frequency  function  for 
the  modulation  induced  in  the  CASE 
experiment,  reproduced  from 
Reference  [32  3  . 
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The  source  signal  power  spectrum  of 
the  modulation  used  in  the  CASE 
experiment/  reproduced  from 

Reference  [32].  247 

The  error  in  percent  versus  number 
of  runs  in  determining  the  bandwidth 
of  the  modulation  of  the  CASE  data. 

The  shaded  portions  represent  those 

runs  in  which  v  “  A.  248 

• 

The  time  series  and  histogram  of  A 

is  plotted  with  the  best  fit  of 

Equation  (1.7)  for  a  representative 

run  of  the  CASE  experiment  in  which 

A  »  v.  255 

The  time  series  and  histogram  of  $ 

for  the  same  run  shown  in  Figure  71 

is  plotted  along  with  the  best  fit 

of  Equation  (A8)  which  assumes  no 

modulation  (a)  and  the  best  fit 

of  Equation  (1.176)  which  includes 

the  effect  of  source  induced 

modulation  (b) .  256 

Histogram  of  v  for  all  the  runs 
analyzed  in  the  CASE  experiment. 

The  interval  is  .025  rad/sec  and 

the  black  bars  represent  33  Hz,  the 

clear  bars  15  Hz.  259 

Range  vs  frequency  for  <$z>  •  2tt. 

The. portion  of  the  graph  above  and 
to  the  right  of  the  line  is  the  fully 
saturated  region  based  on  internal 
wave  models.  The  CASE  experiment  is 
indicated  by  the  dotted  lines  and  the 
Eleuthera  data  by  the  circles.  260 

A  plot  of  v  vs  range  rate.  The 

circles  represent  15  Hz,  the  X's  33  Hz.  262 

A  plot  of  v  vs  range.  The  circles 
represent  15  Hz,  the  X’s  33  Hz.  263 
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A  plot  of  v  vs  time  of  day.  The 
circles  represent  15  Hz,  the  X's 
33  Hz. 

A  plot  of  v  vs  receiver  for  15  Hz. 
The  lines  connect  values  of  v 
obtained  simultaneously  at  each 
receiver.  The  start  times  for  each 
run  (ZULU)  are  indicated  to  the 
left  of  the  receiver  1  values.  The 
N  indicates  there  was  zero  range 
rate  for  that  run. 

Plot  of  v  vs  receiver  for  33  Hz. 

The  lines  connect  values  of  v 
obtained  simultaneously  at  each 
receiver.  The  start  times  for  each 
run  (ZULU)  are  indicated  to  the 
left  of  the  receiver  1  values.  The 
N  indicates  there  was  zero  range 
rate  for  that  run. 

A  plot  of  v  vs  range  rate  for  the 
33  Hz  signal  measured  at  receiver  1. 

A  plot  of  v  vs  range  rate  for  the 
33  Hz  signal  measured  at  receiver  2. 

A  plot  of  v  vs  range  rate  for  the 
33  Hz  signal  measured  at  receiver  3. 

A  plot  of  v  vs  range  rate  for  the 
15  Hz  signal  measured  at  receiver  1. 

A  plot  of  v  vs  range  rate  for  the 
15  Hz  signal  measured  at  receiver  2. 

A  plot  of  v  vs  range  rate  for  the 
15  Hz  signal  measured  at  receiver  3. 

The  computer  generated  histogram  for 
N»2  paths,  plotted  with  the  exact 
pdf,  Equation  (C8) ,  and  the 
limiting  pdf  (N  ■*  *)  , 

Equation  (A2) . 
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C2  The  computer  generated  histogram 

for  N  -  3  paths,  plotted  with  the 
exact  pdf,  Equation  (C23),  and  the 
limiting  pdf  (N -♦'«),  Equation  (A2)  . 

C3  The  computer  generated  histogram 

for  N  -  4  paths ,  plotted  with  the 
limiting  pdf  Equation  (A2) . 

C4  The  computer  generated  histogram 

for  M«2  paths,  plotted  with  the 
exact  pdf.  Equation  (C9) ,  and 
the  limiting  pdf  (N  ♦  «) , 

Equation  ( A1 ) . 

C5  The  computer  generated  histogram 

for  N*4  paths,  plotted  with 
the  exact  pdf,  Equation  (C22) , 
and  the  limiting  pdf 
Equation  (Al) . 

C6  The  computer  generated  histogram 

for  N-4  paths,  plotted  with 
the  limiting  pdf.  Equation  (Al) . 

C7  The  computer  generated  histogram 

for  N«2  paths,  plotted  with 
the  exact  pdf,  Equation  (CIO) ,  and 
the  limiting  pdf  (.V  ■*>  ®)  , 

Equation  (A3)  . 

C8  The  computer  generated  histogram 

for  A/ -  3  paths,  plotted  with  the 
exact  pdf,  Equation  (C24) ,  and 
the  limiting  pdf  (N +  ») , 

Equation  (A3)  . 

C9  The  computer  generated  histogram 

for  N»4  paths,  plotted  with  the 
limiting  pdf,  Equation  (A3) . 

CIO  Computer  generated  histogram  for 

N»2  paths,  plotted  with  the  exact 
pdf.  Equation  ,  (C21) ,  and  the 
limiting  pdf  (M  •*■  «)  ,  Equation  (1.4). 
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Comparison  plot  of  the  effect  of 
amplitude  variation  on  the  pdf  for 
p.  Equation  <A1)  is  plotted  with 
Equation  (D3)  with  various  levels  of 
uncertainty,  W-5  paths. 

A  comparative  plot  of  the  effects 
of  uncertainty  in  the  single  path 
amplitude  on  the  density  of  the 
level.  Equation  (A3)  no  uncertainty, 
is  plotted  with  Equation  (D5)  with 
various  values  of  uncertainty. 
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LIST  OF  SYMBOLS  AMD  NOTATION 


maximum  excursion  of  the  carrier  in  Hz, 

2A  is  therefore  the  bandwidth  of  the 
modulation 

the  inverse  of  the  mean  of  x  (a  -  1/u) 
bandwidth  of  the  signal 
the  inverse  of  2c|v 

product  of  the  modulation  amplitude  and 
modulation  frequency  for  sinusoidal  phase 
modulation 

gamma  function 

incomplete  gamma  function 

incomplete  gamma  function 
[y (a,z)  -  T(a)  -  T(a,z) ] 

coefficient  of  skew 

coefficient  of  excess 

Euler's  constant  (»  .5772...) 

=  1* 3 •  5*  •  •  (2n-l)  ;  (-1)1!  =  1 

parabolic  cylinder  function 

defined  as  the  expected  value  of  random 
variable  x 

the  exponential  integral 

complete  elliptic  integral  of  the  second 
kind 

defined  as  10  log1Q  e  »  4.34... 

Gauss's  hypergeometric  series 
variables  of  integration 
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G(Pq) 

G(V 

9n 

Jv(z) 
Kv(2) 
K[tt/2,  y] 
L 
LO 

LE 

A 

• 

A 

M(t) 

M(t) 

M  (w) 

A 


amplitude  crossing  rate  at  axis  level  equal 
to  pQ,  or  the  conditional  expection  of 
|p|  given  p  -  pQ  (*  E[p|p  =  pQ] ) 

phase  crossing  rate 

single  path  phase  for  the  nth  path 

Bessel's  function  of  order  v 

the  modified  Bessel  function  of  order  v 

complete  elliptic  integral  of  the  first  kind 

number  of  sources 

number  of  sources  with  different  intensity 
and  different  mean  square  phase  rates  at 
the  receiver 

number  of  sources  with  equal  intensity  and 
equal  mean  square  phase  rates  at  the 
receiver 

level  in  dB's  of  the  short  time  average 
mean-square  pressure  (A  *  10  log1Qx) 

time  rate  of  change  of  A 

temporal  modulation  function 

time  rate  of  change  of  M(t) 

characteristic  function  of  random 
variable  x 


M  (w,a)  two-dimensional  characteristic  function  of 
xl,x2  the  joint  pdf  of  random  variables  Xt 

and  x2 

M(a,5,z)  the  confluent  hypergeometric  function 

[sometimes  written  <?(a,6,z)] 


N 


SI 


number  of  groups  of  sources,  each  group  with 
LE^  (i  ■  1,2 , . . . ,N)  sources  per  group 

number  of  independent  propagation  paths 
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U  the  long  time  average  mean  square  pressure 

at  the  receiver  for  a  single  source 
(y  -  2a*  =  E[X  1) 

v2  the  single  path  mean-square  phase  rate 

(v2  5  E[S2]V  n) 

P  (x)  probability  density  function  (pdf)  of 

x  random  variable  x 


P^  v  (Xi  »X?)  the  joint  pdf  of  random  variables  x-i 
xl'x2  1  1  and  x2 

P  ( x)  the  pdf  of  random  variable  x  for  N  groups 

NX 

(P  (s)  the  Laplace  transform  of  the  pdf  of 

x  random  variable  x 


r 


single  path  amplitude  at  the  receiver 


P 


short  time  average  root-mean-square  pressure 


p  time  rate  of  change  of  p 

a2  variance  of  phase  in  Gaussian  phase 

modulation 


a 


z 

1 


W 


v,u 


(2) 


X,Y 


one  half  of  the  long  time  average 
mean-square  pressure  at  the  receiver 
(a2  -  Nr2/2) 

Whittaker's  function 

quadrature  components  of  the  signal  envelope 


•  • 
X,Y 


time  rate  of  change  of  the  quadrature 
components 


#,(ot,8,Y;x,yy  degenerate  hypergeometric  series  in  two 
variables 


X  short  time  average  mean-sauare  pressure 

(X-  P2) 

X  time  rate  of  change  of  x 

$  multipath  phase 


„»*■ 

■*..  rax* 


multipath  phase  rate 
Euler  psi  function  ['Mx)  = 
digamma  function  [ ^ 1  ( x)  =  3/dx 


4> 

4»(x) 
'J'1  (x) 
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INTRODUCTION 


Acoustic  signals  propagated  to  long  ranges  in  the 
ocean,  tens  to  hundreds  to  thousands  of  kilometers,  via 
all  modes  of  propagation  including  surface  ducts,  the 
deep  sound  channel,  or  sea-surface  and  ocean-bottom 
reflections  exhibit  fluctuations  in  amplitude  and  phase 
which  are  now  recognized  to  be  dominated  by  the  multipath 
interference  of  the  acoustic  field.  The  statistics  of 
these  fluctuations  as  well  as  their  relationships  to  the 
dynamics  of  the  ocean  has  been  one  focus  of  recent  research 
in  under standing  this  important  physical  process.  The 
optimum  design  of  sonars  (e.g.,  the  receiver  operating 
characteristics) ,  underwater  communications  devices,  and 
in  fact  any  system  which  operates  via  acoustic 
transmission  in  the  sea  depends  upon  the  knowledge  of  the 
statistical  behavior  of  these  transmissions. 

The  recognition  of  the  dominance  of  the  multipath 
structure  on  the  statistics,  or  the  assertion  that  long 
range  multipath  acoustic  propagation  in  the  ocean  can  be 
modelled  as  a  phase  random  process  has  been  established 
only  within  the  last  ten  years,  although  the  phase  random 
process  or  random  walk  problem  has  been  under  study  since 
Rayleigh  (1880)  [1] ,  and  is  one  of  the  classical  problems 

of  mathematics  and  physics.  Bergmann  (1946)  [2]  was  among 


the  first  to  speculate  that  observed  fluctuations  in 
signal  intensity  might  be  attributed  to  the  interference 
of  many  paths  summing  in  random  phase.  Dyer  (1970)  [3] 
formally  applied  the  theory  of  a  phase  random  process  or 
random  walk  problem  to  long  range  acoustic  multipath 
propagation  in  the  ocean,  and  was  the  first  to  investigate 
the  statistics  of  log  transformed  variables.  Dyer  also 
showed  that  even  in  the  presence  of  scattering  randomness 
multipath  interference  would  dominate  the  statistics. 

This  research  in  fact  indicated  a  basic  shift  from  the 
scattering  models  of  earlier  research  which  are  more 
appropriate  for  high  frequencies  and  short  ranges  when 
multipath  effects  are  less  important.  Dyer  also  proposed 
a  model  of  distant  shipping  noise  based  upon  the  precepts 
of  phase  random  multipath  propagation  and  continued  this 
research  in  a  later  paper  (1973)  [4] ,  and  most  recently  in 

Mikhalevsky  and  Dyer  (1978)  [5] ,  results  of  the  latter 

being  included  as  part  of  this  thesis.  This  model, 
appropriate  to  distant  shipping  noise,  assumes  the  noise 
in  a  band  is  dominated  by  narrowband  lines  discrete  in 
frequency. 

Mark  (1972)  (6]  investigated  the  statistics  of  the 
multipath  propagation  of  finite  bandwidth  signals. 

Employing  a  systems  approach,  he  derived  general  expressions 
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for  the  mean  and  variance  of  the  received  signal  energy 
in  terms  of  the  correlation  function  of  the  time  varying 
impulse  response  of  the  medium.  He  also  showed  that  the 
gamma  or  Erlang  probability  density  function  is  often  a 
good  approximation  to  the  real  pdf  of  the  received 
energy.  Much  earlier,  Nakagami  [7]  had  noted  the 
utility  of  an  appropriately  transformed  Erlang  pdf  in 
approximating  the  densities  of  received  HF  electromagnetic 
radiation  undergoing  rapid  fading,  indicating  the 
broadness  of  scope  of  the  phase  random  model  and  its 
general  applicability.  All  these  efforts,  however, 
concentrated  on  the  amplitude  (or  related  quantities)  of 
the  signal  and  did  not  address  the  amplitude  rate  or 
phase  rate  of  the  signal. 

Longuet-Higgins  (1975)  [8]  in  connection  with  research 

on  random  sea  surface  waves  (another  phase  random  process) 
introduced  the  joint  pdf's  of  amplitude,  amplitude  rate, 
phase,  and  phase  rate  as  well  as  the  marginal  densities  to 
the  growing  body  of  knowledge  of  phase  random  processes . 

It  remained  for  Hamblen  (1977)  [9]  to  formally  extend  the 
phase,,  random  analysis  to  the  multipath  acoustic  propagation 
process,  incorporating  the  results  of  Longuet-Higgins  [8] 
and  also  S.O.  Rice  [10]  whose  extensive  research  on  noise 
statistics  were  also  applicable  to  the  long  range  acoustic 
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propagation  problem.  Hamblen  [9]  also  established  the 
dependence  of  the  probability  density  functions  for  a 
single  narrowband  acoustic  source  on  the  two  parameters 
a£,  one-half  the  mean  square  pressure  at  the  receiver, 
and  v2,  the  single  path  mean  square  phase  rate.  He  also 
verified  the  basic  results  with  data  from  an  ocean 
acoustic  propagation  experiment. 

Concurrent  with  the  development  of  the  phase  random 
model  of  multipath  acoustic  propagation,  much  research 
was  and  is  being  conducted  on  another  important  aspect 
of  the  problem,  namely  to  discover  what  ocean  dynamic 
processes  are  the  driving  mechanisms  and  how  parameters 
of  these  ocean  dynamic  processes  are  related  to  the 
parameters  of  'the  acoustic  field,  and  v2.  Most 
notable  perhaps  is  the  recent  research  of  Dyson,  Munk, 
and  Zetler  (1976)  [11]  who  have  proposed  a  theoretical 

model  relating  the  dynamics  of  internal  waves  in  the  ocean 
to  the  fluctuation  of  the  acoustic  field.  However,  this 
model  appears  to  have  serious  limitations  at  low 
frequencies.  In  fact,  little  research  has  been  reported 
on  the  low  frequency  cases. 

In  the  following  paragraphs  I  will  introduce  the 
research  reported  on  in  this  thesis.  There  are  three 
basic  areas  in  which  significant  progress  has  been  made  in 


understanding  the  statistical  behavior  of  multipath 
acoustic  propagation:  (1)  the  single  narrowband  source, 

(2)  finite  bandwidth  and/or  modulated  sources,  and 

(3)  multiple  sources  or  receptions.  An  extensive 
analysis  of  data  from  acoustic  experiments  in  the  ocean 

as  well  as  a  computer  simulation  of  phase  random  multipath 
acoustic  propagation  not  only  increase  our  confidence  in 
this  new  understanding,  but  reveal  new  information  for  low 
frequency  signals  about  the  driving  mechanism  of  v2,  the 
single  path  mean  square  phase  rate. 


The  Single  Narrowband  Source 

I  derive  for  a  single  narrowband  source  the  pdf's 
for  the  time  rate  of  change  of  the  sta  (short  time  average) 
mean  square  pressure,  ;<,  and  the  time  rate  of  change  of 
the  level  in  decibels.  A,  results  which  are  unique  to  this 
thesis. ■  The  pdf  for  A  is  independent  of  a2  and  is  a 
function  only  of  v2.  This  result  is  of  particular 
importance  as  it  affords  a  method  of  measuring  vJ  from 
ocean  acoustic  data  without  error  due  to  uncertainties  in 
the  signal  carrier.  Included  in  this  analysis  are  the 
joint  densities  P  . (x,x)  and  P,  ;  (A, A)  and  their 
characteristic  functions,  also  unique  to  this  thesis. 

This  analysis  in  fact  completes  the  family  of  first  order 


and  joint  pdf's  for  a  single  narrowband  source. 

Appendix  A  contains  a  list  of  all  the  first  order  pdf's, 
their  characteristic  functions,  means,  and  variances  as 
found  by  me,  and  earlier  as  found  by  others,  so  that  a 
complete  set  can  be  referred  to. 

Finite  Bandwidth  and/or  Modulated  Sources 

Many  acoustic  signals  of  interest  received  in  the 
ocean  have  bandwidths  which  are  not  narrow,  and  carriers 
that  may  not  be  stable  or  may  be  deliberately  modulated. 

In  applying  the  phase  random  model  of  multipath  acoustic 
propagation,  it  is  necessary  to  assume  that  the  spectrum 
of  the  received  signal  is  narrow,  and  that  in  homodyning 
the  signal  the  spectral  mean  is  zero  [8,9].  Clearly, 
the  signals  mentioned  in  the  beginning  of  this  paragraph 
would  violate  these  assumptions.  I  show,  however,  that 
the  amplitude  and  amplitude  rate  statistics  (including 
P^(A)  are  independent  of  finite  bandwidth  and/or  carrier 
modulation  effects.  The  multipath  phase  rate,  <j>,  is 
sensitive  to  these  effects,  and  it  is  the  pdf  for  <j>  which 
must  be  modified. 

The  pdf  for  0  is  in  fact  a  function  of  v2,  as  well  as 
the  bandwidth,  and/or  parameters  of  the  modulation. 
Bandwidth  effects  can  be  neglected  when  B  <<  2v,  which  is 


-29- 


therefore  a  criterion  for  what  is  meant  by  "narrowband" 

in  the  ocean.  Furthermore,  this  analysis  reveals  a 

method  of  separating  and  understanding  the  source  induced 

modulation  independent  of  the  ocean  induced  modulation 

or  vice  versa.  In  addition  to  solving  for  the  pdf's  of 
• 

<j)  in  the  presence  of  bandwidth  and/or  modulation,  I  solve 
for  the  crossing  rate  statistics  of  phase  for  these  cases 
as  well. 

Multiple  Sources  or  Receptions 

The  statistics  of  the  received  signal  amplitude  and 
amplitude  rate  when  there  is  multiple  source  structure 
depends  on  the  exact  nature  of  the  received  multi- source 
signal  and  the  analysis  performed  by  the  receiver.  I 
consider  two  basic  cases.  First,  I  assume  the 
receptions  (one  per  source)  are  disjoint  in  frequency  in 
the  analysis  band  and  can  be  separated  and  summed 
incoherently  (that  is',  without  concern  for  phase)  .  I 
assume,  therefore,  that  each  of  the  receptions/sources 
are  independent,  thus  the  analysis  band  should  not  include 
harmonics  of  a  signal  already  in  the  band.  This  type 
of  analysis  is  motivated  by  the  structure  of  distant 
shipping  noise  [4,5].  For  the  second  case,  I  consider 
the  receptions  to  be  at  or  so  close  in  frequency  that  they 
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must  be  summed  coherently.  Note  that  the  dividing  line 
between,  case  one  and  case  two  depends  upon  the  resolution 
o£  the  receiver.  Most  all  of  the  pdf's  I  derive  for  the 
amplitude  and  particularly  the  amplitude  rate  for  these 
cases  are  unique  to  this  thesis.  Where  I  have  been 
unable  to  derive  the  exact  pdf  analytically  or  in  which 
the  exact  solution  is  extremely  time  consuming  to  obtain 
even  with  the  aid  of  a  computer,  I  have  in  most  cases 
found  an  approximation  based  upon  Edgeworth's  series. 

This  approximation  is  shown  to  be  excellent  in  the  main 
lobe  of  the  density  but  performance  is  degraded  in  the 
tails. 

I  use  the  analysis  of  coherent  sources  to  model  the 
effect  of  ocean  ambient  noise  on  the  pdf's  for  a  single 
narrowband  source.  The  pdf's  are  expressed  in  terms  of 
the  SNR  (signal  to  noise  ratio) .  In  light  of  this 
analysis,  I  am  able,  as  well,  to  extend  with  only  slight 
modification  all  the  frequency  disjoint  multiple  source 
solutions  to  include  the  cases  when  both  coherent  and 
disjoint  source  structure  is  present  in  the  analysis  band 
or  for  noise  which  is  continuously  distributed  over  the 
passband. 

The  reader  is  forewarned  that  the  sections  of  this 
thesis  on  multiple  sources  (see  Table  of  Contents)  are 
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lengthy  and  in  many  cases  quite  tedious  due  to  the 
complexity  of  the  analysis.  Table  III  is  provided  (page  143) 
to  aid  in  understanding  the  organization  of  this  material 
and  to  aid  the  reader  in  finding  that  analysis  which 
is  most  pertinent  to  his  problem.  The  analysis  of  the 
statistics  of  multiple  sources  which  is  not  immediately 
motivated  by  any  current  ocean  acoustics  problem  of 
interest  is  presented  for  completeness  with  potential  for 
application  to  future  problems,  even  perhaps  in  areas 
unrelated  to  acoustics.  For  example,  the  solutions  for 
P  (x) #  using  the  terminology  of  Queuing  Theory,  are  in 

A 

fact  the  pdf's  for  the  interarrival  times  of  cascaded 
Poisson  processes  which  are  unique  to  this  thesis  and, 
to  my  knowledge,  not  to  be  found  in  Queuing  Theory 
literature. 

Computer  Simulation  and  Data  Analysis 

A  computer  simulation  of  phase  random  multipath 
propagation  is  developed  to  assist  in  and  to  check  the 
theoretical  analysis.  The  simulation  demonstrates  the 
independence  of  the  exact  nature  of  the  pdf  for 
(the  single  path  phase)  to  the  statistics  when  there  are 
at  least  four  propagation  paths  from  the  source  to  the 
receiver.  Of  importance  is  the  confirmation  of  the 
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generally  accepted  criterion  that  .V  _>  4  paths  results  in 
an  almost  fully  saturated  phase  random  process,  and  pdf's 
for  the  limiting  case  as  W  -*■  «°  suffice.  The  simulation 
allows  me  in  Appendix  C  to  compare  the  theoretical 
pdf's,  many  which  I  derived,  for  N  <  3  to  the  computer 
generated  histograms  with  excellent  results. 

I  analyze  data  from  two  ocean  acoustic  experiments. 
One  was  conducted  near  Eleuthera  in  which  two  CW  signals 
at  220Hz  and  406Hz  were  transmitted  approximately  300  km 
northeast  towards  Bermuda  and  received  by  drifting 
sonobuoys  [12] .  These  data  support  the  theoretical 
pdf's  derived  in  Chapter  1.  Furthermore,  the  modulation 
theory  uncovered  a  heretofore  unrecognized  modulation 
in  the  phase  rate  data  due  to  errors  in  the  sonobuoy 
tracking  system  that  has  dramatic  results  on  the 
statistics  for  the  phase  crossing  rates.  These  data  are 
also  consistent  with  a  model  of  phase  random  multipath 
propagation  resulting  from  the  interaction  of  the  acoustic 
field  with  internal  waves  [11] . 

I  also  analyze  data  taken  in  the  Pacific  in  1973 
known  collectively  as  the  CASE  experiment  [13]  in  which 
CW  signals  at  15Hz  and  33Hz  were  propagated  to  ranges 
varying  from  250  km  to  450  km.  These  data  were 
deliberately  modulated  and  the  predictions  derived 
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theoretically  in  Chapter  1  are  once  again  confirmed.  The 
analysis  I  used  was  capable  of  determining  the  bandwidth 
of  the  modulation  from  the  received  multipath  signal 
with  an  average  error  of  8%.  Where  previously  there  had 
been  some  problems  in  the  consistency  of  the  data  with  the 
phase  random  model,  the  modulation  theory  successfully 
lays  these  problems  to  rest.  These  data,  however,  are 
not  consistent  with  the  internal  wave  model  and  this 
analysis  has  brought  sharply  into  focus  a  shortcoming 
in  our  current  understanding  of  the  driving  mechanisms  of 
v2,  particularly  for  low  frequency  signals.  Correlations 
I  have  derived  appear  to  support,  though  tenuously  without 
additional  research,  either  a  rough  scattering  or  tidal 
mechanism  to  account  for  the  fully  saturated  nature  of 
the  CASE  data. 

In  Appendix  D,  the  effects  of  amplitude  parameter 
variation  are  discussed.  This  analysis  is  aimed  at 
uncovering  the  effects  of  temporal  variations  in  the 
totai  energy  of  the  signal  during  the  observation  period 
on  the  statistics  of  the  received  multipath  signal. 
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CHAPTER  1 


ANALYTICAL  TREATMENT 


1.1  Single  Narrowband  Source 

For  a  narrowband  signal  the  phase-random  model  of 
multipath  acoustic  propagation  predicts  that  the  sta  root- 
mean-square  pressure/  p,  is  a  Rayleigh  distributed  random 
variable  [3/9]/  its  rate,  J,  is  Gaussian  [9],  the  multipath 
phase,  <t> ,  is  uniform,  and  its  rate,  p,  is  distributed 
according  to  a  density  first  given  by  M.S.  Longuet-Higgins 
[8,9].  The  transformation  x  “  P2  yields  the  density  for 
the  sta  mean-square  pressure  which  is  exponential  [3 ] , 
and  following  Dyer  [3 ]  the  transformation  A  »  10  log  x 

i  a 

gives  the  density  for  the  level  in  decibels  which  is 
Log-Rayleigh.  To  complete  the  family  of  first  order 
densities  for  a  single  narrowband  source,  I  have  derived 
the  densities  for  x  and  A. 


1.1.1  Derivation  of  P; (x)  and 

For  phase  random  multipath  propagation  the  joint 
density  of  the  sta  rms  pressure  p  and  its  rate  p  is  given 
by  [8,9], 


P  •  (o,P) 

p ,  p 


- - -  exp  |  -  - 

ai:  L 


2j^: v2 


(1.1) 


and  N  is  the  number  of  propagation  paths,  r  is  the  single 
path  amplitude,  and  u  is  the  long-time  average  mean-square 
pressure. 

For  x  38  P2  we  have  x  =  2p £>.  To  find  the  pdf 
(probability  density  function)  for  x  1  first  solve  for 
the  cumulative  distribution  function  0.4!  of  x  by 
integrating  over  the  joint  density  of  p  and  $, 

Equation  (1.1). 


*<*<*> 


J  P 


PfP 


(p,£)dpd&  ,  x  >  0 


(1.3) 


p  —oo 


The  pdf  for  x  will  be  given  by  P<^(x).  Combining  this 
result  with  Equations  (1.1)  and  (1.3)  differentiating 
under  the  integral  sign  and  integrating  once,  I  obtain 
the  result 


4c,  v 


exp 


2C1ZV; 


-oo<  v  <® 


(1.4) 
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The  absolute  value  results  from  the  analysis  for  x<  0. 
Equation  (1.4)  is  the  Laplace  pdf  Q.4]  .  The  variance  of 
is  given  by 


a.2  »  Sa.^v2 
X  1 

The  pdf  for  A  is  solved  in  the  same  manner.  First  I 
find  the  pdf  for  y  where  y  *  lnp2  and  thus  J  *  2p /p, 
therefore. 


_3_ 

3<r 


CO 


0 


2 

P  .  (p,£>)dpd£> 

P  /  P 


V  >  0 


(1.5) 


As  before,  I  differentiate  under  the  integral  sign  and 
perform  the  remaining  integration  to  obtain 


(1.6) 


Making  the  final  trams formation  A  *  where 
s  *  10  log^e  »  4.34...  (this  notation  will  be  adhered  to 
in  the  remainder  of  the  thesis) ,  I  obtain  the  final  result 
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PA(A)  =  ?  Py(*> 


£  =  A/e 


—  00  <  \  <  oo 


(1.7) 


This  density  is  in  fact  of  the  same  form  as  Pr(<£)  (see 

<P 

Equation  A8  in  Appendix  A) .  The  second  moment  of  this 
density  is  infinite.  Physical  insight  of  this  phenomenon 
can  be  obtained  if  one  visualizes  the  random  walk  problem. 
The  amplitude  of  the  vector  is  p  and  its  phase  <j>.  When 
the  amplitude  of  p  goes  to  zero  as  in  a  deep  fade  it  is 
easily  seen  that  the  phase  can  undergo  very  rapid  changes. 
Likewise  because  A  *  e2(J/P  it  is  also  clear  that  A  can 
assume  very  large  values  when  p  is  small.  Thus  both  A 
and  i  are  governed  by  the  same  form  of  pdf  and  their 
variances  are  infinite.  It  is  also  noteworthy  that 
while  the  pdf's  for  5  and  x  depend  upon  both  and  v. 
Equation  (1.7)  depends  only  on  v. 

1.1.2  Joint  Densities  of  Amplitude  and  Amplitude  Rate 
an&  .Grossing  Rate  Statistics 

In  order  to  complete  the  family  of  joint  densities 
of  amplitude  and  amplitude  rate  for  a  single  source,  I 

have  derived  P  . <X/X)  and  P.  :(A,A).  This  analysis  also 

X  •  A  A ,  iV 

provides  an  alternative  method  for  deriving  P.  (x)  and 

X 

P^(A)  to  check  the  calculations  of  the  previous  section. 

In  addition,  I  derive  the  two-dimensional  characteristic 


functions  CL4]  which  will  be  of  significant  importance  when 

I  investigate  multiple  sources. 

For  x  =  P 2  we  have  p  =  x^2  and  P  =  ^  XX  •  In 

order  to  solve  for  P  .  (x,x)  I  make  a  two  dimensional 

transformation  on  P  - ( p , p )  (Equation  1.1)  as  follows: 

P  /  P 


3  2 


PX,X(X'*}  =  3x3x 


1/2  1  -1/2 

;  f7*x 


P„  . (p ,p) dpdp 

P  /P 


(1.3) 


By  differentiating  under  the  integrals  in  Equation  (1.8) 
I  obtain 


P,  *(x,x)  -  - — — : 

x'x  4a^3  v/2ttx 


[-  -*■ 


JL _ i 


exp  t - -  - — i  / 

2a,  2  q,"t  2"2 


8xpx  v' 


X  >  0,  and  -«<x<® 


(1.9) 


The  two  dimensional  characteristic  function  or  two 


dimensional  Fourier  transform  can  be  defined  as 
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M  (w,a)  = 
X^X2 


—  GO  —  00 


in»Xt  +  iaX2 

PXl,X2{Xi*X2)e  dXldx2  (1.10) 


For  px,x(x'x) 


M  •  (u,a)  = 

X  /  X 


P  ^  (X/X)  e1‘°x  +  iax  dxdx 
X  •  X 


(1.11) 


Performing  the  double  integration  yields 


M  -<u,a)  =  - - - 

1  -  i2a. 2u  +  4a.4v2a2 


(1.12) 


Recalling  that  y  =  lnX ,  we  have  x.  =  and  X  =  £sy. 


v 


Thus , 


•y  .Jey 


Py,y(y^>  ’  W 
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Pv  ;(X#X)dXdx 

X '  X 


(1.13) 


As  before  I  differentiate  under  the  integrals,  and  make  the 
final  transformations  A  =  sy  and  A  =  tv  to  obtain, 
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’a,A(A'A) 


exp  ( 
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exp 


-exp(^)  [— 

£  -I  _  2 
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2  o. 


8e  2c, 2v2 


— c0< A< 00 /  and  -«<A<a»  (1.14) 


The  two  dimensional  characteristic  function  of  ^(A,A) 
(Equation  1.14)  is 


MA,A(u,a)  *  2ev 


.  .  .(1.15) 


where  K  (x)  is  the  modified  Bessel  function  of  order  z. 

.  z 

Performing  the  integrals  over  x  and  A  in  Equations 
(1.9)  and  (1.14)  respectively  yields  the  marginal  densities 
of  x  and  A  which  are  given  by  Equations  (1.4)  and  (1.7) 
respectively  as  expected. 

Following  Rice  [10]  and  using  Equation  (1.14)  I 
derive  the  mean  axis  crossing  rate  for  the  amplitude  of 
the  signal  expressed  in  decibels. 


G(Ao)  =2 


•  • 


APArA(VA)dA  5  EtA|A  =  Ao] 


(1.16) 
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Thus, 


G(Aa) 


A«xp(^A0> 

■  _  4e2o. Jv/nr 
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-exp  (~)  [-i- 
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G(AC) 


Ivexp  (-*fo 
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exp(^n 


(1.17) 


If  we  transform  Equation  (1.17)  to  determine  the  axis 
crossing  rate  for  pQ  =  exp(A0/2e) 


G<po> 


^VpO 

a^/27 


exp 


(1.18) 


Equation  (1.18)  was  previously  obtained  by  Dyer  and 
Shepard  [15],  and  Hamblen  [9]. 

1.2  Multiple  Sources 

In  this  section  the  statistics  of  the  amplitude  and 
amplitude  rate  variables  are  analyzed  when  the  source 


structure  consists  of  many  independent  acoustic  generators 
distributed  at  various  frequencies  across  the  bandwidth 
of  the  receiver.  It  is  assumed  that  propagation  from  each 
source  is  phase  random  and  multipath.  An  initial  assump¬ 
tion  is  also  made  that  each  source  radiates  at  a  different 
frequency.  This  problem  is  applicable  to  the  case  of  deep 
ocean  ambient  noise  due  to  distant  shipping  [4] .  For 
this  case  the  spectrum  of  ship  radiation  is  assumed  to  be 
dominated  by  lines  disjoint  in  frequency,  so  that  each 
ship  contributes  as  many  independent  sources  as 
there  are  lines  within  the  observational  bandwidth. 

(Note,  the  observational  bandwidth  cannot  be  so  large  as 
to  include  harmonics.)  Henceforth  line  and  source  will 
be  used  interchangeably.  Furthermore,  the  analysis 
assumes  that  the  Fourier  components  of  the  received 
signal  over  the  observational  bandwidth  are  squared  and 
summed.  Thus,  by  Parseval's  theorem  it  is  the  square 

of  the  signal  amplitude,  in  this  case  x  *  P2  which  is 

L 

summed  for  each  source,  or  Xtota]_  *  xn  ^or  L  sources* 
As  noted  in  Reference  [4],  this  model  breaks  down  when 
sources  cannot  be  separated  in  frequency,  and  then  the 
model  must  be  modified  to  include  the  effect  of  two  or 
more  sources  which  may  be  radiating  at  the  same  frequency. 
This  coherent  problem  is  treated  in  Section  1.2.3. 


1.2.1  Amplitude  Densities 


Dyer  [4]  first  investigated  the  amplitude  statistics 
for  the  multicomponent  case  and  developed  the  framework 
for  the  problem  which  will  be  followed  here.  Three 
categories  cure  defined:  (a)  all  line  components,  L  in 
number,  arrive  with  the  same  long-time  average  intensity; 

(b)  all  line  components,  M  in  number,  arrive  with  different 
intensities;  and  (c)  N  groups  arrive,  each  with  equal 
intensity  components.  As  an  example,  Case  (c)  was 
applied  in  Reference  [4]  to  noise  as  might  be  measured 
at  low  frequencies  in  deep  water  near  Bermuda. 

As  pointed  out  in  Reference  [4] ,  and  as  will  be 
shown  in  Section  1.2. 1.3,  use  of  the  models  describing 
line  component  noise  often  entails  considerable 
computational  tedium.  This  complexity  often  motivates 
adoption  of  approximate  methods,  which  will  be  discussed 
in  Section  1.2. 1.4.  Cases  (a)  and  (b)  are  treated  in 
Reference  [4].  Exact  solutions  for  N =  2  and  N *  3  of 
Case  (c)  are  derived  in  Sections  1.2. 1.1  and  1.-2. 1.2 
respectively.  It  is  true  that  the  exact  solutions  of 
N  *  2  and  N  =  3  may  be  of  little  practical  value  as  most 
cases  of  interest  will  probably  contain  many  more  than 
three  groups.  The  analysis  is  performed,  however,  to 
build  the  theoretical  framework  of  this  problem,  perhaps 
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enabling  someone  else  to  solve  it  exactly  in  all 
generality,  and  more  importantly  in  relation  to  this 
thesis  research  these  exact  solutions  will  provide  a  check 
for  the  approximate  solution  which  is  presented  in 
Section  1.2. 1.4. 

1.2. 1.1  Exact  Solution  for  N»2,  Case  (c) 

Consider  L  line  components,  each  with  the  same 
long-time  average  mean  square  pressure  u.  The 
probability  of  the  sta  mean  square  x  is  P  A  ] 


P. 


'■LE 


(X) 


aLvL"1e"ax 

fLJS _  ® _ 

(L-1) i 


x  >  0 ,  a  >  0 


L  a  1,2,... 


(1.19) 


where  a  *  1/U/  the  mean  u  =  Lu  ,  and  variance 

XL  x 

a2  ■  Ly2  *  y2  /L.  Equation  (1.19)  is  the  Erlang  or 
XL  XL 

gamma  pdf  (see  Appendix  B) . 

Let  X]_  a^d  X2  136  distributed  according  to 

Equation  (1.19)  but  with  different  y^  and  different 

(i*l,2).  Further,  let  2X  ■  X]_  +  X2*  Then 

X 

P  (x)  *  /  Pi  (?)  P,(x-?)d?.  If  we  make  the  change  of 
2X  0  1  2 
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variables  Z  3 


(X) 


alLla2L2e”a2X  L2+L1‘1 

r  )  r  a2)  X 


j  tLl  (l-t)  2  etx(a2“al)dt 


.  .  .(1.20) 


The  integral  is  the  confluent  hypergeometric  function 
defined  in  its  integral  form  [16/17  ]  : 


M  (a, 6 , Z) 


n<5) 

r  (<5-a)  T  (a) 


1 

'  e2t  ta-l (x.t) 
0 


dt 


(1-21) 


This  is  a  well  tabulated  function  [IS  ,19,20]  alternatively 
defined  by  Kuramers  series  (absolutely  convergent) : 


M(a,<5,Z) 


(cl>n  Zn 


(1.22) 


where  (a)  *  o(a+l) (a+2) . . . (a+n-1) ,  (a)  =  1. 

n  o 

Applying  Equation  (1.21)  to  Equation  (1.20)  I  obtain 
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the  desired  result: 


a1Lla,L2XL2+Ll‘1e'xa2 

P  „  (X)  - - -  M  [L1,L.+L,,x(a2-a1)l  (1.23) 

2*  112  i  1 


with 


“Xj.  +  "x2  *  L1U1  +  L2“2 

0x\  +  *  LA  *  L2u2 


If  a1  »  a2  then  as  expected,  Equation  (1.23)  reduces  to 
Equation  (1.19)  with  L  »  +  L2,  upon  noting  from 

Equation  (1.22)  that  M(a,o,0)  *  1. 

To  obtain  the  density  of  the  level  in  dB  I  make  the 
now  familiar  transformation  A  =  e  lnx  in  Equation  (1.23)  : 


a  Lla  L2 

P  (A)  »  — i - - -  M  [L,  ,L.+L2,  (a,-a.)e(A/£)] 

2A  eTt^  +  Lj)  1  1  2  2  1 


expjj-  <L2  +  L1}  “  a2  e(A/£)] 


(1.24) 
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The  statistics  for  A  are: 
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where  M  »  +  L2»  and  y  =  Euler's  constant  ■  .5772...  . 

Note  that  if  a^  *  a2 ,  Equations  (1.25)  reduce  to 
Equations  (4)  of  Reference  [4]  as  expected. 

For  non-integer  values  of  and  L2 ,  Equations  (1.25) 
can  be  expressed  as  follows: 


uA/e  -  lnu2  + 


> 

al 

Li  «  (L. ) 

r  in 

a2"al 

N 

L  n 1 

n»0  n* 

a2 

y 

['P  (I^+L^+n)  ] 


(1.26) 


Li  -  (L. ) 

r  in 

n 

a2 
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[tp  (L2+L^+n)  - 


lna2l 2 


+  tp 1  (L^+L2+n) 


where 


ip(x)  is  the  Euler  Psi  function  [21] 


and  where 


Pl  (X>  *  35  <P  (X)  • 
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Equation  (1.23)  can  also  be  expressed  in  terms  of  the 
incomplete  gamma  function  [22,23]  which  is  a  special  case 
of  the  confluent  hypergeometric  function. 


1.2. 1.2  Exact  Solution  for  N»3,  Case  (c) 

I  now  consider  three  groups  with  (i*  1,2,3)  sources 
in  each  group.  Let  3;<  *  X1  +  X2  +  X3  *  2x  +  :<3  where  the 
X^  (i  *  1,2,3)  are  distributed  according  to  Equation  (1.19) 
(with  different  means  and  L^)  and  2X  is  distributed  • 
according  to  Equation  (1.23).  Then: 


X 

p  (x)  »  p  Y  (?)  PY  (x-Odc 
3x  2X  *3 

0 


As  before,  let  c  *  and 


(X) 


Li  L 2  Ll+L2+L3“l  — a^x 

al  a2  a3  x  e 


r(L1+L2)  T(L3) 


(1.27) 


and 
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I 


L2-l  L.+L--1  tx (a.-a-) 

(l"t)  t  e  •  M  [L1,L1+L2,tx (a^a^) ] dt 


.  .  .(1.28) 


Equation  (1.28)  is  evaluated  by  expanding  the  confluent 
hypergeometric  junction  in  its  series  form  (Equation  1.22), 
and  integrating  term  by  term.  When  I  combine  this 
result  with  Equation  (1.27)  I  obtain  the  result 


P  Y(X) 

3X 


L.  L.  L.  +L-5+L,-l  — a«Y  _  _ 

L±i  in >a 

n=0 


r(i(1+L2+L3) 


l2  al 
n! 


(L1^  n 

(L1+L2+L3)n 


•  M[L1+L2+n,  L1+L2+L3+n,  x  (a3  — a2 ) 1 


(1.29) 


with 


“jX  *  Llul  *  L2u2  *  L3u3 
0*X  *  LlUl  *  L2“2  +  L3u3 

Deriving  numerical  resu^^s  from  Equation  (1.29)  is 
straightforward  by  rewriting  the  summation 


dHUii 


In  these  I  define  m  =  M(a+n,  b+n,  z)  and,  a  =  L, +L- , 

n  i  z 

b  *  Li+L2+L3'  and  z  =*  x(a3-a2^* 

From  Equation  (1.29)  the  density  of  the  level  in  dB 

is  obtained  as  before: 

P  .(A)  -  i  eyP  v(ey)  (1.33) 

3a  e  3x  y  *  A/e 

The  statistics  of  the  level  are: 
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and  independent  random  variables 


(1.37) 


and  thus. 


(X) 


N 

n 

n=l 


‘s  +  a 


-)  n  esx  ds 


n 


(1.38) 


Though  Equation  (1.38)  is  attractively  compact,  it  is 
computationally  tedious  for  most  cases  of  interest,  i.e., 
large  N  and  large  Ln  as  you  might  expect  with  many  ships 
and  many  lines.  To  carry  the  analysis  a  bit  further,  I 
make  a  partial  fraction  expansion  of  Equation  (1.37), 


N 

(  n 

n*l 


n=l 
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l 

k=l 


'nk 


(s  +  an) 


a. 39) 


It  is  clear  from  this  expression  that 
sum  of  weighted  Erlang  distributions. 


P  (x)  is  a  linear 
Nx 

and  the  problem  is 


I  take  the 


reduced  to  finding  the  coefficients  C  ^ 
inverse  transform  of  Equation  (1.39)  by  applying  the 
calculus  of  residues  to  obtain 


j 


(1.40) 


and  finally. 


P  v(x) 
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=  (  n 

n=*l 


n, 


N 

l 

n*l 


Ln 

l 

k*l 


CnkXk-le'anX 

(k-i) ! 


(1.41) 


The  generally  intractable  nature  of  Equation  (1.40)  leads 
to  the  introduction  of  approximate  methods  discussed  in 
the  following  section. 


1.2. 1.4  Edgeworth's  Series  Approximation,  Case  (c) 
As  before, 

N 

NX  =  jJiXi 

where  the  are  independent  random  variables  distributed 
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according  to  Equation  (1.19)  and  in  general  with 

different  long  term  average  means  and  different 

number  of  line  components  L..  Although  P  (x)  cannot  be 

Nx 

solved  exactly  for  N  >  3  without  considerable  computational 
tedium,  it  can  be  approximated  in  the  main  lobe,  and  as 
we  shall  see,  quite  accurately  and  easily  by  an 
Edgeworth's  series  04]: 


P  (X)  -  tt- 
NX  > 


ZU>  -  jr  ysz<3)  (?)  +  Jr  vez(4)  (5) 


♦#*<««,  - 


fr  Y3YeZ<7)  151 


JV  —  Z<5)  (?) 

ax 


280  »  ( 9 ) 

-  -s-r-  ;  (?)  +  ...  terms  in  higher  order 

3  •  s 


moments 


.  .  .(1.42) 


where 


Z  ‘  (£)  81  n  derivative  of  the  Gaussian  pdf, 


yv  K3 

=  coefficient  of  skew  =  — =  — 


X 


(1.43) 


y  =  coefficient  of  excess  = 


u 

— *-3 

cr* 

X 


|4 


(1.44) 


n 

u  =  n  central  moment,  and 

A. 

K  =  nth  order  cumulant  or  semi- invariant . 
n 

The  point- by- point  error  in  the  approximation  of 
Equation  (1.42)  is  of  the  same  order  as  the  first  term 
neglected  [14]  . 

A  computationally  more  efficient  expression  of 
Equation  (1.42)  can  be  derived  by  actually  taking  the 
derivatives  of  Z(5)  as  indicated  with  respect  to  £  and 
collect  terms  in  powers  of  £.  I  have  used  only  the  first 
four  terms  in  Equation  (1.42)  and  as  will  be  demonstrated 
later,  this  will  be  entirely  adequate  for  practical 
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computations  . 
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TT~  2(5) 

Nx 


U  + 


3Ye  54  Ys  , 

TT  '  IT  ■  T  5 


'5i  \ 

.  2  .  YS  _  J 

'  y  5 y2 
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(1.45) 


It  should  be  noted  that  because  Equation  (1.45)  is  an 
approximation,  negative  values  may  be  obtained  for  some 
regions  in  the  tails.  In  fact,  the  Edgeworth's  series 
performs  best  in  the  main  lobe  of  the  density  and  wor.se 
in  the  tails.  The  mean  and  variance  is  respectively 
(evoking  the  properties  of  sums  of  independent  random 
variables  14] )  : 
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(1.46) 
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For  the  gamma  pdf  [14] : 

K  -  unLr(n)  (1.48) 

n 

Furthermore,  for  a  sum  of  independent  random  variables  [14]: 


n 


lKn 


+  2Xn  +  3Xn 


+  NKn 


(1.49) 


I  use  Equations  (1.46)  through  (1.49)  to  obtain  the 

coefficients  of  skew  and  excess  for  P  (x) : 

Nx 


N 

l  M-.uJ 

i-1  1  1 


(1.50) 


1  N 

4  l 

a  i-1 


6Vi 


(1.51) 


y  and  y  are  zero  for  the  Gaussian  pdf  and  attain  maximum 
values  of  2  and  6  respectively  for  the  exponential  pdf 
(which  is  the  case  of  N  -  1  and  L-  1)  . 

Before  examining  the  nature  of  y  and  y  more  closely, 
let  me  return  to  Equation  (1.42)  and  make  the  log 
transformation  to  obtain  the  Edgeworth's  approximation  for 


i 
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the  level  in  dB: 
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(1.52) 


or  alternatively, 
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(1.53) 


where 


ncr  ■  exp  (A/e)  -  y  . 

Nx  N* 

As  a  result  of  the  fact  that  the  transformation 
A  ■  10  log^gX  is  nonlinear,  the  statistics  of  A  cannot  be 
found  without  exact  knowledge  of  P  . (A) : 
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N  >  „  Nx 
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Nx 


(1.54) 


Using  Equation  (1.41)  in  the  previous  section  and  inter¬ 
changing  the  order  of  summation  and  integration  I  obtain: 


Evaluating  Equations  (1.55)  yields  [241: 


\ 

\ 


) 


N  _  N  n  C  . 
.I'm  r  r  nk 


y  A  -  (e  n  a”11)  l  l  -Sp  [ip(k)  -  In  a  ] 
Na  n-1  n  n-1  k=l  a_  n 


(1.56) 


V  L. 


y*A  +  w*A  *  (£2  11  ann)  ^  “X  {['Mk)  "  ln  an] 

N  N  n-1  n=l  k-1  a_  ^ 


2  +  y1  (k) 


•nr' 


t 


i-  .■Mr' 


i 


-62- 


Though  Equations  (1.56)  are  exact,  they  are  predicated  on 
the  knowledge  of  the  C^’s  given  by  Equation  (1.40), 
which  in  turn  motivates  an  approximate  solution. 
Unfortunately,  direct  integration  of  the  Edgeworth's 
series  is  not  possible  because  term  by  term  the  integrals 
diverge.  An  approximation  can  be  made,  however,  using 
the  first  term  only,  that  is,  to  integrate  the  log 
transformed  Gaussian.  Thus, 
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(1.57) 


The  integrals  in  Equations  (1.57)  are  not  straightforward 

but  can  be  evaluated  as  follows.  First,  the  square  in  the 

exponential  is  expanded,  and  the  constant  term  removed 

from  under  the  integral  sign.  Second,  the  change  of 
1/2 

variables  x  *  u  ' 


is  made: 


(1.58) 
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Third,  the  radical  is  removed  from  the  In  and  the 
exponential  of  /u  is  expanded  in  its  power  series. 
Fourth,  the  order  of  integration  and  summation  are 
interchanged  and  finally  I  have, 
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(1.59) 
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The  integral  in  Equations  (1.59)  can  be  evaluated  [24]  and 
the  final  result  is  obtained: 
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Although  the  Edgeworth's  series  will  provide  an  accurate 
approximation  to  the  density  we  seek,  the  best  approxima¬ 
tion  for  the  moments  of  A  are  those  of  the  log  transformed 
Gaussian  given  by  the  statistics  of  Equations  (1.60). 
Though  not  immediately  derivable  from  Equations  (1.60) 

when  a2  is  small.  Dyer  [4]  obtains,  by  taking  an  expansion 
X 

around  the  peak  of  the  density, 

V  =  !lnux  •  <0x/2llx>u 


(1.60b) 
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Returning  to  the  Edgeworth's  series,  and  to  gain  more 
insight  into  the  way  the  coefficients  of  skew  and  excess 


( y  and  y  )  behave,  depending  upon  the  number  of  groups  N 
and  number  of  line  components  in  each  group  ,  I  have 
plotted  ys  vs  L  in  Figure  1  for  N=1  and  N  _>  3 .  In  like 
manner,  Ye  is  plotted  in  Figure  2.  Each  was  constructed 
as  follows:  Each  group  N  has  the  same  number  of  line 
components  L,  and  the  long  term  average  mean  y^  of  each 
group  is  quantitized  in  3dB  steps,  i.e.: 


'll1"1 


relative  units 


L.  *  L 


i  =  1,2,3,. . .  ,  N 


For  the  case  N=l,  Equations  (1.50)  and  (1.51)  yield: 
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In  the  limit  as  Equations  (1.50)  and  (1.51) 

converge  to: 


Ys  *  [1.43]  L_1/2 
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groups  ( N>oJ  as  a  function 
of  line  components  in  each  grou 


- 1»7  ■ 


Fig. 2  The  coefficient  of  excess,  ye,  for  one  group  (  N= I ) 
and  three  or  more  groups  (N>3)  as  a  function 
.  of  L  ,  the  number  of  fine  components  in  each  group. 


I 
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Ye  =  [3. SOIL  ^ 

The  convergence  is  very  rapid.  As  evidence,  the  exact 
value  for  N  =  3  is  only  about  10%  greater  than  the  limiting 
value . 

From  Figures  1  and  2  it  can  be  seen  that  when  N  =  1 , 
the  coefficient  of  skew  and  excess  closely  approximate 
Gaussian  values  when  L>_6.  When  N_>3  they  closely 
approximate  Gaussian  values  for  L_>3.  More  justification 
for  this  conclusion  is  given  in  the  next  several 
paragraphs . 

For  further  illustration  of  significance  of  the  value 
of  y  and  y  I  consider  the  case  N  =  1  and  L=2.  For 

5  6 

simplicity,  I  let  u  =  l.  Then, 

Px (X)  -  Xe"X  (1.61) 

-  2 

a2  -  2 
X 

ys  =  n 

Ye  *  3 

Equation  (1.61)  is  plotted  with  the  Edgeworth's  series , 
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Equation  (1.45),  and  its  corresponding  Gaussian  in 
Figure  3.  Note  that  the  large  values  of  y  and  y&  would 
predict  the  Gaussian  to  be  a  poor  approximation,  and  so 
it  is.  The  Edgeworth's  series,  on  the  other  hand,  is 
quite  good  except  at  the  tails. 

Taking  the  log  transformation  of  Equation  (1.61)  for 
the  same  simple  case  (N  =  1  and  L=  2)  ,  I  have  the  density 
of  the  level  in  dB: 


VA)  *  e  exp 


{¥- 


exp  (A/e) 


(1.62) 


Equation  (1.62)  is  plotted  with  the  Edgeworth's  series  as 
given  by  Equation  (1.53)  in  Figure  4.  Again,  the  values 
for  y  and  y  suggest  the  log  transform  of  the  Gaussian 
to  be  a  poor  approximation,  and  indeed  it  is.  But  the 
Edgeworth's  series  for  the  density  of  the  level  is 
remarkably  close.  Thus,  I  conclude  than  when  y  and  y 

S  c 

are  large,  the  Gaussian  is  not  a  useful  approximation, 
but  the  Edgeworth  is,  especially  when  dealing  with  the  pdf 
of  the  level.  Also,  it  is  interesting  to  compare 
estimates  of  the  mean  level.  Equation  (1.60 a)  gives 

=  2.28dB  while  the  exact  value  given  by  Equation  (1.25)  , 


P  v(X) 


A,  LEVEL  OF  SHORT-TIME  AVERAGE  MEAN  SQUARE  PRESSURE  dB  (arbi  ref.) 


Fig. 4  Probability  density  of  the  level  for  the  case  plotted  in  Fig.  3. 
Shown  with  the  exact  density  are  its  transformed 
Edgeworth  approximation  and  its  transformed  Gaussian. 


1.83dB. 


setting  =  L2  =  a.^  =  a2  =  1,  is  = 

As  a  second  example,  I  consider  Case  (c)  of  Figure  2 
of  Reference  [4],  with  the  specification  N  =  2,  a^  =  0.3, 
a-  *  0.6,  and  L,  **  L-  *  2.  For  this  we  have  y  *  1.14 
and  Ye  =  2.04.  The  log  transformed  Edgeworth  approximation, 
Equation  (1.53),  is  plotted  with  the  exact  density  in 
Figure  5.  Here,  Equation  (1. 60a.)  gives  u  ^  *  9.16d3  while 
the  exact  value  given  by  Equation  (1.25)  is  =  9.39dB. 

The  log  transform  of  the  Gaussian  is  not  a  good  approxima¬ 
tion,  but  here  again  the  Edgeworth  and  the  exact  density 
are  for  all  practical  purposes  identical. 

For  a  final  example  I  examine  a  three  group  problem 
(N »  3)  using  the  results  of  Section  1.2. 1.2  to  compare 
with  the  Edgeworth's  approximation  and  the  Gaussian. 

In  this  example,  a^  *  1,  a2  =  2,  a^  =  8,  =  2,  L2  =  4, 

and  =  16.  Using  Equations  (1.49)  and  (1.50),  I  find 
Y  *  .86  and  y  =  1.28.  Equation  (1.53)  and  the  exact 

3  S 

density  of  the  level  as  given  by  Equation  (1.33)  are 
plotted  in  Figure  6  with  the  transformed  Gaussian.  Once 
again.  Equation  (1.53)  is  exact  for  all  practical  purposes 
while  the  Gaussian  assumption  will  result  in  some  error. 

Thus,  from  our  examples  for  y  and  y  not  too  large, 
i.e.,  <_  1.5  and  3  respectively,  the  first  four  terms  of 
the  log  transformed  Edgeworth  series  adequately  represents 


LO 

TRANSFORMED 

GAUSSIAN 


LOG 

TRANSFORMED 

EDGEWORTH 


A  LEVEL  OF  SHORT-TIME  AVERAGE  MEAN  SQUARE  PRESSURE, 

dB  (arb.  ref.) 

Fig,  5  Probability  density  of  the  level  for  the  case  N=2, 
L(=L2  =  2  ,  /x(  =  10/3,  and  =  5/3  ,  /s=  1.14  and 
ye  =  2.04  .  Shown  with  the  exact  density  are  its 
transformed  Edgeworth  approximation  and  its 
transformed  Gaussian. 


II 


A  LEVEL  OF  SHORT-TIME  AVERAGE  MEAN  SQUARE  PRESSURE, 

dB  (arb.  ref.) 


Fig.  6  Probability  density  of  the  level  for  the  case  N=3, 
L(  =  2,  L2  =  4  ,  L3  =  I6  ,  [i |  =  l  ,  =  1/2  ,  and 

fiy  1/8.  ys  =  .86  and  ye •  1.28  .  Shown  with  the 
exact  density  are  its  transformed  Edgeworth 
approximation  and  its  transformed  Gaussian. 


the  levels,  while  the  log  transformed  Gaussian  requires 
Y  %  y  ^  0  for  it  to  be  adequate. 
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I  now  apply  the  foregoing  analysis  to  the  N=  15 
problem  of  Reference  [4].  This  problem  involves 
estimates  of  ocean  traffic  in  the  North  Atlantic. 

Components  from  each  of  the  ships  are  grouped  with  3dB 
quantitization  in  N *  15  steps,  with  components  within  each 
group  numbering  as  high  as  =  57,  as  is  shown  in  the 
first  three  columns  of  Table  I.  With  the  use  of 
Equations  (1.50}  and  (1.51)  I  find  ys  *  1.025  and 
Ye  ■  1.743.  In  Table  I,  the  contributions  of  each  of  the 
groups  to  the  mean,  variance,  and  coefficients  of  skew 
and  excess  are  also  tabulated.  As  is  clear  from  the 
table,  groups  9-15  contribute  very  little  to  the  overall 
density,  since  the  variance,  skew,  and  excess  do  not 
change  (to  within  three  decimal  places)  beyond  N=  8. 
Reduction  to  an  eight  group  problem,  however,  is  not 
much  of  an  improvement  over  the  15  group  one.  But, 
Equations  (1.46),  (1.47),  (1.50),  and  (1.51)  are  very 
simple  for  any  N,  and  in  comparing  y  and  Y-  for  this 
case  with  those  of  our  previous  examples,  we  can  expect 
the  density  of  the  mean  square  pressure  to  be  closely 
approximated  by  Equation  (1.45).  Similarly,  Equation  (1.53) 
should  yield  the  density  of  the  level  in  dS,  which  for  all 
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practical  purposes  would  be  exact.  Equation  (1.45)  for 
this  case  is  plotted  along  with  the  Gaussian  in  Figure  7, 
and  Equation  (1.53)  is  plotted  with  the  transformed 
Gaussian  in  Figure  8.  The  Gaussian  pdf  is  seen  to  depart 
significantly  from  the  Edgeworth  pdf,  even  for  a  case 
having  a  very  large  number  of  line  components  such  as 
may  be  appropriate  to  an  actual  oceanic  situation. 

However,  with  reference  to  the  examples  above  and 
Figures  3-6,  we  are  justified  in  expecting  that  the 
Edgeworth  pdf  is  virtually  exact. 

As  noted  earlier,  an  analytical  expression  has  not 
yet  been  derived  for  the  statistics  of  the  leg  transform 
for  N  >  3.  However,  I  can  estimate  the  mean  from 
Figure  8  to  be  (supported  by  a  numerical  integration  of 
Equation  1.53) 

UA  ^  7.2dB 

The  Gaussian  assumption  used  in  Reference  [4]  led  to  a 
slightly  higher  value  for  the  mean  (%  7.5dB)  as  a  glance 
at  Figure  8  would  explain.  Equation  (1.60a)  gives 

%  7.51d3  affirming  the  result  obtained  by  Dyer  [4].  The 
standard  deviation  as  computed  in  Reference  [4]  of  1.2d3 
appears  reasonable,  again  by  inspection  of  Figure  S,  while 
Equation  (1 . 60a) yields  - 


1 . 40dB . 


PvIX) 
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Ftg.7  The  Gaussian  and  Edgeworth  approximations  of  the 
probability  density  of  the  short-time  average  mean 
square  pressure  of  the  noise  that  may  be  sensed 
in  deep  water  near  Bermuda  in  winter,  for  a  l/3-oct 
band  at  60  Hz  and  an  omni-directional  hydrophone. 
The  ships  have  been  grouped  in  3 d B  steps  in 
fifteen  groups,  with  overall  mean  ,  variance  ,  and 
coefficients  of  skew  and  excess  as  specified  in 
Table  I . 


/ 


s- 


A  LEVEL  OF  SHORT-TIME  AVERAGE  MEAN  SQUARE  PRESSURE, 

dB  (arb.  ref.) 


Fig. 8  The  transformed  Gaussian  and  the  transformed 


Edgeworth  approximation  of  the  probability  density 

of  the  level  for  the  case  plotted  in  Fig.  7.  ! 

\ 

t 
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We  may  conclude  that  when  the  number  of  groups  which 
add  significantly  to  the  overall  density  is  greater  than 
3,  then  use  of  the  Edgeworth  series  approximation  can  be 
quite  valuable  in  estimating  the  pdf  of  the  received 
signal.  While  it  may  be  tempting  to  use  a  Gaussian  pdf 
for  values  of  <  1  and  y£  <  2,  significant  differences 
in  the  pdf  for  the  mean  square  pressure  or  the  level 
should  dissuade  us  from  this  course.  It  is  true  that 
the  mean  and  standard  deviation  of  the  levels  are  less 
sensitive  to  the  differences  between  the  Edgeworth  and 
the  Gaussian  approximations,  but  the  Edgeworth  is  not 
much  more  difficult  to  use  and  is  thus  to  be  recommended. 

The  Edgeworth's  series  approximation  therefore 
provides  an  easily  implemented  method  of  (1)  estimating 
the  statistics  of  the  level  in  d3  for  even  the  most 
complicated  realization  of  Case  (c) ,  and  (2)  revealing 
the  extent  to  which  the  Gaussian  assumption  is  a  valid  one. 
This  is  particularly  valuable  because  one  cannot  merely 
assume  on  the  basis  that  N  is  large  that  the  Gaussian 
assumption  will  be  a  valid  one.  The  critical  factor  is 
the  amount  cf  energy  in  each  group.  If  one  or  two  groups 
contain  most  of  the  energy,  then  the  governing  density 
will  be  significantly  different  from  the  Gaussian  and,  in 
fact,  will  mere  closely  resemble  the  density  that  would  be 
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associated  with  the  most  energetic  group. 

If  the  tails  of  the  pdf  are  the  primary  regions  of 
interest,  then  exact  computer  solution  may  be  a  viable 
alternative  for  a  complicated  realization  of  Case  (c) . 

It  should  be  noted,  however,  that  more  terms  in  the 
Edgeworth's  series  can  be  taken  to  obtain  any  arbitrary 
accuracy  desired,  or  alternative  methods  such  as  the 
Chernoff  bound  or  "tilting"  the  density  can  be  applied 
[25] .  These  methods  will  not  be  discussed  in  this  work. 

1.2.2  Amplitude  Rate  Densities 

In  this  section  of  the  thesis  I  will  derive  the 
pdf's  of  p,  x,  and  A  for  the  multiple  source  case. 

These  rate  variables  are  dependent  variables  with  respect 
to  the  amplitude  variables  except  for  certain  special 
cases.  This  fact  introduces  a  great  deal  more 
complexity  than  has  been  encountered  up  to  now.  Solution 
for  the  joint  densities  of  amplitude  and  amplitude  rate 
are  generally  required  before  the  marginal  rate  densities 
themselves  can  be  found.  Solution  of  the  rate  density 
for  one  variable  does  not  lead  by  simple  transformation 
to  the  solution  for  the  rate  densities  for  the  other  two 
variables  as  was  the  case  for  the  amplitude. 

In  general,  the  rate  variables  depend  upon  both  c£, 
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a  measure  of  the  intensity,  and  v2 ,  the  single  path  mean 
square  phase  rate.  This  implies  that  there  are  a  plethora 
of  different  possible  combinations  of  multiple  sources 
with  different  or  same  and/or  different  or  same  v2. 

The  breakdown  into  the  same  cases  employed  in  Section  1.2.1 
will  be  followed  here,  though  their  definitions  must  be 
expanded  to  include  v2.  In  Section  1.2. 2.1  I  examine 
Case  (a) ,  the  case  of  multiple  sources  or  lines  in  which 
cr£  and  v2  are  the  same  for  all  source/receiver  pairs. 

This  would  apply,  for  example,  to  noise  which  is  flat 
across  the  passband  of  the  receiver  from  a  small 
geographical  area  or  sector.  In  Section  1.2. 2. 2  I 
investigate  various  special  cases  when  the  a£'s  and  the 
v2  may  be  different  for  each  source/receiver  pair,  with 
more  general  applicability,  including  Cases  (b)  and  (c) . 

1.2 .2.1  Multiple  Components  of  Equal  Intensity 
and  Equal  Single  Path  Phase  Rate  *” 

I  will  first  solve  for  P.  (x)  which  is  the  pdf  for 

Xt/e 

X  from  L  sources  or  components  of  equal  a?  and  v2 .  From 
L  L 

*total  ”  E,  V  u  is  clear  that  *total  *  L,  V  Thus 
n—  i  n= 1 

for  this,  and  only  this  variable,  we  can  employ  all  the 

very  nice  properties  of  sums  of  independent  random 

variables.  For  the  pdf  of  £ 

LtHi 
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P.-  (x) 

XLE 


—  CO  —00 


X 

PX{yl)PX(y2_yl) ' * *P^(X“yL)dy- * -dYL 

—  00 

.  .  .(1.63) 


or 


1  f 

2tt 


b2L 

(b2  +  u2)L 


dui 


(1.64) 


by  making  use  of  M.  (oj)  (Equation  A6)  ,  where  b  =  l/2a£v. 

X 

For  x  ?  °»  therefore. 


Pv  (x) 
*LE 


h2L 

1  d1"1 

1 

•X 

3 

•r| 

<D 

1 _ 

"  £  Tl  ± 

2ir 

|_(u)+ib) 

id  =  ib 

(1.65) 


Applying  contour  integration  around  the  lower  half  plane 
for  x<  0  yields,  as  expected,  the  identical  result  given 
by  Equation  (1.65)  because  x  is  symmetric.  Applying 
Equations  (1.63)  and  (1.65)  for  L  =  1-6,  I  obtain 

pv  ('x)  s  I  exp(-bjx|) 

X1E 
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?.  (x)  =  j  {b2ix|+b}  exp  ( -b  |  x  | ) 

X2E 


Pv  ('x)  =  rk  I  x  1  2  +  3b2  Ixl  +  3b}exp(-b|x|) 

X3E 

PX4E(x)  =  ¥6  I Xi 3  +  6b3 |x| 2  +  15b2 |x|  +  15b>exp(-b|x|) 
<X)  =  {b5|xr+10blf|x|3  +  45b3|xi2+105b2|'x| 


+  105b}  exp(-b|x|) 


Pv  (X)  “  5 hj  <i>6!x|S  +  15b5|xP  +  105bMx|3  +  420b3[xl 


+  945b2  |x|  +  945b}  exp  (-bp  | ) 


.  .  .  (1 . 66a-f ) 


Upon  inspection  of  Equations  (1.66)  the  general  form  of 
the  density  for  arbitrary  L  emerges  (which  will  be  proven 
below)  ; 


<x) 


- — z  e*b|ii  i  l 


(L-l)  ! 2 


(L+k-2)  i 


k*l  (k-1)  I (L-k)  i  2 


| x I L“kbL~k+1 


.  .  .(1.67) 
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Clearly,  y.  -  0,  and  in  general 
XLE 


1 

(L-l) 12L“1 


N 

I 


k«l 


(L+k-2)  >bL~k+1 
(k-1) l (L-k) >2k“1 


^L-k+n  e-bx 

o 


.  .  .(1.68) 


which,  for  n  odd  is  0,  and  for  n  even: 


Et* 


1 _ 

(L-l) !2L_1bn 


(L+k-2)  t  (L-k+n)  1 
(k-1) I (L-k) 


(1.69) 


Using  the  independence  property  H4I  I  get 


E[* 


a?  *  8La?v2 


(1.70) 


This  result  is  identical  .to  Equation  (1.69)  when  n=2. 


A  much  more  elegant  solution  for  P.  (ft)  will  now  be 

XLE 

developed.  The  approach  yields  the  complete  solution  for 


the  equal  a |  equal 


case  including  P.  (p)  ,  P'.  (A)  ,  and 

°LE  * lE 


all  the  joint  densities  of  amplitude  and  amplitude  rate. 
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I  begin  with  Equation  (1.12).  Because  xT^  =  1  X  / 

L  “  n-1  n 

Xtc  *  I  X„  /  the  X- ' 3  are  indeoendent  of 
LiZ*  n  n 

n—  l 

each  other,  and  the  y  ' s  are  likewise  indeoendent  of  each 

n 

other  the  joint  density  P  .  (x,x)  can  be  expressed  as 

X  /  XLE 

follows : 


X/X 


LE 


(X/X)  * 


•  \ 

i 

i 

00 

f 

k  j 

4  IT  2  j 

1 

•  00  - 

-iojX 


dwda 


[w  +  i  (^r  +  2a?  v2a2)  ] 


(1.71) 


i*  h 

which  is  just  the  inverse  transform  of  the  L  power  of 
Equation  (1.12).  I  use  the  calculus  of  residues  and 
perform  the  integration  over  u  first,  which  has  a  pole  of 
order  L: 


P  . 
X,XLE 


(X/X) 


*  y 

i 

2tri 

00 

t  j  L- 1  r—  t 

i2oij 

(L-l)  !4tt2  ; 

w  dojL-1  L  J 
•00 

e 


-iax 


a 


0} 


2c£v2a2 


(1.72) 


and 


iMuyiiSS 
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(x*x) 


CO  1 

■  j£L-  f  e'(W  + 


2a?v2a2) x  -  iaft 

da 


.  .  .(1.73) 


Performing  the  final  integration  I  obtain 
result,  the  joint  density  of  the  sta  mean 
X  and  its  rate  x  for  L  sources  with  equal 


the  desired 
square  pressure 
and  equal  v: 


X/X 


LE 


(X,X)  = 


vL-f 

(L-l) !2^+1aJL+1v/2r 


exp 


it? 


3a^vzx 


.  .  .(1.74) 


To  find  the  marginal  density 
in  Equation  (1.74): 


P.  (x)  I  integrate  over  x 
*LE 


p*  (*) 

*LE 


ao 


•00 


p  . 

X'XLE 


(x»x)<ix 


(1.75) 


to  obtain: 
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Pv  (*} 

*LE 


L  1 


f  s 

L 

r  i 

2  4 

1 

1 

.  j 

JL_ 

•  Y 

lx  1 

r  (L)  /2ira£~2— 

U2J 

Kl-I 
L  2 

.2alv- 

.  .  .(1.76) 


where  T (L)  is  the  Gamma  function  and  K^(Z)  is  the  modified 
Bessel  function  of  order  v.  If  we  make  use  of  the 
identity  [20] 


K  .  (Z) 

L‘I 


e-Z  j  _ (L-2+k)  I _ 

k=l  (k-1) ! (L-k) ! (2Z)k_1 


it  is  easily  proven  that  Equation  (1.76)  is  identical  to 

Equation  (1.67) .  When  x  is  integrated  out  in 

Equation  (1.74) ,  I  obtain  P  ( x)  which  confirms  the 

XlE 

result  first  obtained  by  Dyer  [4]. 

The  next  step  is  to  solve  for  the  joint  density  of 
p  and  p.  This  is  easily  accomplished  by  the  following  two 
dimensional  transformations  of  Equation  (1.74) : 


Pp^LE(p'^)  "  3p3p  J 


f2p£> 


X/X 


LE 


(X,x)dxd/ 


(1.77) 


o 
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which  yields: 


P 


(p,£>) 


2L-1 

P 

(L-l)  !  2Li-1a^L-iv/2*f 


(1.78) 


This  result  is  rather  remarkable.  One  can  see  immediately 

that  p  and  p  are  independent  as  they  are  for  the  single 

source  case  and  furthermore  that  p  is  independent  of  L, 

and  P.  (5) ,  the  density  for  L  equal  receptions  is  in  fact 
PLE 

identical  to  P.(£>),  the  density  for  the  single  source. 
Integrating  over  $  and  p  respectively  in  Equation  (1.73) , 

I  obtain: 


P 

PLE 


(P) 


o2L~l 

(L-l) !2L-1a2L  eXF 


(1.79) 


The  mean  and  variance  are  respectively 


/2t  c1  (2L-1)  i I 
PLE  2"  (L-l)  1 
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where  (2L-1)  !  !  =  1  •  3  •  5  •  •  •  (2L-1) 


2a2L  - 


U 


2 


PL£ 


and 


Px  (|5) 

PLE 


-  exp 

m oi~z- 


2  a2v2 


u 


0  and  a?  =  <?fv2 
PLE  1 


(1. 80) 


To  complete  the  statistics  for  this  case,  I  make  the 
final  transformation  A  =  e  In  x  and  A  *  ex/X  in 
Equation  (1.74) : 


Mt 


A, A 


(A, A) 


LE 


3  2 

3A3  A" 


A  A/e 
re 


Pv  ^  (X/X)dxdx 
X'XLE 


(1.81) 


to  obtain: 
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S7  (L) 


L-l 

I 

w=l 


,  s2(i) 


0 


Again  this  result  was  previously  obtained  by  Dyer  [4]. 
Integrating  over  A  in  Equation  (1.82),  I  obtain 


1  ( 2L-1) ! ! 
£  2L (L-l)  1 


(1  + 


l/2v 

A2  :  L-ri/2 
4e2  vz 


(1.84) 


The  moments  of  A  for  the  equal  reception  case  are 
interesting: 


0  ,  n  odd 

n  vn2n“L+!j /t  (n-1)  II  (2L-n-2)  11  n 


(L-l)  ! 


,  2  <  L 


(1.85) 


,  a>  L 


Note  that  the  single  source  case  yields  an  infinite  second 
moment  as  reported  in  Section  1.1,  however,  for  L  _>  2 
the  variance  always  exists. 
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1.2. 2. 2  Multiole  Comoonents  with  Different 


The  most  general  case  of  this  problem  is  Case  (c)  of 
Section  1.2.1,  i.e.,  there  are  many  groups  N,  and  within 
each  group  there  are  receptions  with  identical  a2  and 
v2.  Because  there  are  now  two  parameters  which  can  vary 
for  each  source,  this  implies  that  in  general  for 
multiple  sources  the  number  of  groups  will  be  larger  for 
the  rate  variables  than  the  amplitude  variables.  I 
define  and  v?  to  be  and  v2  for  the  jth  source 
receiver  pair.  Considering  the  most  general  case  and 
using  Equation  (1.12)  I  find  the  characteristic  function 


Unfortunately,  I  have  been  unable  to  evaluate  this  integral; 
thus  a  retreat  for  the  moment  from  the  most  general  case 
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is  in  order.  Much  more  progress  can  be  made  in  solving 
for  the  density  of  x  because  here  we  are  dealing  with  a 
sum  of  independent  random  variables,  and  it  is  not 
required  that  we  know  the  joint  density.  This  analysis 
will  ultimately  lead  to  a  very  simple  and  useful 
approximation  of  the  density  for  x  under  the  most 
general  case.  Upon  completion  of  this  analysis  I  will 
return  to  the  problem  of  the  joint  densities  and  the  pdf's 
of  $  and  A. 

1.2. 2. 2.1  Solutions  for  P^(v) 

First  I  will  consider  the  exact  solution  when  the 
product  is  different  for  each  source/receiver  pair, 

Case  (b)  ..  I  apply  Equation  (1.63)  for  L=2: 

b,  -bJxl  b_  -b2lxi 

?•  ( x )  *  -*r  e  1  *  -4-  e  2  (1.88) 

X2D  2  2 


where  b^'  =  l/2a^v^,-  and  the  *  denotes  convolution. 
Performing  the  convolution 
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Rewriting  Equation  (1.89)  in  a  slightly  different  way. 


P* 

X2D 


Cft> 


(1.90) 


b,  -b3|x| 

Now  P .  (x)  *  P.  (x)  *  -5-  e  .  Examination  of 

X3D  X2D  * 

Equations  (1.88  -  1.90)  reveals  that  successive 
convolutions  can  be  done  by  inspection: 


P*  <X> 

X3D 


b! 


b!~bi 


blb3 

b!“bi 


b3e 


“bi  lx 


ble 


-bjxn 


bl 


b2”bi 


b2b3 

b|-b| 


b^e 


-b,lx 


b2e 


-^ixn 


(1.91) 


Rewriting  this  result  so  that  the  method  of  inspection  can 
be  applied  again. 
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*  i 

^  I 


1 


m  ’ 


blb3 


(b|  -  b* )  (b|-  b|) 


b2  “b2lxl 

T  e 


bibi 


b*)  (b| 


and  thus. 


(1.92) 


P*  (X) 
X4D 


b!b5bS 


(b|  -  b*)  (b*  -  b£)  (bj  -  b|) 


bl  _bi I  X I 
T  e 


blb3b4 


(b*  -  b£)  (b^  -  b|)  (bj  -  b^) 


b2  “b2!xl 

~Y  ® 


bib2bl 


(bl-bj)(bl-bl)(bj.bl) 


b3  -b3lx| 
T  e 


blb2b3 


(bj  ~  hi5  <b4  "  b2}  (b4  ~  b3J 


b4  “b4  I  X  I 

T  6 


(1.93) 


6t^C  •  •  •  • 

Thus  it  can  be  seen  that  P.  (ft)  is  a  weighted  sum  of 

XLD 

the  individual  single  source  pdf's  of  x.  It  is  also 


i 


L 


J _  _ 
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evident  from  the  L»  2-4  cases  that  in  general 


P. 

XLD 


(X) 


L 

n 


b* 


n  (bl-  b?) 
m-1  m  i 
m^i 


(1.94) 


or  more  simply  in  terms  of  and 


XLD 


(X) 


aiivi 


L 

n 

m*l 

m^i 


'”u»l  - 


almvra* 


exp 


(1.95) 


Solving  for  the  moments. 


kl 


L 

I 

i-1 


si 


L 

n 

m«l 

m^i 


L 

n 

j£L 


b5 


(bm 

zn 


-  b|) 


i 

Tc+T 


k  odd 


k  even 


(1.96) 


Making  use  of  the  properties  of  sums  of  independent 


random  variables, 


L 

l 

i-1 


•ii"i 


(1.9 


Equation  (1.96)  with  k=2  confirms  this  result. 

An  exact  solution  for  P.  (x)  for  the  most  general 

XN 

case  considered  in  the  beginning  of  this  section  would  be 
an  N-fold  convolution  of  Equation  (1.67)  or  (1.76)  with 
itself.  An  exact  solution  for  N=2  will  now  be  presented. 


»  .*<*> 


2X 


l 


(L. +k-2) 1  L.-k+l 

1  b.  1 


(L1-1)!2L1  k*l  (k-1)  1  (LL  -  k)  !  2k”1  1 


(Lj+ni-2)  I  Lj—m+l 

b. 


(L,  -  1)  1 2L2  m*l  (m-1)  ML,  -  m)  1 2ra“ 


1  “2 


*bl>l  -b2|*-y|  .Vk  _,V" 


e 


yl 


i£-y|  dy  (i-9 


I  will  now  evaluate  the  integral.  I  consider  x>  0 

(we  know  the  result  for  x  <  0  must  be  identical) ,  y  >  0 ,  and 
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then  y  <  0  s 

For  5c  >  0  and  y  >  0 , 

let  y  *  j(t 


I^+Lj-k-m+l  r*  “b1xt  -b2x|l-t|  I^-k  ,L2~m 


e 


1-t  dt 


.  .  .(1.99) 


To  remove  the  absolute  values  this  integral  must  be  further 
subdivided  for  t<  1,  and  t>  1: 


k  -  (b.  -  b7)  j(t  L,-k  L,-m 
e  i  ^  t  x  (1-t)  *  dt 


r  L.+L-”k— m+1 

*  X 

-b2X  f 
e 

* 

0 

k 

b2X  f  -(b1+b2)it  L2-m 


+  e  j  e 
1 


t  *  (t-1)  dt) 


(1.100) 


The  first  integral  in  Equation  (1.100)  is  the  confluent 
hypergeometric  function  encountered  in  Section  1.2. 1.1. 


/ 
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The  second  integral  is  a  degenerate  hyper geometric 
function  known  as  Whittaker's  function  [18,24].  Performing 
a  similar  analysis  for  y<  0  also  yields  Whittaker's 
function.  Simplifying  to  as  much  an  extent  as  possible, 

I  obtain  for  the  final  result: 


1  L1  L2  (L1+k-2)!(L2+m-2)!b?b^|x|Ct'hS"1 

P  - -  Li+t  l  l  - - - - - — - r— 7 

2X  r (L  ) r (L2) 2^1  L2  it-1  m-i  (k-D  ! (L-k) ! (m-1)  ! (L2~m) ! TT*  l 


exp(-b2ljcl)  M[a,Of6t  0>2“  bL)  ixl  1 


«P[-  \  (bi.“b2)  l^lHO^+bj)  |xl  ]  2  ' 


l-a-5  ICbi  +  b2)|X|]  +  (-1)4"1  exp[-  \  (bj-b^  |*|  ] 


2  2 


t<Vb2>ix|]  ‘W1* 11 

2  ’  2 


(1.101) 
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where , 


a  ■  -  k  +  1 

6  ■  Lj  -  » +  1 

M(a,5,2)  is  the  confluent  hypergeoroetric  function, 

end  W  (z)  is  Whittaker's  function. 

V  9  u 

In  order  to  check  Equation  (1.101),  if  I  let 
*  Lj  ■  1#  then  I  should  obtain  Equation  (1.89).  For 

exp(-b2l*|)M(l,2,  (b2-b1)[x[] 

+  «xp[-  \  (bx-  b2)  Me  I  ]  l(bx+  b2)  |  x  1 3  ”1  wo,-ist(bi+  b2lxM 

* 

♦  exp[-  \  (bj-bi)  U|]  K^  +  bj)  Ull"1  W0^[(b1+b2)  U|]  ■ 

.  .  .(1.101a) 


P  *(*} 
2X 


bLb2 ! X I 


/ 
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I  now  make  use  of  the  following  identities  [18  ,24] 


M(l,2,2z)  *  sinh  z  , 
z 


Vx(z)  *  w,.^(z)  ' 


U#“X 


WQ>x(z)  »  /z Trr  Kx(|-)  ,  and 


K,  (w)  »  /ir/2z  e 
1 


-w 


where  K^(z)  is  the  modified  Bessel  function  of  order  v, 
I  obtain 


W0,-*S(z)  *  W0,^(z)  *  ex^“ 


and  applying  these  results  to  Equation  (1.101a) 


P  *(*) 
2X 


blb2 


S~-~~g  [exp(-b1|  %\)  -  exp(-b2|x|)] 


'2  "1 


+  E~y  b2  [•*PC-b1Ul)  +  exp(-b2U|)] 


>?■-  ‘j 

J 
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It  is  easily  shown  that  this  result  is  identical  to 

Equation  (1.89).  However,  though  Equation  (1.101)  may  be 

exact,  it  is  of  significant  complexity  and  thus  of 

limited  engineering  value.  For  arbitrary  N  it  is 

therefore  clear  that  P  . (x)  is  analytically  out  of  hand 

NX 

at  this  time.  The  complexity  of  the  result,  however, 
once  again  motivates  an  approximate  solution.  An 
Edgeworth's  approximation  turns  out  to  be  quite  simple 
and  very  useful  for  even  the  most  general  case.  Recall 
the  Edgeworth's  series  is  given  by  Equation  (1.42).  In 
the  present  case,  however,  because  x  is  symmetric  about 
the  origin,  the  coefficient  of  skew,  y  ,  is  identically  0. 
Thus  it  remains  only  for  me  to  solve  for  y  ,  the  coeffi- 

C 

cient  of  excess.  First  I  solve  for  X4,  the  fourth  order 
cumulant  or  semi -invariant  0-4]  of 

K4  -  E[x“]  -  3 (E t^2l ) 2  (1.102) 

Performing  the  required  calculations  on  P. (x)  given  by 

X 

Equation  (1.4)  I  obtain 

K4  -  192ajv"  (1.103) 

Recalling  Equation  (1.49)  and  considering  N  groups  of  L. 
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identical  sources  each,  the  most  general  case,  I 


N 


K  4  *  192  I  Viivi 
tr  i*i  1X11 


Prom  Equation  (1.70)  I  have  for  N  groups  of 


N 


Nx  i=l 


T  8L.a.\.v? 


i"  liwi 


Applying  Equations  (1.44),  (1.104),  and  (1.105) 
the  result 


N 


V.  -  3 


Y  L.  a?.v? 
i-1  1  11  1 

< X  VVj’ 


and  the  Edgeworth’s  series  is  then: 


**<*>  "o17{z('''  +  TT  V<4>  (n)} 
N  Nx 


or 


have : 

(1.104) 

sources, 

(1.105) 

,  I  obtain 

(1.106) 

(1.107) 
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P  «<*>  - 


i-  S(n) 


y  Y 

e  n2  4.  « 

T  n  +  4 7 


(1.108) 


where  n  =  x/c  .  and 


2 (n)  =  -i-  exp  (-  i  n2) 

/I? 


It  is  clear  from  Equation  (1.106)  that  as  N  becomes  large 

or  the  L.  become  large,  P  . (£)  approaches  a  Gaussian. 

1  Nx 

The  maximum  value  of  y  is  3  for  a  single  source/receiver 
pair  and  for  L  identical  ones  y&  goes  to  zero  as  1/L.  In 
Figure  9  I  have  plotted  the  exact  density  for  N=l,  L=1 
given  by  Equation  (1.4)  with  its  Edgeworth's  approximation 
given  by  Equation  (1.108).  This  is  clearly  the  worst 
case.  I  have  let  a£v  =  1  for  convenience.  In  order  to 
gain  some  understanding  of  the  behavior  of  y  for  this 
case,  in  Table  II  I  have  listed  y  and  P^,(0)  from  the 
exact  density,  and  its  Edgeworth's  approximation  for 
various  values  of  L  identical  sources.  The  percent  error 
is  also  tabulated.  From  Figure  9  it  is  evident  that  the 
error  shown  in  Table  II  will  be  the  maximum  error  of  the 


6- 


i 


i 


X  volts  rod/sec 

Fig.  9  Comparison  of  Edgeworth’s  approximation  for  X  when 
N  s  L  *  I  given  by  equ.  (1.108)  to  the  exact  density  for 
X  given  by  equ.  (1.4).  This  is  the  worst  case. 


107 


TABLE  II 


The  maximum  point  by  point  error  of  the  Edge¬ 
worth’s  approximation  to  P.  (x)  for  various 
values  of  L.  *LE 

Also  listed  is  the  value  of  y  ,  the  coefficient 
of  excess  for  each  L.  ~ 


Error 


3.00 

.1939 

.2500 

22 

1.50 

.1184 

.1250 

5 

1.00 

.0916 

.0938 

2 

0.75 

.0771 

.0781 

1 

0.60 

.0678 

.0684 

1 

0.50 

.0612 

.0615 

.  5 

0.30 

.0463 

.0464 

.  3 
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approximation  in  the  main  lobe.  From  the  table 
convergence  to  the  Edgeworth's  is  very  rapid  and  it 
appears  quite  reasonable  to  assume  that  the  Edgeworth's 
approximation  will  be  a  very  good  one  when  v  £  1.  It  is 

C 

true  that  the  Gaussian  assumption  will  also  be  reasonable 
for  ve  <  1,  however,  because  the  corrections  required  by 
Equation  (1.108)  are  trivial,  accuracy  need  not  be 
sacrificed  for  expediency. 


1.2. 2. 2. 2  Solutions  for  the  Joint  pdf's, 

Pg(6)  and  P^(A) 

In  this  section  I  will  first  solve  the  problem  of 
two  different  source/receiver  pairs,  i.e.,  the  and  v2 
are  not  the  same.  I  will  then  generalize  the  analysis 
and  solve  for  N=2  of  Case  (c) .  Finally,  I  will  indicate 
the  analysis  required  for  arbitrary  N. 

For  two  independent  pairs,  the  joint  density  of 
Xi'X2.'*2  and  ?iven  bY  taking  the  product  of 

Equation  (1.9)  with  itself  with  different  and  v2: 


*Y  Y  V  ^X1»X1»X2»X2^  * 

X1»X1»X2»X2  1122  2a2 


1*2  /2TT4Xia^2  /2tt4 x2aj2vf 


exp[- 


-]  exp  [- 


8viivi  8x2aLv! 


(1.109) 


wn 


Now,  the  pdf  we  are  after  is  the,  joint  density  of  x  and 


\  where  X  *  X^  +  X2  an<*  X  “  /(1  +  %2m  1  note  that  in 

Equation  (1.109)  an<*  &2  are  Gaussian.  Because 

the  density  of  their  sum  is  given  by  convolution,  the 
following  can  be  done  by  inspection: 


P 

X^»X2 


*  exp 


*2  *  AV 

“i0! 

+  ^X2312V2) 

*1 

N 

X 

*2 

-  2ali 

2C12 

2(4XloJiVj  +  *X2^2V^)_ 

(1.110) 


The  final  result  is  obtained  using  the  convolution  again: 


*  (X.X)  - 


x.x 


2D 


exp( - ,  ) 

2oz 


exp 


rxi  (dr  -  dr 

L  2®u  2ai: 


12  0 


4allCT12  /27T[4xlallvl  +  4(x-x1)aJ2vJ] 


■)  - 


U  *[\ahvl  +  (x-x^^viu 


2  .,2 1 


dx.  (1.111) 


Xt: 


I  now  let  Xi 
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p  *  Cx.'x) 

*’X2D 


X^exp(-x/2aJ2) 


0 


exp 


-t(- 


x(«u  -  »n> 

aoL.°12 


■)  - 


Sx°12v2  +  8xt(0UVl",12',2>-i 


2  »\2  _  «2  .,2> 


dt  (1.112) 


I  have  not  been  able  to  evaluate  the  remaining  integral  in 
Equation  (1.112),  however,  as  will  be  demonstrated  below, 
numerical  integration  is  very  simple.  Making  the 
transformation  to  p  space  I  obtain: 


o,P 


2D 


(P.ft) 


p3exp(-p2/2a£2) 

2a^2/H 


t  (a 


v?  - 


2T 


11  l 


exp 


n  p2(<j2iz  ‘  °2u) 

-t( — — — — ) 

2aiial2 


2a!2V2  +  2t(<J"V? 


'11J1 


r2  v2) 
12  l- 


dt  (1.113) 


Integrating  Equation  (1.113)  over  p  I  obtain  the  marginal 
density  of  jS  for  two  different  sources: 


~1T. - 
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PA  (ft) 
$2D 


a!iCTi2 


m 


0  [au  +  t(aU  -  °ll)] 


exp 


~ _ af _ I 

-  2ff^2  +  2*AA  -  <J12V2- 


dt 


(1.114) 


The  integral  in  Equation  (1.114)  can  be  evaluated  by 
expanding  the  exponential  in  its  power  series  and 
integrating  term  by  term.  After  significant  labor 
I  obtain: 


P6  (ft) 
P2D 


Iftkf.o? 


iri2 


AA  - 


IK  -  0j2v|) 


l  (-l)n(n+l) 
n-0 


— Ttl 


2o?, v?  -  2o?„vJ 


^“11  1 


12  2-1 


r(-n  -  f, 


2anvi 


-)  - 


rc-a-|.  j2 


20? „V? 


12  2 


.  .  .(1.115) 


where  T(a,z)  is  the  incomplete  gamma  function  [23/24].  The 

appearance  of  |ft|  does  not  mean  that  P.  (0)  ■  0  because 

P2D 

the  limit  as  ft  ■+  0  of  the  incomplete  gamma  function  is 


infinite.  If  the  incomplete  gamma  function  is  expanded  in 


its  power  series,  I  get: 


*  6  «> 

“2D 


Idalzkl  I  !k|!  f  (. 
■/5?  (<,uvi-°i2i)" 


°12  “  ail 


2  ”!Q 


^iivi 


2  _ 


_**  .  .2 

a!2V 


l)a(tr+l) 


Q-0 


I 

m*0 


Sill 


m 


a! (a-n-^) 


*2  ' 

«V> 


2n-2nrt-l 

a12V2  ~  G11  1 


2n-2o+l 


.  .  .(1.116) 

In  obtaining  numerical  results  from  Equation  (1.116)  I 
encountered  overflow  problems  in  computing  successive 
terms  of  the  series  before  an  accurate  result  could  be 
obtained.  In  order  to  avoid  this  problem.  Equation  (1.116) 
should  be  rewritten  as  follows: 


'2D 


tf)  - 


•'5F(0Uvl-°l2',2)1 


o. »v_  l  (-l)n(n+l) 
^a-0 


(o?,  -  aMo^v21 


12 


n  12  2 


ailVl’C12V2 


Szir 


a«0  m!  (a-n--) 
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-  a  .v.  I  (-l)n(tr+'l) 

■L  n«0 


1j212  -  -  (.d°  (  £_  » 

ailVl  *  C12V2  -  m“°  “!  2aUVl 


.  .  .(1.117) 


Values  for  P.  (jS)  were  computed  by  numerical  integration 
P2D 

of  Equation  (1.114)  and  from  Equation  (1.117).  Depending 
upon  the  value  of  p  desired,  the  numerical  integration 
of  Equation  (1.114)  to  three  place  accuracy  was  about 
40  times  faster  than  use  of  Equation  (1.117).  The  value 
for  £  -  0  is  obtained  from  Equation  (1.114)  analytically 
and  exactly  to  aid  in  these  comparisons.  Performing  the 
integration  and  simplifying  as  much  as  possible: 


(0)  - 


1  FilVl  -  CT12V2 


m  °u'i  -  0i2vi 


fu&gul  -  g?2V2> 

<°i2  -  g(l>2^  <ailvl  -  °nv2) 


g?2-jii 
^ 2  -  °Uvl> 


pll°i2Vi  -  2gUvj  *  gnJLV2  -  Jl2(g?2  ~  gU)2g12v2  £ 

lL2gIlVi  -  °ilCTl2V2  -  gU°12Vl  -  gil(gi2  -  gU)2gUvlP 

\  1 


fulfil 

kgM1. 

ah  -  ah 


.  .  .(1.118) 


The  mean  of  p  is  zero  and  the  variance  of  p  is  obtained 
by  interchanging  the  order  of  integration  in 
Equation  (1.114),  integrating  over  $  first,  then  t: 


.  .  .(1.119) 

Despite  the  complexity  of  Equations  (1.114)  and  (1.117) 

it  appears  as  though  a  Gaussian  approximation  is  an 

exceptionally  good  one.  For  values  of  a^,  v£,  a*2,  and 

\>2  taken  from  data,  which  will  be  discussed  in  more 

detail  in  Chapter  III,  the  values  of  P.  (p)  from 

P2D 

Equations  (1.114)  and  (1.117)  deviate  from  the  Gaussian, 

using  the  value  of  from  Equation  (1.119),  only  after 

P2D 

the  fourth  significant  figure.  A  more  formal 
calculation  with  equivalent  results  is  obtained  if  the 
coefficient  of  excess,  y e,  is  solved  for.  I  apply  the 
definition  of  y  from  Equation  (1.44  ) ,  perform  the 
two  integrals  required,  and  obtain  yQ  for  $  for  two 
different  sources. 
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[3°^ 


3cuvi +  6(ffLvi 


a2  \i2)a?  ]  -  3c^ 

12  2  P2D  P2D 


. (1.120) 


Applying  the  values  of  the  parameters  used  above  and 

Equation  (1.119)  in  Equation  (1.120)  I  obtain  for  this 

example  that  ■  2.4  x  10-7  which  certainly  warrants  the 

Gaussian  assumption!  Coupled  with  the  result  of 

Section  1.2. 2.1  in  which  it  is  proven  that  P,  (p)  remains 

15  LE 

identically  Gaussian  and  independent  of  the  number  of 
sources,  it  appears  quite  reasonable  to  state  that  even 
for  the  most  general  cases,  P^fJ)  will  for  all  practical 
purposes  be  Gaussian. 

Making  the  transformation  A  »  elnx  in  Equation  (1.112), 
I  obtain 


PAA  (A’A) 
a»a2D 


5  eA/c  1 

~  25T71  ’ 

8£‘°u9u  ^  o 


*  £(9uvi  - 


r  f.A/c<9i2  -  w  . 

L  l 


8e2<7^2v2  +  8£2t(a21v2  - 


9uv2>J 


.  .  .  (1.121) 
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Integrating  over  A  I  obtain 


PA  (A) 
A2D 


j 

■I 


I2ta;2v|  +  t(»;^  -  o^)i  dt 


+  '<3*7  -  4;,Hi",nu2  +  - 


11 


12 


ai2V2> ] 


.  .  .(1.122) 


Making  an  enlightened  change  of  variables  allows  the  final 
integration  to  be  performed/  and 


PA  (A) 
A2D 


A 


2v2  -  C  A2/e2  +  4v2 

-4c2xvJ  +  B  +  -±5 -  E[1  + - ? - ~] 


; 


(^2  +  ^V2)372 


C  -  2v2  A2/e2  +  4v2 

-4aJ2v2  +  B  +  — E[1  + - y— H 


.  .  .(1.123) 


where : 


A 
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12 


1 


a 


2 

11 
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2  ,allvl 
7  a12 


C12V2, 


gllVl  *  a12v2 

a12all 


4  (v2 


a12v2 


\  2 

V25^- 


4(°^  •  g»VV 

glla12 


*ffiivi  - 


_ 2  _4  _  _4  _ 2 

alla12  allff12 


.a  .  2 

-r  ■=■  (at 


llvl 


2  2  \ 

-  a12V2JF 


A2  i  r 

Te7  “  2  o 


(al!Vl  ‘ 


lla12(ai2 


a12v2} 


Though  arithmatically  messy  obtaining  the  exact  solution 
from  Equation  (1.123)  is  approximately  50  times  faster 
than  numerical  integration  of  Equation  (1.122)  which 
results  in  three  place  accuracy.  The  mean  of  A  is  of 
course  zero  and  integrating  Equation  (1.121)  first  over 


A  and  then  t,  I  obtain  the  variance. 
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2s2 


aUai2 


CTiiCT52v! 


°12allu2 


a!2ailV2  " 


°U  -  ° 


11 


12 


-  2ail^2  + 


21n 


11 


12 


.  .  .{1.124 

Unlike  the  single  source  case,  the  variance  is  finite 

though  the  next  nonzero  moment  E  [  A^]  is  infinite. 

The  rule  derived  for  equal  source/receiver  pairs  appears 

to  apply  to  unequal  ones  as  well.  Namely;  EtA^g]  exists 

only  for  n  even  and  n  <  2L,  where  L  is  the  number  of 

• 

sources.  An  example  of  P:  (A)  is  plotted  in  Figure  10 

iV2D 

using  values  for  c2^,  a12'  v2  taken  from  data.  The 

Gaussian  using  the  variance  given  by  Equation  (1.124) 
is  also  plotted.  The  Gaussian  assumption  is  clearly  not 
warranted  in  this  case,  and  in  fact  because  densities 
for  A  for  all  cases  are  related  to  the  Longue t-Higgins 
type  density,  the  Gaussian  assumption  will  be  a  poor  one 
in  general  even  when  the  densities  are  appropriately 
normalized  (see  Figure  1,  Reference  [9.]),  in  which  case 
the  variance  of  the  Gaussian  will  not  be  relevant  to  any 
of  the  source  parameters. 

I  will  now  generalize  the  approach  used  to  solve  the 
two  source/receiver  case  to  solve  the  K  *  2  case.  That  is 
two  groups  of  pairs,  with  pairs  each,  with  the  same 


l 


I 


Fig.  10  The  pdf  of  A  for  two  different  sources,  equ.  (1.123) 
is  plotted  with  the  Gaussian  using  the  variance 
given  by  equ.  (1 .124). 


I  start  by  considering  the  join 


<?^  and  v:  (i  =  1, 2)  . 
density  of  xL  Ef  XL  E'  XL  E#  and  *L  E  given  by  the 
product  of  Equation  (1.74)  with  itself  with  different  / 
a*  ^  and  v?  (i«l,2),  rewriting  it  slightly  to  make  the 
first  operation  obvious: 


CL,E'  xL,E'  xL~E' 


kL2E 


(  X1'X1»X  2*X  2  ^ 


V1  V1 

<1  x2 


(L1“l)  !  (L2-l)  !2Ll+L2/27r4x'1aJ1v2  ✓^4x2ffJ2v|  a ^ 


•  exp 


X1 

2all 


•  2 

Xo 


23u 


8ffiiuixi 


So12v2x2J 


(1.125) 


As  before,  I  seek  the  joint  density  of  X/  and  x  where 
X  *  Xl^j.  +  and  X  =  XL  E  +  XL  E»  A  double  convolution 

will  yield  the  desired  density  and  because  xT  o  and  x_  _ 
are  jointly  Gaussian  as  before  the  first  convolution 
can  be  done  by  inspection: 
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Pv  Y  UvXJc)  = 

XL12'XL2E  1  1 


Ll-1  L2-l 


(Ll)  r  (L2)  2Ll+L2a^1aJ22  /2tt  (4x]_a^1v|  +  4x2a£2v|) 


exo  - 


rll  2c12  2(4xlallvl  +  4X2^£2^2)- 


Setting  up  the  final  convolution  after  making  a  change  o; 
variables  leaves  the  following  result: 


P  v  f.  (X/X)  = 


VV!  X 

x  exp(- 

r(L,)r(LJ2Ll+L2+1  /m  a?h  o?i2 


1'  '2 


l  L,-l  L--1 

.  t  1  (1-t)  2 _ 

J0  +  tCa*p*— *77|T 


11  12 


X(a12  "all) 

’ exp  -t( - — — -ii) 

2allc12 


Sxal2V2  +  8xt(allvl  -  ai2v2)-l 


. (1.126) 


This  is  a  key  result  and  is  essentially  the  solution 
of  the  entire  N = 2  problem  because  this  joint  pdf 
contains  all  the  information  required  for  the  solution 
for  the  pdf's  of  x,  p.  A,  x,  0/  A.  As  a  quick  check, 
if  I  integrate  Equation  (1.126)  over  x  (which  can  be 
done  by  inspection),  I  obtain  Equation  (1.20)  and  hence 
all  the  amplitude  densities  for  N=2  follow  directly. 
Integrating  over  x  in  Equation  (1.126), 


P  *<*> 
2X 


r(L1)r(L2)2^  /v  cgl  ag2  8lM+1"2 


L.-l  L.-l 
t  1  (1-t)  * 


r—7T7T 


+  (,"+1)/2 [I5^+  '<5^-  1 


iiL 


,2 

'll 


a,‘,  +  t  (a*2  -  a?, ) 


'’ll' 


CaT7aTTvT  +  t(ailai2vl  "  a12allv2)  J 


12  11  2 


]2'tdt 


.  .  .(1.127) 


where  M 


The  remaining  integral  is  clearly  not  a  straight- 
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forward  one,  however,  we  have  already  solved  for 

P  • (x)  in  the  previous  section  given  by  Equation  (1.101) 
2X 

and  its  complexity  has  been  noted.  Clearly  the  x 

densities  were  handled  by  much  more  efficient  means  in 

the  previous  section.  If  I  transform  Equation  (1.126) 

into  p,i>  and  A,  A  space,  than  I  can  obtain  expressions  for 

the  densities  P  .(p)  and  P  ‘.(A)  and  their  moments  which 
2 p  2a 

have  been  unattainable  until  now. 


2L1+2L2-1 


,P,P 


<P,(5) 


exp(“  2C^") 
za12 


r(L1)ra2)2Ll+L2~1  a^122L 2 


V1 

t  1  a-t) 


V1 


o  ^2 +  t(ciivr-"  turv 


_ 

P2 

[  2°lla12 

2ouv2  +  2t(au',l  ‘ 

.  .  .(1.128) 


Integrating  over  p  I  obtain: 
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<L  +1.  -1)!  a^2c£l 

?  6»)  -  — — - ~rL- 

2P  ra  )ra2)  JS 


l  L-l  L-l 

.  f _ t  a-t)  _ 

J ,  +  t(a^  -  o-2vJ)  [a21  +  t 

.  .xS- - - lit 

l  20‘2v|  ♦  2t(^uvJ  -  o^|)J 


-  a2u)lLl^ 


(1.129) 


As  a  check,  if  I  let  L1 = L2 » 1  I  obtain  Equation  (1.114) 
as  expected.  Numerical  integration  of  Equation  (1.129)  is 
not  difficult.  Solving  for  the  variance,  I  obtain 


2L,  2Lx  L  -l  L--1 

3(L1  +  L2~1)!CT11  au  rl  t1  (1-t)  ^2V2  +  t(ailVl  ~  a!2V2: 

2/2  r(L1)r(L2)  1<J11  +  t(ai2  -  ffU)]Ll+L2 


[ail  +  t(a12  -  ail)] 


l1+l2 


.  .  .(1.130) 


To  obtain  the  final  rate  density  for  this  case,  I 


transform  to  log  variables: 


P2A,A<A’A> 


A  eA/E 

£2r(L1)ra2)2Ll+l2_1  a^1 


I 


L  -!  L  -1 
e  1  (1-t)  * 


•  ^al2V2  +  t(CrUVl  “ 


•  exp's 


« 

-  "U>1 

* 

A  2  A/t 

A  e 

L  2ouaL  - 

8eZai2V2  +  8e2c<anvi  "  J12V2}’ 

dt 


.  .  .(1.131) 


and  integrating  once  more  over  A,  I  obtain  the  desired 


result 


ra,  +  l,  +  4) 


V(A) 


er(L1)ra2)2Ll+L2+1  m  agi  og*  J0 


12^ 


2al2 


+  t 


2ailff12 


8e2a^  +  8e2t(oLvl  '  ci2v2> 


-LrL2-2 


di 


.  .  .(1.132) 

Equations  (1.131)  and  (1.132)  reduce  to  Equations  (1.121) 

and  (1.122)  for  L1=*L2  =  1.  The  integral  for  the  variance 

of  A2  is  not  straightforward  and  cannot  be  performed 

unless  and  L2  are  known.  This  completes  the  analysis 

for  the  rate  densities  for  Case  (c)  when  N=2. 

The  solutions  for  arbitrary  N  are  very  complicated. 

Fortunately,  an  Edgeworth's  series  approximation  has  been 

derived  for  P  .  (x)  (Equation  1.108)  and  it  has  been  shown 
Nx 

that  P  ,($)  will  be  accurately  approximated  by  a  Gaussian. 
Np 

Unfortunately,  P  ! (A)  cannot  be  successfully  approximated 

Na 

by  either  an  Edgeworth's  series  or  the  Gaussian.  The 
coefficients  of  skew  and  excess  cannot  be  found  for  A 
for  arbitrary  N,  thus  ruling  out  the  Edgeworth's  series, 
and  though  one  can  certainly  fit  a  Gaussian  to  A,  as 


/ 
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noted  earlier,  the  fit  will  not  be  good,  and  the  variance 

obtained  from  the  Gaussian  will  not  be  relevant  to  any 

physical  parameters  of  the  problem.  For  realizations 

with  N > 2  in  which  P  1(A)  is  needed,  and/or  approximations 

N 

for  P  .  (£>)  or  P  .  (£)  are  not  good  enough,  the  following 
Np  Nx 

procedure  can  be  used  for  arbitrary  N  provided  computer 
time  is  available. 

For  arbitrary  N,  following  the  approach  used  for 
N  ■  2,  I  can  write: 


^L1E,aL1Ej 


■Vv 


l'Al» 


,xN^N) 


L-l 

J  _ Xn 

n-1  r<Ln)2Ln  ,2ir4xna^' 


'n 


8a?  v2y 
In  nV 


(1.133) 


i 

As  before,  the  sum  of  the  *  '  s  involves  the  convolution  of 

n 


Gaussians  which  can  be  done  by  inspection  and 


_N. 


exp 


/8F7J1  xiciivI 
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•—  8  .2,  x  ■  o  ?  •  v 
1=1  *i  li  i 


L  -1 

N  *n 

U  _ 2 _ 

n-1  r(Ln)2Lnc^n 


exp 


'n 


2*L-i 


(1.134) 


Unfortunately,  this  is  as  far  as  the  analysis  can  be 

taken  without  considering  a  specific  N.  Thus,  the 

procedure  would  be  to  numerically  compute  N  convolutions 

with  respect  to  the  x„'s  to  obtain  P  *(xrX).  Then  one 

n  ^jX  t  X 

more  integration  over  x  yields  the  exact  density  for 

Xyj?  likewise,  integration  over  p  and  A  after  simple 

transformation  yields  the  exact  density  for  and  AN 

respectively.  Perhaps  an  easier  route  to  P  .(X/X) 

NX/X 

can  be  taken  by  numerically  performing  the  two 
integrations  required  by  Equation  (1.87)  . 
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In  summary,  the  solutions  for  the  rate  densities 
for  Case  (c)  with  N>  2  are  of  sufficient  analytical' 
complexity  that  approximate  methods  or  numerical 
techniques  must  be  employed.  The  densities  for  even  the 
most  complicated  realizations  of  Case  (c)  for  £>  and  x 
are  well  in  hand  with  the  Gaussian  and  Edgeworth’s 
series  approximations  respectively.  For  A,  however, 
complicated  multiple  source  configurations  must  be 
handled  numerically. 

1.2.3  Coherent  Source  Addition 

In  this  section  of  the  thesis,  I  investigate  the 
statistics  of  the  received  signal  when  two  or  more  sources 
are  radiating  at  the  same  frequency  or  so  near  in  frequency 
that  the  receiver  is  unable  to  resolve  them.  The 
application  of  this  analysis  is  therefore  dependent  upon 
parameters  of  the  receiver;  i.e.,  its  resolution  R,  which 
is  a  function  of  its  averaging  time  T  (R  ^  1/T) ,  and 
the  source  bandwidths,  specifically  R>Af,  and  R>B^,V^, 
where  Af  is  the  frequency  separation,  if  any,  of  the 
sources  and  is  the  source  bandwidth  of  the  ith  source. 
With  these  conditions,  the  procedures  used  in  previous 
sections  are  not  applicable,  and  coherent  addition 
of  the  signals  must  be  considered.  This  analysis  will 


i 
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model  random,  narrowband,  sea  noise  from  biologies, 
weather,  distant  shipping,  etc.  that  overlap  in  frequency 
in  a  small  band  of  interest. 

The  multipath  signal,  p(t),  under  the  assumptions 
of  phase  random  propagation  can  be  written  for  a  single 
source  as 


W 

p(t)  «  r  l  cos  (tut  +  0  )  (1.135) 

n-1  n 


where  r  is  the  single  path  amplitude,  M  is  the  number 
of  paths,  and  8n  is  the  single  path  phase,  which  is 
distributed  uniformly  between  0  and  2tt.  For  many 
sources  at  or  very  near  in  frequency  as  indicated 
above , 


p(t) 


l  r,  7  cos(ut  +  9_ ) 
i-1  1  n-1  n 


(1.136) 


where  L  is  the  number  of  coherent  sources.  If  this 
c 

analysis  is  applied  to  random  sea  noise,  then 
Equations  (1.135)  and  (1.136)  imply  that  each  radiator 
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is  CW  with  essentially  constant  mean  output  levels.  This 
is  perhaps  not  too  bad  an  assumption  if  integration 
times  and  record  lengths  are  not  very  long.  It  is 
assumed  that  AT  ^4,  V^,  such  that  phase  random  propagation 
is  obtained  for  each  source.  Thus,  Equation  (1.136) 
reveals  that  the  problem  is  one  of  solving  for  the 
statistics  of  a  random  vector  which  is  the  sum  of  L 

c 

random  vectors,  each  with  Rayleigh  distributed  magnitudes 
with  different  means,  and  uniformly  distributed  phases. 
Forming  the  envelope  of  Equation  (1.135),  I  obtain  the 
quadrature  components  for  one  source,  and  applying 
Equation  (1.136),  I  obtain  for  the  quadrature  components 
of  the  total  vector: 

L 

c 

X  *  pcos$  -  F  p  cos$ 
n-1  n  n 

and  (1.137) 

Lc 

Y  -  psin$  -  l  p  sinO 
n-1  n  n 

where  p  is  the  rms  amplitude  of  the  total  vector,  i  is  the 

total  multipath  phase,  pn  is  the  rms  amplitude  of  the  nth 

source,  and  P  the  multipath  phase  of  the  nth  source, 
n 


The  quadrature  components  for  the  single  source  are 
Gaussian  (8,9],  thus  from  Sqn  (1.137)  it  is  clear  that  X 
and  Y  are  Gaussian  as  well: 


P..(X)  -  — = -  exp 

x  /IWT*- 


Py(Y)  *  PX(Y) 


2a: 


(1.138) 


where 


Z  "2 
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and 


PX(X) 


exp 


X  z 
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IX 


'I 


P*(Y)  -  PX(Y) 


(1.139) 


where 


IX 


n»l 


a2  v 2 
In  n 


Note  that  the  short  hand  notation  for  Y  and  Y  is  meant  to 

show  that  the  form  of  the  densities  is  identical.  The 

quadrature  components  are  in  fact  independent.  Eqns  (1.138) 

and  (1.139)  illustrate  that  no  conditions  need  be  placed  on 

the  strengths  of  the  sources  or  relative  values  of  the 

v2  * s.  At  this  coint,  however,  I  should  caution  that  if 
n 

one  or  two  sources  are  propagating  energy  via  3  or  fewer 
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paths,  and  they  are  strong  relative  to  the  other  sources, 
dominant  path  affects  must  be  considered.  Please  refer 
to  Appendix  C,  and  for  dominant  paths  plus  diffuse 
(Rayleigh)  noise  see  Reference  (26] . 

Continuing  from  Equations  (1.138)  and  (1.139),  the 
analysis  is  straightforward.  The  resulting  densities  are 
in  fact  identical  in  form  to  the  densities  for  the  single 
narrowband  source  as  given  in  Appendix  A.  The  functional 
difference  is  that  wherever  a2  and  v2  appear  in  Equations 
(A1-A8 )  they  must  be  replaced  by  a2  and  cr^/a2  respectively. 


e.g.  , 


exp 
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p  >  0 


(1.140) 
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— ^r 
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0 


(1.141) 


The  most  important  difference,  however,  is  that  now  the 
densities  for  A  and  $  depend  upon  the  a£n's  while  for  a 
single  source  they  do  not.  It  should  be  noted,  however, 
that  if  v2  is  the  same  for  all  the  source/receiver  pairs. 
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whsn  ones  dedin  A  and  $  cljts  iricLspsncicnu  of  tins  diaplitude 

parameters,  c£n's,  and  are  in  fact  equivalent  to  the 
•  • 

densities  of  A  and  4>  for  a  single  source.  Even  though  all 
the  sources  are  at  the  same  or  near  in  frequency  the 
likelihood  that  v2  is  the  same  or  almost  the  same  may  not 
be  a  good  one  given  v2,s  range  dependence  [11]  (when  the 
internal  wave  model  applies,  see  Section  3.2). 

I  now  consider  a  narrowband  experiment  in  which  random, 
background  sea  noise  is  present.  From  the  analysis  above, 

I  can  write: 


”  "l  +  I,  °ln 
n^x 


and 


(1.142) 


r2  =  c?v“  +  T  a?  v2 
II  1  S  In  n 

n=l 


where  Lc  is  now  the  number  of  noise  ,  sources  in  the  analysis 

band  and  and  v2  are  the  signal  parameters,  and  the  cr2n 

and  are  the  noise  related  parameters.  I  now  make  the  . 

assumption  (perhaps  bad)  that  v2  is  a  constant  for  the 

signal  and  the  noise.  Thus,  =  a2v2,  and  I  define 

Lc 

a2  =*  a2  ,  which  is  one  half  the  sta  mean  square  noise, 
Nq  n*l  In  * 

and  therefore: 
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or 


CI  =  °1  +  aN 


all  -  ^  +  aN  ^ 


°i  -  al  {sk  +  X) 


0II  Clv2  (SNR  +  15 


(1.143) 


(1.144) 


where  SNR  = 


'N 


long  time  average  mean  square  signal 
long  time  average  mean  square  noise 


Applying  Equations  (1.144)  to  the  densities  and  statistics, 
e.g.,  Eqns.  (1.140)  and  (1.141),  the  effect  of  background  noise 
can  be  accounted  for.  It  is  clear  that  SNR  will  be  a  function  of 
the  analysis  bandwidth,  and  that  noise  levels  should  be 
measured  in  the  absence  of  the  signal.  By  narrowband  I  mean 
the  signal  bandwidth  B  should  be  <<2v  which  in  turn  should  be 
<<  1/T  where  T  is  the  observation  time  and  the  analysis  band¬ 


width  Ba  should  not  be  so  large  that  the  frequency  dependence 
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Finally/  if  one  is  analyzing  over  a  wide  band  and  the 
components  are  disjoint  but  each  is  a  sum  of  one  or  more 
sources  at  that  frequency,  then  all  the  results  of  the 
frequency  disjoint  analysis  can  be  applied  with  c|  replaced 
by  a2  and  v2  replaced  by  a^c2  where  these  parameters  are 
defined  as  given  in  Equations  (1.138)  and  (1.139). 

For  noise  which  is  continuous  over  the  analysis  band, 
the  statistics  will  be  a  function  of  the  total  bandwidth 


and  the  receiver  resolution.  All  the  analysis  of  coherent 


and  incoherent  (or  more  correctly,  frequency  disjoint  and 


independent)  sources  apply  except  now  the  number  of  groups 
N  is  given  by  B^/R  where  is  the  analysis  bandwidth  and 
R  is  the  receiver  resolution.  Because  the  sources  within 
each  group  are  coherent  c2  and  c2x  must  be  found  for  each 
group  using  Equations  (1.138)  and  (1.139).  The  groups  can 
now  be  considered  for  purposes  of  applying  the  results  of 
Sections  1.2.1  and  1.2.2  as  individual  source/receiver  pairs, 
with  the  a2,s  given  by  the  a2's  and  v2's  given  by  the 
aII//al's*  In  conc^us^on'  therefore,  the  statistics  of  the 
received  signal  are  a  function  of  receiver  parameters  as 


well  as  source  parameters.  For  large  observation  or  analysis 
time  T,  and  consequently  higher  receiver  resolution,  the 
analysis  of  Sections  1.2.1  and  1.2.2  will  most  probably  apply. 
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As  T  becomes  smaller,  however,  coherent  effects  must  be 
considered  as  developed  in  this  section.  Careful  regard  for 
the  inequalities  stated  in  the  beginning  of  this  section 
must  be  applied,  case  by  case,  to  the  data  set  of  interest 
and  the  receiver  characteristics,  in  determining  which 
analysis  applies. 


1.2.4  Crossing  Rate  Statistics 
In  previous  sections  I  have  solved  for  the  joint 
densities  of  amplitude,  and  amplitude  rate  for  many  multiple 
source  cases .  This 

allows  me,  therefore,  to  solve  for  the  theoretical 
amplitude  crossing  rates  for  these  cases  as  well. 

Following  Rice  [10]  and  Dyer  and  Shepard  [15] ,  the 
mean  crossing  rate  for  the  sta  rms  pressure  p  can  be 
defined  as 

00 

G(pq)  =  |  |p|Pp,. (po,5)d5  (1.146) 


where  pQ  is  the  axis  crossing  level.  For  L  equal 
source/receiver  pairs,  Equation  (1.78)  is  used  in 
Equation  (1.146)  and  I  obtain: 
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Equation  (9)  of  Reference  [15]  is  obtained  for  L*l. 
(Note  that  uQ  =  aj,  u2  =  s2v2.)  For  two  different 
sources/receptions  Equation  (1.146)  is  applied  to 
Equation  (1.113) : 


G2D(0o) 
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(a  +  6t)^  exp(-ct)dt 


(1.148) 


where  a  =  °£2V2  ' 


aiivi  **  aLy2  '  and 


2,1  1  ' 
P°  2aii  2ak 


If  I  make  the  change  of  variables  Z  =  a  +  St  the  integral  in 
Equation  (1.148)  can  be  evaluated: 
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where  y(v,u)  is  the  incomplete  gamma  function  [23,24]. 
Finally,  for  Case  (c)  N=2  I  apply  Equation  (1.146)  to 
Equation  (1.128)  and 


2L1+2L2-1  exp(_ 


^12J 


g2(pq) 


r(L1)r(L2)2Ll+L2“2  /2?  c2^  a2^2 


(l  L.-l  L--1  ,, 

t  (1-t)  (a  +  6t)'*  exp(-ct)dt 

0 


(1.150) 


where  a,  <5,  and  c  are  the  same  as  above. 
Integrating  I  obtain  [24]: 
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where  r (x)  is  the  gamma  function  [24],  and  (a, 8 , y;x,y)  is 
a  degenerate  hypergeometric  series  in  two  variables  [24] . 

For  coherent  sources,  phase  crossing  rates  as  well 
as  amplitude  crossing  rates  can  be  found.  As  with 
amplitude,  for  the  multipath  phase  [10,15] 


G(»0) 
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(1.152) 


Following  Dyer  and  Shepard  [15]  but  allowing  for  the 
effect  of  many  sources  at  the  same  frequency,  it  is  easily 
shown  that 
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Thus, 
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(Stt)  15  ^  exp (-  J^r) 


(1.153) 


(1.154) 


(1.155) 


where  a*  and  are  given  below  Equations  (1.138)  and 
(1.139)  respectively.  Allowing  for  the  normalization 
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employ  ed  in  Reference  [15]  ,  and  the  parameters  a2  and  a2^, 
the  form  of  Equations  (1.153  -  1.155)  are  identical  to 
Equations  (9)  and  (10)  in  Reference  [15] .  It  is  also 
noteworthy  that  Equation  (1.153)  remains 
unchanged  whether  the  v2,s  for  all  the  coherent  sources 
are  the  same  or  not.  In  terms  of  the  SNR  then 


(8tt)% 


2 

SNR  ,  Po  a  ,  Po  ,  SNR  ,  , 
1  +  SNR  cr1  expt~  2a1  1  +  SNR 


(1.156) 


I  now  consider  the  ratio  of  Equation  (1.155)  with 
Equation  (10)  of  Reference  [15]  (the  limit  of  Equation  1.156 
as  SNR  -*•  «)  , 


Thus  the  effect  of  noise  on  the  ratio  of  amplitude 
crossing  rates  to  phase  crossing  rates  is  a  multiplicative 
factor  given  by  Equation  (1.157)  which  depends  on  both  the 
SNR  and  the  axis  crossing  level  selected. 

Finally,  if  one  desires  the  crossing  rates  for  the 
mean  square  pressure  x  or  the  level  in  decibels  A,  the 
simple  substitutions  qq  =  and  pQ  *  exp(AQ/2s),  respectively 
is  all  that  is  required. 
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Because  of  the  many  cases,  the  notational  difficulties, 
the  complexity  of  the  equations,  and  the  plethora  of 
random  variables  involved,  I  will  attempt  a  summary  in 
the  form  of  a  table  of  the  multiple  source  section  to  aid 
the  reader  in  gaining  a  little  perspective  on  what  exactly 
has  been  accomplished.  In  Table  III,  I  have  compiled 
the  overall  results  of  the  analysis  for  the  various  pdf's 
of  interest  in  the  multicomponent  and  single  source  cases. 

On  the  left  is  the  breakdown  into  cases  based  on  source 
structure  and  signal  analysis.  Across  the  top  are  the 
random  variables  whose  pdf's  we  seek.  When  a  number 
appears  alone  in  a  box,  it  indicates  the  equation  number 
of  the  pdf  and  also  that  the  result  is  unique  to  this 
thesis.  Superscripts  (circled  numbers)  appear  in  boxes 
for  comments  below,  and  numbers  in  brackets  refer  the 
reader  to  those  references  in  which  the  equation  appeared 
previous  to  this  work. 

1. 3  Finite  Bandwidth  and/or  Modulated  Source 

In  many  actual  oceanic  situations  the  source  does  not 
exhibit  stability  in  frequency  but  in  fact  oscillates  or 
wobbles  about  a  center  frequency  which  can  be  characterized 
by  either  frequency  or  phase  modulation.  The  effect  of 


(L)Tf.i!*  tlir.iiily  can  U:  .tcciuutuly  approx  i  Stated  l>y  a  Ca'issi  on  tor  an  arbitral y  number  of  sources,  or  a  rou|»s. 
i  /)b  ini|>l«:  li  .ms for. ration  of  the  X  densities. 

I'UM.'ia.iiown  yi  i cal iy  until  the  inlujials  in  equs.  (l.8/)>r  (1 . 1  •4/can  be  evaluated.  Numerical  methods  must  ba  applied, 
at*  £in*jL<;  Eunice  with  oJ,  and  u*  /o*  icplaeiivj  u*  and  \>2  respectively. 


this  modulation  on  the  family  of  single  source  densities 
is  analyzed.  Even  "narrowband"  signals  have  finite  though 
small  bandwidths.  I  have  determined  a  criterion  for 
smallness,  as  well  as  the  bandwidth  effect  when  this 
criterion  is  not  met.  Furthermore,  the  analysis  presented 
in  this  section  reveals  a  method  by  which  finite 
bandwidths  and/or  source  induced  modulation  can  be 
determined  from  the  received  signal. 

1.3.1  Amplitude  and  Amplitude  Rate  Densities 

In  the  absence  of  modulation  we  can  write  the  multi- 
path  signal,  p(t),  under  the  assumptions  of  the  phase 
random  model  as  given  by  Equation  (1.135).  When  the 
source  is  frequency  or  phase  modulated,  we  can  write 
p(t)  as 

SI 

p(t)  *  r  l  cos[uit  -  M(t)  -  8  ]  (1.158) 

n=l  n 

where  M(t)  is  a  function  of  time  which  may  be  random  that 
represents  the  modulation.  As  indicated  in  Equation 
(1.158),  it  is  assumed  that  source  induced  modulation 
will  be  path  independent  which  implies  that  any  change  in 
signal  propagation  characteristics  (i.e.,  path  structure, 
or  volumetric  absorbtion)  will  be  independent  of  the 
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" instantaneous"  carrier.  This  is  apt  to  be  so  unless  the 
modulation  is  extreme.  This  assumption  also  applies  to 
the  finite  bandwidth  effect.  The  bandwidth  cannot  be  so 
large  that  different  propagation  characteristics  obtain 
for  the  extremities  of  the  signal.  Further,  it  is 
assumed  that  none  of  the  energy  in  the  signal  is  rejected 
because  the  "instantaneous"  carrier  is  outside  the 
bandwidth  of  the  receiver.  Likewise,  for  a  finite 
bandwidth  source,  it  is  assumed  that  the  entire  signal 
bandwidth  is  within  the  bandwidth  of  the  receiver.  Note 
that  depending  upon  the  specific  temporal  dependence  of 
M(t)  the  modulation  would  be  classified  as  either  frequency 
modulation  or  phase  modulation.  This  distinction, 
however,  does  not  alter  the  analysis  to  follow. 

I  perform  quadrature  demodulation  on  p(t)  and  obtain 
the  quadrature  components : 

N 

X  «  pcos[M(t)  +  $J  ■  r  l  cos[M(t)  +  0  ] 

n-1  n 

(1.159) 

N 

Y  »  psin[M(t)  +  $]  ■  r  [  sin[M(t)  +  9  ] 

n-1  n 

where  o  and  *>  are  the  amplitude  and  phase  respectively  of 
the  complex  envelope,  or  alternatively  in  the  terminology 
of  phase  random  acoustic  propagation  the  sta  rms  pressure 
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and  multipath  phase.  Taking  the  derivatives  with  respect 
to  time, 


M 

X  -  -r  l  [M ( t)  +  I  ]sin[M(t)  +  9  ] 
n-1  n  n 

(1.160) 

N  . 

Y  »  r  T  [M(t)  +  b  ]cos[M(t)  +  8  ] 
n-1  a  n 


In  terms  of  the  quadrature  components,  Equations  (1.159) 
and  their  derivatives.  Equations  (1.160),  I  can  write: 


P  -  (X2  +  Y2)^ 

(1.161) 

C  *  (XX  +  YY)  (X2  +  Y2) 


As  given  in  Section  1.1,  ;<,  x.  A,  and  A  can  be  expressed 
in  terms  of  p  and  p.  By  applying  Equations  (1.159)  and 
(1.160)  to  Equations  (1.161)  and  making  use  of 
trigonometric  identities,  it  is  easily  proven  that 
Equations  (1.161)  are  independent  of  M(t)  and,  in  fact, 
are  equal  to  the  result  obtained  when  M(t)  =M(t)  =•  0. 
Thus,  the  amplitude,  and  amplitude  rate  variables  are  in 
fact  independent  of  the  modulation.  It  also  follows, 
therefore,  that  the  joint  densities  of  amplitude  and 
amplitude  rate  are  independent  of  the  modulation. 
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1.3.2  Multipath  Phase  Rate  Densities 
For  the  multipath  phase  rate,  the  modulation  plays  a 
critical  role.  From  Equations  (1.159)  it  is  clear  that 
the  multipath  phase  with  modulation,  <J>„,  is  given  by: 

-  tan”1  |  =  M  ( t)  +  «f>  (1.162) 

and,  therefore: 


4>m  "*  M(t)  +  $ 


(1.163) 


This  result  can  also  be  obtained  from  the  single  path 
variables  alone  [the  extreme  right-hand  side  of 
Equations  (1.159)]  from  which. 


t  H  ,V 

|  l  cos[M(t)  +  9n]  l  [M ( t)  +§n]cos9n 


Af  M  x 

+  ^sinfM(t)  +  0nJ  J  [M(t)  +  9n]sin0nj  (1.164) 

Again,  using  trigonometric  identities,  I  find  that 
Equation  (1.164)  is  equivalent  to  Equation  (1.163),  in 
which  is  given  by  Equation  (1.164)  with  M(t)  =  M(t)  =  0. 
Note  that  from  Equations  (1.159),  Equation  (1.164)  is- 
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independent  of  r. 

M(t)  can  be  either  a  known  deterministic  function  of 
time  or  a  random  process  governed  by  a  pdf  P^(M)  .  In 
the  latter  case,  because  M(t)  is  independent  of  ocean 
parameters.  Equation  (1.163)  reveals  that  is  the  sum 
of  two  independent  random  variables.  Therefore  [14], 


P*M(b)  =  PA(M)*P-(*) 


(1.155; 


where  *  denotes  convolution.  The  mean  =  y* ,  and  as 

with  an  unmodulated  source,  the  variance  a*  =  ®. 

rM 

If  M(t)  is  a  known  deterministic  function  of  time, 
I  define: 


(1  + 


l/2v _ 

-  M(t)  ]  ^  3/ 


-  dt 


(1.156) 


where  T  is  the  length  of  the  time  series  (not  to  be 
confused  with  T,  the  averaging  time  of  the  receiver) . 

(b)  is  a  continuous  histogram  and  has  all  the  properties 
of  a  pdf,  i.e.,  it  is  always  positive  and  integrates  to 

one.  This  function  or  pseudo  pdf  can  be  employed  when 

• 

M ( t )  is  deterministic  but  not  periodic  for  a  given  ensemble 
of  time  series.  Equation  (1.166)  also  applies  for 
periodic  deterministic  M(t).  However,  for  periodic 


i 


deterministic  M(t)  in  which  there  is  exactly  one  or  n 
integral  number  of  periods  M(t)  can  be  randomized,  treated 
as  if  it  were  a  random  variable,  and  its  pdf  found 
enabling  use  of  Equation  (1.165).  If  there  are  many 
periods  in  the  record,  then  an  integral  number  is  not 
required;  however,  some  error  will  be  introduced.  As  will 
be  demonstrated  by  examples  below  and  in  Section  1.3.3 
for  periodic  modulation  functions,  Equations  (1.165)  and 
(1.166)  yield  identical  results.  For  many  interesting 
problems  in  the  ocean,  M(t)  may  be  deterministic  but 
unknown,  the  real  (nonrandom)  parameter  estimation 
problem  [25]  .  In  these  situations  one  will  obtain 
experimental  realizations  of  H*  (i)  from  which  one  is 
able  to  learn  characteristics  of  M(t) ,  as  will  be 
demonstrated  in  Chapter  3.  Also,  as  will  be  shown  in  the 
figures,  some  deterministic  modulation  functions  will  have 
easily  recognizable  histograms,  H*  ($) ,  and  in  fact 
knowledge  of  H*  (<$)  is  by  itself  a  valuable  piece  of 
information  to  have. 

I  -shall  now  consider  three  analytical  examples 
illustrating  the  effect  first  of  sinusoidal  phase 
modulation,  second  of  uniform  frequency  modulation,  and 
third  of  Gaussian  frequency  modulation  on  the  pdf  for  ?. 

For  sinusoidal  phase  modulation, 
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M ( t)  =  2sin(at  +  o  ) 

s 

and  thus,  (1.167) 

« 

M(t)  =  30 cos(at  +  <J>  ) 

s 

I  have  randomized  the  phase  with  <j>  (uniformly  distributed 
between  0  and  2r) ,  which  indicates  uncertainty  in  initial 
conditions.  I  obtain  the  pdf  for  M(t)  [3] , 

P*(M)  *  7  - “7 - T77  '  |M|  <  Sa  (1.168) 

^  17  (0202  -  M2)1/2 

Combining  Equations  (1.168)  and  (A8)  in  Equation  (1.165), 

I  find 


3a 

_ *1 _ 

iSa/Tv2  +  (x  -  y) 2 ) 3 (82a2  -  y2) 


(1.169) 


An  analytical  expression  has  not  been  found  for  this 
integral;  however,  numerical  integration  is  straightforward. 
Applying  Equation  (1.166)  with  T  =  2?rn/c  (n  is  any 

9 

integer)  and  M(t)  as  given  by  Equation  (1.167)  with  $  *  0 
(don't  forget  to  exploit  the  symmetry  of  the  cosine) 
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yields  Equation  (1.169)  as  well.  When  x=0  in 
Equation  (1.169),  the  integral  can  be  evaluated: 


(0) 


tt/v2  +  82a2 


_ So _ 

tt/2,  /v2  +  62a2~ 


(1.170) 


where  E[tt/2,  k]  is  the  complete  elliptic  integral  of  the 
second  kind. 

It  is  possible  to  make  some  progress  in  solving  for 
P:  (x)  if  I  make  use  of  the  convolution  property  of 
Fourier  transforms  or  characteristic  functions.  For 
independent  random  variables  [14] , 

M’  (oj)  =  M^(<j)M*(u)  (1.171) 

where  Mx(oj)  is  the  Fourier  transform  or  characteristic 
function  of  the  pdf  of  random  variable  x  [14].  M^(w)  is 

given  in  the  Appendix  under  Equation  (A8) ,  and  for  the 
sinusoidal  density  Equation  (1.168)  [24], 

M^j(co)  =  Jo(oj8  a) 

where  J  (z)  is  the  Bessel's  function  of  zeroth  order. 


/ 


Thus,  by  exploiting  symmetry  and  taking  the  inverse 


transform  of  Equation  (1.171), 


(x) 


IT 


uK.  (u>v)J  (uga)cos  ojx  dx 
1  o 

a 


(1.172) 


where  K^(z)  is  the  modified  Bessel  function  of  order  one. 
Expanding  the  cosine  and  integrating  term  by  term  [24] , 
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2n  +  3  2n  +  1 


(n  +  ±)  [  ( 2n  -1)M] 


8  2o2  ■ 


(1.173) 


where  (2n-l)l!  =  1*3*5  . ...  (2n-l);  (-1)11  =  1,  and 
F(a,b;c;z)  is  Gauss's  hypergeometric  series  [24]. 
Unfortunately,  Equation  (1.173)  converges  only  for 
Bo  <  v  and  x<v.  When  3 o=»v. 
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As  with  Equation  (1.173),  Equation  (1.174)  converges  only 
for  x<v.  P^($)  is  clotted  in  Figure  11  for  the  case  of 
no  modulation.  Equation  (A3) ,  and  for  various  values  of 
So  relative  to  v  using  Equation  (1.169).  Applying 
Equations  (1.170),  (1.173),  and  (1.174),  when  applicable, 
revealed  that  the  error  of  the  numerical  integration  of 
Equation  (1.169)  is  approximately  1%. 

For  uniform  frequency  modulation,  M(t)  =  Mt  and 
M(t)  = M,  where  M  is  a  uniformly  distributed  random 
variable: 


1/2A 


Pjj(M)  = 


M|  <  A 


otherwise 


(1.175) 


and  A  is  the  maximum  excursion  from  the  carrier  in  Hz, 
thus,  2A  is  the  bandwidth  of  the  modulation.  This 
characterization  of  5l(t)  would  apply,  for  example,  to  a 


deterministic  M(t)  such  as  a  saw  tooth.  As  before,  I 
apply  Equation  (1.165)  using  Equation  (1.175)  and  (A8) ,  to 
obtain: 


(x) 


x  +  A 

+  (x  +  A)  2 


x  -  A 

AT2  +  (x  -  A) 2 


(1.176) 


Equation  (1.176)  is  plotted  in  Figure  12  with  no  modulation 
(Equation  A8)  and  for  various  values  of  A  relative  to  v. 

For  the  final  example,  Gaussian  frequency  modulation, 


p.  (M)  =  - — 

tM  /2^2 


M2 


(1.177) 


G 


where  a2  is  the  variance  of  the  modulation.  As  before, 

I  convolve  Equation  (A3)  this  time  with  Equation  (1.177) 
and 


(x) 


■  _3 

(v2  +  y2)  ^  exp [- 


2ad 


(x  -  y) 2]dy 


.  .  .(1.178) 


As  with  sinusoidal  phase  modulation,  I  have  been  unable  to 
solve  the  convolution  integral  analytically.  However, 
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where  W  (z)  is  Whittaker's  function,  and  F(a,b;c;z)  is 

V,  u 

Gauss's  hyper geometric  series. 

Finally,  Equation  (1.176)  will  also  model  the  effect 

of  a  nonmodulated ,  but  finite  bandwidth  source  if  the 

conditions  stated  in  the  beginning  of  this  section  hold. 

It  is  clear  that  as  long  as  the  energy  is  uniformly 

distributed  on  the  average  between  f  -A  and  f  +  A. 

c  c 

Where  f  is  the  carrier  frequency,  then  Equation  (1.176) 
c 

applies,  and  the  bandwidth,  B,  is  given  by  B *  2A.  A 
glance  at  Figure  12  reveals  that  the  effects  of  the 
bandwidth  on  the  pdf  for  }  can  be  neglected  if  B<<  2v. 

The  above  analysis  also  reveals  that  modulation  coupled 
with  bandwidth  effects  are  additive.  Thus, 


Pi  (<*»)  38  p;  (<*>)*?:  ($)*p!($) 
*MB  *M  41 B  41 


(1.182) 


where  P*  (<j>)  is  the  pdf  of  the  bandwidth  which  we  have 
B  ^ 

assumed  is  uniform  between  ±B/2,  and  P*  (<$)  is  the  pdf 

for  $  when  a  narrowband  signal  is  modulated  and  the 

criterion  B  <<  2v  is  not  satisfied. 

Taking  the  three  examples  of  modulation  used  above, 

and  coupling  them  with  the  bandwidth  effect,  I  obtain  for 

sinusoidal  phase  modulation, 
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dy 


or. 


Pi  (x)  “  K-  (wv)  J  (oiSOcos  ux  sin  ~  du) 

♦mb  Bit  J  1  o  2 
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For  Gaussian  frequency  modulation. 


?:  (x)  =  - : 

9 MB  2B/2WI  J 


G  —oo 


exp  [-  jiy  (x  -  y)  2] 
G 


y  +  b/2 


y  -  B/2 


/v2T  (x  +  B/2)  2  /v2  +  (x  -  B/2)  2‘ 


My 


or 


—  f 

Btt  J 


-Wa2 


(uv)  cosujx  sin 


u>  B 


du 
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(1.184) 


(1.185) 


1.186) 


-160- 


and  for  uniform  frequency  modulation, 


Pi  (x)  »  7—  {/v2  +  (x  +  A  +  B/2)  2  -  /v2  +  (x  +  A  -  3/2)  2 

*MB  4AB  *- 

+  /v2  +  (x  -  A  -  B/2) 2_  -  /v2  +  (x  -  A  +  B/2) 2“J 

.  .  .(1.187) 

Finally,  it  should  be  noted  at  this  point  that  because 
a  finite  bandwidth  signal  is  indistinguishable  from  a 
"narrowband"  signal  which  is  experiencing  extremely  rapid 
uniform  frequency  modulation  as  indicated  above,  the  effect 
of  the  bandwidth,  as  with  modulation,  is  felt  only  by  $, 
the  phase  rate,  and  the  amplitude  variables  remain 
unaffected. 


1.3.3  Crossing  Rate  Statistics 

For  path  independent  source  induced  modulation,  I 
have  shown  that  the  amplitude,  amplitude  rate,  and  joint 
densities  of  amplitude  and  amplitude  rate  are  unaffected 
by  the  modulation.  It  follows,  therefore,  that  G(pQ) 
will  be  independent  of  the  modulation  and  will  be  equiva¬ 
lent  to  the  nonmodulated  single  source  result  given  by 
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Ecuation  (1.18)  . 

For  a  modulated  or  finite  bandwidth  signal,  G(?  ) 
is  as  expected  critically  dependent  on  the  parameters  of 
the  modulation  and  the  bandwidth.  To  find  G  ( )  it  is 
necessary  first  to  solve  for  £($,$)  and  this  in  turn 
is  crucially  dependent  upon  the  exact  nature  of  the 
modulation  or  bandwidth.  For  deterministic  modulation, 

I  consider  Equations  (1.162)  and  (1.153) .  The  joint  pdf 
for  4>  and  4>  without  modulation  is  given  by  [3,9]: 
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4tt  v 
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1 

^  3/2  ' 

— 1 

>  i  2  J 


0  < ■  <p  < 

Id  I  < « 


(1.188) 


Making  the  change  of  variables  given  by  Equations  (1.162) 
and  (1.163)  I  obtain 
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<P  r  <P 


M 


4irvT 


[1  + 


( 0  -  M  ( t )  ) 


3/2 


dt 


(1.189) 


For  the  interval  0  to  T  :  min[M(t)]  <  <S  <  2tt  +  max[M(t)] 

s 

and  is1  <  ®.  Note  again  that  the  H  function  has  all  the 
properties  of  a  pdf.  To  find  the  phase  crossing  rate 


/ 


G  (?  )  I  apply  Equation  (1.152)  and 


gm(0  = 


M  o  4ttvT. 


<t> 


0  -«  [1  + 


- r~T7i  dt  d* 

(<b  -  M(t)  )  J/ 
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•] 


Rewriting  Equation  (1.190),  I  obtain 


Gv(4  ) 
M  O 


4ttT 


[v2  +  M(t)  2  +  2M(t)<?  +  02] 
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Jo  [v2  + M(t) 2  -  2M(t) i +  i2] 3/2 


1  ' 

do  rdt 


Performing  the  integrations  over  $ : 


W 


1 

2ttTs 


/v2  +  M(t) 2  dt 

0 


This  is  as  far  as  one  can  proceed  without  the  exact 
of  M ( t )  .  (Note  that  if  M(t)  -*•  0  I  recover  the  no 


(1.190) 


(1.191) 


(1.192) 


form 


modulation  result.)  Equation  (1.192)  can  also  be  used  to 
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find  the  phase  crossing  rate  when  nonstationarities  make 
v  a  function  of  time  as  well.  Obviously,  however,  v(t) 
must  be  known.  If  M(t)  is  a  random  process  or  the 
modulation  is  periodic  and  exactly  one  or  more  integral 
number  of  periods  are  on  the  record  (if  there  are  many 
periods  the  number  need  not  be  integral) ,  then  the 
probabilistic  approach  used  in  the  previous  section  can 
be  applied  here  as  well.  From  Equations  (1.189),  (1.165), 
and  (1. 166)  , 

a-l93) 

For  sinusoidal  phase  modulation,  I  use  Equation  (1.169) 
in  Equations  (1.193)  and  (1.151)  to  obtain 


So 


W  =  fa 


ig0/F2a2  -  y2  /Tv2+(x-y):) 


dxdy 


.  (1.194) 


Integrating  first  over  x  I  obtain 


3  a 


W  =  7 


A JL 


/TPa2  -  y2  )  (v2  +  y 


2T 


(1.195) 
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The  remaining  integral  is  a  complete  elliptic  integral 
of  the  second  kind: 


W  *  w  /J2<j2  +  v*  E[I  '  rl 


(1.196) 


where 


/B272“2“ 

For  either  6  or  a  0,  Gx,  ( <f>  )  for  no  modulation  is 
obtained. 

To  demonstrate  the  equivalence  of  the  probabilistic 
approach  and  the  deterministic  approach  given  by 
Equation  (1.192),  I  consider  again  sinusoidal  phase 
modulation: 

M ( t)  =  Ssin(at  +  <J>  ) 

and 

M(t)  =  Bacos(at  +  <$>s) 

Assuming  I  have  exactly  n  cycles, 

2ni 

3 

GM(V  =  37  2n7  j  ^2~+Szo2coszGt  dt  (1.197) 

0 
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I  observe  that  the  integral  over  one  cycle  of  a  cos2x  is 
equal  to  four  times  the  integral  over  0  to  tt/2.  Applying 
this  observation  to  Equation  (1.197), 


w 


7r/2a 

1  4na  f 
2tt  2nir  j 

o 


/u2  +  32a2cos2ct  dt 


and  therefore 


G  (d  ) 


7r/2o 

p 


/v2  +  S2cr2cos2at  dt 


(1.198) 


I  make  the  change  of  variables  x  =  at  and  apply  the 
trigonometric  identity  cos2x  =  1  -  sin2x  in 
Equation  (1.198)  and  simplifying. 


W 


TT/2 

/I  -  [j32F27(v2  +  82a2)Jsin2x  dx 
o 


.  .  .(1.199) 


The  integral  in  Equation  (1.199)  is  the  definitive  form 
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of  the  complete  elliptic  integral  of  the  second  kind  and 
Equation  (1.196)  is  recovered  exactly. 

For  Gaussian  frequency  modulation,  Equation  (1.178) 
is  applied  and 


2  J.  „2\  3/2 


io,  +  r) 


x  exp[- 


(x  -  y) 2 1 dxdy 


.  .  .(1.200) 


Performing  the  integration  over  x  yields 
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where  D^fz)  is  the  parabolic  cylinder  function  [24].  I 
now  apply  the  identity  [24] 
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to  Equation  (1.193)  and  I  have 


2v*a, 


•*p(- 


W 


tt/27T  J  (v2+y2) 
o 


72  M(1'5'fe)dy 


(1.202) 


where  M(a,b,z)  is  the  confluent  hypergeometric  function. 
The  final  integral  in  Equation  (1.202)  must  be  performed 
numerically. 

For  uniform  frequency  modulation,  direct  use  of 
Equation  (1.176)  yields  a  value  for  G($  )  =  «.  This 
result  is  not  a  physical  characteristic,  but  rather  a 
consequence  only  of  the  mathematical  form  of 
Equation  (1.176).  The  integral  leading  to  Equation  (1.176) 


PI  (x)  -  tt  — - - 


4A  -a (v 2  +  Y)2>3/2 


(1.203) 


Applying  Equation  (1.203)  to  Equations  (1.193)  and  (1.152) 
yields 
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I  make  the  change  of  variables  t  =  x  -  y  and 


A  « 


a  (*  ,  =  Jd.  f  f  - -  dt  dy 

-a  — y 


Performing  the  integration  over  t  yields 
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|_/v2  +  y2  v2/v2+y2J 


dy 


and  finally  I  obtain  the  result 


W  -  £a  ln 


/v2  +  A2  -  A 


+  -±-  A72^a 

4tt 


(1.204) 


In  the  limit  as  A  -*■  0  Equation  (1.204)  converges  to  the 
no  modulation  result.  Note  if  I  let  A = B/2  in 
Equation  (1.204)  then  I  have  exactly  the  phase  crossing 
rate  for  a  finite  bandwidth  non-modulated  source. 

The  crossing  rates  for  modulated  signals  with  band¬ 
width  have  not  been  solved;  however,  the  procedure  is  quite 
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straightforward,  though  the  integrals  may  not  be. 

I  now  consider  the  ratio  of  amplitude  crossing  rates 
to  phase  crossing  rates.  Dyer  and  Shepard  [15]  obtained 
this  ratio  for  the  narrowband,  no  modulation,  single 
source: 


G<Pq) 

G(4»q) 


(8tt 


v  h 


(1.205) 


Equation  (1.205)  is  independent  of  v  and  depends  only  upon 
cr2  which,  being  a  measure  of  the  energy  in  the  signal, 
is  a  controllable  parameter  unrelated  to  oceanic  phenomena. 
However,  non- stationary  behavior  of  c£  due  to  ocean 
dynamics  will  affect  the  ratio  given  by  Equation  (1.205). 
Likewise,  if  the  source  is  modulated  or  has  a  bandwidth 
which  is  not  <<  2v  then  the  ratio  given  by  Equation  (1.205) 
will  be  affected.  For  these  cases,  the  ratio  G(pQ)/G(0o) 
will  be  a  function  of  c£,  v2,  and  parameters  of  the 
modulation.  Adopting  the  approach  in  Reference  [15],. 

I  obtain 


(  8  TT  )  J  ,  P  (O  ) 
1  0  o 


(1.206) 


where,  for  sinusoidal  phase  modulation. 
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for  uniform  frequency  modulation. 
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for  Gaussian  frequency  modulation. 
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“1 


(1.209) 


where  is  given  by  Equation  (1.202),  and  finally 

for  a  finite  bandwidth,  non-modulated  source  C  is  given 
by  Equation  (1.208)  with  A  =  B/2.  Except  for  the  finite 
bandwidth  result,  great  care  should  be  taken  in  applying 
the  formulas  in  this  section  to  insure  that  the  actual 
modulation  fits  the  kinds  of  modulation  assumed  here. 
However,  the  procedure  developed  in  this  section  can  be 
applied  on  a  case  by  case  basis  to  any  kind  of  frequency 
or  phase  modulation  to  determine  the  phase  crossing  rate 


factor. 
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CHAPTER  2 

COMPUTER  SIMULATION 

A  computer  simulation  of  phase  random  multipath 
propagation  was  originally  developed  for  two  reasons. 

First,  in  the  course  of  analysis  of  the  multiple  source 
cases  it  was  felt  that  a  computer  simulation  would  provide 
confirmation  of  the  rather  complicated  analysis  when 
data  were  unavailable  or  difficult  to  obtain  for  the  case 
in  question.  Second,  when  an  analytical  impasse  was 
reached  the  simulation  could  provide  insight  into  the 
nature  of  the  solution,  thus  aiding  in  the  analytical 
process.  The  simulation  fulfilled  these  two  objectives 
not  only  as  originally  intended  in  the  area  of  multiple 
sources  but  in  all  aspects  of  the  theoretical  development 
presented  in  this  thesis,  as  well  as  providing  confirmation 
of  some  of  the  basic  precepts  of  the  phase  random  model 
of  acoustic  propagation. 

2.1  Computer  Model  of  Phase  Random  Multipath  Propagation 
A  computer  simulation  which  generates  random  samples 
from  a  phase-random  multipath  process  has  been  developed. 
Currently  there  are  two  versions.  The  first,  called 
RANDPHASE,  simulates  up  to  50  equal  intensity,  equal  v: 
components,  including  sinusoidal  phase  modulation,  uniform 


i 
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frequency  modulation/  and  bandwidth  effects  on  phase  rate 
for  a  single  source.  The  second,  called  BURMRAN,  is 
capable  of  handling  an  arbitrary  number  of  groups  of 
unequal  intensity  but  equal  v2,  for  the  amplitude  and 
amplitude  rate  variables.  Both  programs  are  written  in 
FORTRAN  IV  and  were  run  on  an  Interdata  Model  80  computer 
with  an  IMLAC  display  processor.  The  FORTRAN  listings 
for  RANDPHASE  and  BURMRAN  are  contained  in  Reference  [27]. 

The  following  algorithms  are  applied  to  generate 
samples  for  the  amplitude  and  amplitude  rate  variables: 


v  —  y  r? 

'<L  i-l  1 


' 

~SI 

2 

r»  ~f 

y  cose 

n-1  ni_ 

+ 

LnIlSin9-'iJ  ■ 

*L  "  l  2ri 
U  i-l  1 


~N  M 

7  sin6  7  9  cos9 

n-1  ni  n-1  ni  ni 


W  H 

y  cos9  y  8  sin9 

n-1  ni  n-1  ni  ni 


PL  =  XL 


1/2 


1  v  -1/2 


PL  *  2  XL  XL 


-173- 


Al  =  10  log^Q  xL 

Al  *  [10  log10e]  xL/XL  (2.1a-f) 


where  r.  is  the  single  path  amplitude  for  the  ith  source, 

L  is  the  number  of  sources,  M  is  the  number  of  paths, 

9n  the  single  path  phase,  and  the  single  path  phase 

rate.  For  L=l,  the  single  source,  samples  of  the 

multipath  phase  <p  and  the  multipath  phase  rate  $  are 

generated.  Equation  (1.164)  is  used  to  generate  the 

samples  of  £.  The  input  parameters  to  the  program  are 

\)z ,  L,  N,  A  (or  equivalently,  B/2)  ,  3,  c,  and  I,  the 

number  of  samples  desired.  For  RANDPHASE  the  r^  are 

equal  and  are  set  to  one  for  convenience,  though 

initialization  to  some  other  value  is  straightforward.  • 

For  BURMRAN  the  array  r  must  be  specified. 

Each  program  uses  two  random  number  generators 

employing  a  machine  independent  congruence  technique  [28]. 

One  is  used  to  generate  uniformly  distributed  random 

numbers  between  0  and  2t  for  6  ,  and  the  other  is  used  to 

n 

generate  samples  of  It  is  assumed  that  6^  is 

uniformly  distributed  between  ±  / Tv,  thus  in  the  limit 

as  N  -*■  ® ,  E[9  1  =  0  and  E  [  9 2 1  =  v2.  Application  of  the 
n  n 
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program  revealed  that  if  the  number  of  paths  is  greater 

than  three,  phase  random  multipath  properties  are 

obtained.  Because  of  the  central  limit  assumptions 

inherent  in  the  phase  random  model  [9]  (e.g.  W  >  3)  ,  the 

exact  nature  of  the  density  for  is  unimportant. 

RANDPHASE  was  run  with  9  distributed  normally  without 

n 

any  change  in  the  results. 

The  number  of  independent  samples  of  6^  and  9^ 

required  for  any  given  simulation  is  given  by  the  product 

of  L,  SI,  and  I.  For  most  of  the  simulation  runs 

I  =  600  samples.  Where  L  and  SI  are  large,  the  number 

of  independent  samples  of  9^  and  can  become  quite  large. 

One  run  of  BURMRAN  simulating  the  SI  =  15  problem  of 

Reference  [4]  required  ^  400,000  independent  samples  of 

9  and  9  .  The  maximum  integer  that  can  be  accommodated 
n  n 

on  the  Interdata  80  is  32,767.  I  performed  a  run  test 
on  the  system  random  number  generator  (which  employs  the 
linear  congruence  method)  and  for  some  primes  the  maximum 
string  of  unrepeated  numbers  was  8192,  clearly  inadequate. 
In  order  to  increase  the  cycle  lengths  the  machine 
independent  congruence  method  is  used.  This  technique 
artificially  increases  the  maximum  allowable  integer 
number  with  factors  appropriately  segmented  so  that 
numerical  overflows  can  be  computed  without  machine 


overflows.  The  cycle  length  is  greater  than  107,  and  for 

g 

1.5  x  10  samples  generated  by  this  method,  the 
Chi-square  goodness  of  fit  test  [29]  for  uniformity  was 
passed  at  the  a  =  .05  level  of  significance,  though  only 
by  a  small  margin.  It  is  true  that  the  leftmost  digits 
are  the  most  random  and  that  for  large  numbers  of  samples 
some  correlation  occurs  as  more  of  the  rightmost  digits 
are  repeated.  Table  IV,  reproduced  from  Reference  [28], 
gives  the  repeat  characteristics  for  the  two  sequences 
or  generators  used.  Table  IV  indicates  the  randomness 
of  the  generator  as  a  function  of  cycle  length  or  number 
of  samples-.  As  more  and  more  samples  are  taken,  i.e., 
the  cycle  length  becomes  larger,  more  and  more  of  the 
rightmost  digits  are  repeated  which  can  result  in  a 
correlation  that  will  introduce  some  error  in  the 
overall  simulation  results.  Because  the  repeat  cycle 
length  is  so  large,  however,  for  most  simulation  runs 
this  error  is  small. 

When  modulation  or  finite  bandwidth  effects  are 
simulated,  the  Interdata  80  random  number  generator  is 
used  in  addition  to  the  two  described  above.  Because  L=1 
for  these  simulations  the  large  number  of  samples  for  9^ 

and  9  is  not  required. 

n 

When  the  programs  are  executed,  histograms  are 
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TABLE  IV 

Repeat  characteristics  of  the  random  number 
generater  for  the  two  sequences  used. 


Cycle 

Sequence  One 

Sequence  Two 

0 

735776465527g 

311037552421g 

81 

062221556427g 

113326416521g 

82 

650107434527s 

653503553421g 

33 

617512155527g 

514633562421g 

34 

65527g 

52421g 

85 

465527g 

552  421g 

a6 

646  5527  g- 

7552421g 
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produced  by  separating  the  range  of  samples  for  each  of 
the  variables  specified  by  Equations  (2.1)  and  (1.164) 
into  25  equally  spaced  bins.  In  RANDPHASE  each  histogram 
is  then  plotted  along  with  its  respective  theoretical 
density  and  the  Chi-square  goodness  of  fit  test  is 
applied.  The  area  under  the  theoretical  curve  is 
computed  numerically  using  the  trapezoidal  rule,  each 
of  the  twenty  five  bins  being  further  subdivided  into 
16  intervals  each.  The  expected  frequency  is  then  equal 
to  the  number  of  samples  times  the  area.  Bins  in  the 
tails  are  grouped  such  that  the  minimum  expected 
frequency  is  7.  Thus  the  number  of  independent  class 
intervals  (a  function  of  the  variable,  number  of  sources, 
paths,  etc.)  usually  varies  between  15  and  25  which  is 
close  to  the  criteria  specified  by  Bendat  and  i'iersol  [29] 
for  v  600  samples.  The  Chi-square  statistic  (;<2)  is 
computed  and  compared  to  the  pass/fail  value  (;<£:  .05), 
where  n  indicates  the  number  of  degrees  of  freedom  of 
the  statistic  and  .05  is  the  level  of  significance.  The 
hypothesis  is  accepted  if  x2  <  X^:  .05.  The  number  of 
degrees  of  freedom  n  is  obtained  by  taking  the  number  of 
class  intervals  and  subtracting  1  +  5  where  3  is  the  number 
of  independent  parameters  that  are  varied  to  fit  the  pdf 
to  the  histogram  (e.g.  3  =  0  for  the  uniform  pdf  and  3  = 
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for  the  Rayleigh,  Gaussian,  and  Longuet-Higgins  pdf's). 

It  should  be  noted  that  the  theoretical  densities 
were  not  least- squares  fit  to  the  histograms.  Given  the 
input  parameters,  we  know  a*  from  Equation  (1.2),  and 
v  is  obtained  by  taking  the  expected  value  of  6^  from 
the  computer  generated  3n's  (which  is  always  within  a 
few  percent  of  the  input  value  of  v) .  The  theoretical 
densities  are  plotted  and  the  samples  fall  where  they  may. 

In  addition  to  the  output  on  the  graphics  display 
a  printout  for  each  variable  for  each  run  of  the  simulation 
includes:  the  variable  simulated,  the  number  of  sources, 

paths,  samples,  number  of  expected  and  observed  frequencies 
for  each  bin',  including  an  overflow  bin,  the  observed 
and  theoretical  mean  and  variance,  the  range  of  values 
obtained,  and  the  actual  Chi-square  statistic  and  the 
pass/fail  value. 

2. 2  Simulation  Results  for  a  Single  Narrowband  Source 

Examples  of  the  results  of  the  computer  simulation 
RANDPHASE  for  a  single,  narrowband,  non-modulated  source 
are  given  in  Figures  13-36.  The  figures  are  plotted  by 
variable  for  N  =  4,  5,  and  12  paths.  For  all  these  runs, 
the  input  value  of  v  is  .007  rad/sec  and  I  =  700  samples. 
Table  V  compiles  the  Chi-square  statistics  ( x 2 )  for  each  of 
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Pig.  13  Computer  generated  histoqrom  O'  X,  Pig.  16  Computer  generated  histogram  of  X 

thesta  mean  sauore  pressure  for  N«4  *ith  ^  =  4  paths  and  equ.  ( 1.4). 

paths  compared  with  eau.  (A2). 


X  for  N  *5  paths  and  equ.  (A2).  for  N*5  paths  and  equ.  It. 4), 


0.09  r 


Pig.  IS  Computer  generated  histogram  of 
X  for  N  =12  poths  and  equ.  (A2). 


Pig.  18  Computer  generated  histogrcm  of  X 
for  N  *12  poths  and  equ.  U.4). 


Pa(A>  Pa(A> 
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Fig.  25  Computer  generated  histogram  of  A,  Fig.  23  Computer  generated  histogram  of  A, 

the  !*vel  m  d8's  of  the  s  t a  mean  square  with  N»4  paths  ana  equ.  ( 1.7). 

pressure  for  N*4  paths  compared  with  equ.  (A3). 


Fig.  26  Computer  generated  histogram  of  A  Fig.  29  Computer  generated  histogram  of  A 

with  N  *5  paths  and  equ.  (A3).  with  N*5  paths  and  equ  (1.7). 


Fiq.  27  Computer  generated  historgram  of  A  Fig.  30  Computer  generated  histogram  of  A 

withN*l2  paths  and  eau.  ( A3).  with  N  *  12  paths  end  equ.  (1.7). 


I 


5  O'— 

i 

' ±  401- 


a.  0. 1 0 


9  rad 

rig.  31  Computer  generoted  histogram  of 
the  multipot  h  phase,  for  N>4  paths  compared 
with  equ.  (A4). 


9  rad /sec 

Fig.  34  Computer  generated  histogram  of 
with  N»4  paths  and  equ.  (A8). 


9  rad 

Fig.  32  Computer  generated  histogram  of 
with  N  *5  paths  and  equ.  (A4). 


9  rad /sec 

Fig.  35  Computer  generated  histogram  of  ♦ 
with  N*5  paths  and  equ.  (A8), 


H  n  H 


Fig.  33  Computer  generated  histogram  of  ♦, 
with  N  •  12  paths  ond  equ.  (A4\. 


9  rad /sec 

Fig.  36  Computer  generated  histogromof  $ 
with  N*I2  paths  and  equ.  (A8). 
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Results  of  the  hi-square  test  for  the  RANDPHASE 
simulation  of  a  sincrle  narrowband,  non-modulated 
source.  For  x2  <  x«:  -05  the  theoretical  pdf 
passes  the  test  at  the  .05  level  of  significance. 
The  units  are  the  same  as  in  the  figures. 


Fig. 

Variable  x2 

X2:  .05 
n 

H 

Mean 

Std. 

Dev. 

Theory 

Data 

Theory 

Data 

13 

X 

19.36 

25.00 

4 

4.00 

3.64 

4.00 

3.24 

14 

X 

18.59 

25.00 

5 

5.00 

4.85 

5.00 

4.31 

15 

X 

8.62 

25.00 

12 

12.00 

11.19 

12.00 

11.34 

16 

X 

27.44 

23.68 

4 

0.00 

0 .00 

0.04 

0.03 

17 

• 

X 

19.93 

23.63 

5 

0.00 

0.00 

0.50 

0.46 

13 

X 

7.84 

23.68 

12 

0.00 

0.00 

0.12 

0.11 

19 

p 

24.02 

28.87 

4 

1.77 

1.71 

0.93 

0.84 

20 

p 

27.68 

28.87 

5 

1.98 

1.98 

1.04 

0.97 

21 

p 

12.85 

23.  87 

12 

3.07 

2.96 

1.  60 

1.56 

22 

• 

p 

14.37 

18.31 

4 

0.00 

0.00 

0.01 

0.01 

23 

• 

p 

9.79 

23.68 

5 

0.00 

0.00 

0.01 

0.01 

24 

• 

p 

6.50 

26.30 

12 

0.00 

0.00 

0.02 

0.02 

25 

A 

21.29 

22.36 

4 

3.51 

3.30 

5.57 

5.63 

26 

A 

11.24 

22.36 

5 

4.48 

4.55 

5.57 

5.43 

27 

A 

9.00 

22.36 

12 

8.28 

7.99 

5.57 

5.46 

23 

A 

10.61 

26.30 

4 

0.00 

0.01 

00 

0.20 

29 

A 

15.36 

26.30 

5 

0.00 

0.00 

JO 

0.11 

30 

A 

25.96 

26.30 

12 

0.00 

-0.01 

00 

0.11 

31 

P 

20.64 

36.42 

4 

7T 

3.16 

1.  81 

1.76 

32 

P 

37.43 

36.42 

5 

IT 

3.21 

1.  81 

1.78 

33 

P 

13.07 

36.42 

12 

77 

3.20 

1.  81 

1.  80 

34 

P 

23.64 

26.30 

4 

0.00 

0.00 

oo 

0.01 

35 

P 

24.  60 

26.30 

5 

0.00 

0.00 

00 

0.01 

36 

P 

9.36 

26.30 

12 

0.00 

0.00 

00 

0.02 
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the  figures  and  the  x^:  -05  in  which  n  is  the  number  of 
degrees  of  freedom,  which  is  equal  tc  the  number  of 
class  intervals  minus  one  (recall  that  the  densities 
were  not  fit),  and  .05  is  the  level  of  significance. 

For  x2  <  X^;  .05  the  theoretical  density  passes  the 
Chi  square  test.  Also  shown  in  Table  V  is  the 
theoretical  value  of  the  mean  and  standard  deviation  and 
the  measured  values  from  the  computer  generated  "data". 

First,  note  that  all  the  theoretical  densities  match 
the  histograms  very  well  as  indicated  by  results  of  the 
Chi-square  test.  Particularly  note  that  the  pdf's  for 
X  and  A  given  by  Equations  (1.4)  and  (1.7)  respectively 
are  indeed  supported  by  computer  simulation.  Of  the 
24  examples  shown,  only  two  fail  the  Chi-square  test. 

All  the  pdf's  in  Figures  13-36  are  the  limiting  densities 
as  M  -*■«.  The  results  of  the  simulation  for  V  *  2  and  3 
paths  (including  further  analysis  of  this  case)  are 
given  in  Appendix  C.  Upon  comparison  with  the  results  of 
Appendix  C,  it  is  quite  acceptable  to  assume  that  for 
N  _>  4  paths  phase  random  multipath  propagation  is 

obtained.  It  is  true  that  x  for  W  =  4  paths  did  fail 
the  Chi-square  test  but  this  is  only  one  out  of  eight 
variables.  It  should  be  noted  that  for  .V  =  2  and  3  paths 
the  limiting  densities  for  seven  out  of  eight  of  the 
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variables  (the  multipath  phase  <?  is  always  uniform) 
failed  the  Chi-square  test  (see  Reference  [27]),  thus  the 
justification  for  choosing  SI  _>  4  as  the  required  number  of 
paths  is  evident.  The  failure  of  $  for  the  5  path  case 
indicates  perhaps  a  lack  of  randomness  or  a  correlation 
for  that  particular  cycle  length  in  the  random  number 
generator.  Though  the  results  on  the  average  get  better 
as  N  becomes  large,  as  indicated  in  Table  V,  (except 
curiously  enough  for  A) ,  after  12  paths  the  results 
are  approximately  constant  indicating  that  the 
correlation  in  each  of  the  random  number  generators  is 
now  the  limiting  factor.  This  assertion  follows  from 
consideration  of  M  *  30  [27]  in  which  by  the  wav  x2  for 
A  is  4.49,  and  x£s  >05  =  26.30.  In  addition  to 
containing  the  other  simulations  not  shown  in  the  figures 
nor  mentioned  above  Reference  [27]  also  includes  the 
more  detailed  printout  mentioned  earlier.  It  is 
interesting  to  note  that  even  when  phase  random  conditions 
are  met  for  the  single  source  the  pdf  for  p  always  appears 
to  be  among  the  worst  performers  while  p  is  always  among 
the  best.  This  result  was  first  noted  by  Hamblen  [9] 
when  he  investigated  ocean  acoustic  data.  Because  this 
result  is  repeated  in  the  simulation,  it  appears  to  be 
a  function  of  the  statistical  nature  of  the  variables,  as 
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well  as  the  non-stationarity  of  when  analyzing  real 

data. 


2.3  Simulation  Results  for  Multiple  Sources 

Figures  37-42  are  the  results  of  a  simulation  run 
with  20  equal  intensity,  equal  v2  receptions,  5  paths, 
and  I  =  600  samples.  The  theoretical  densities  are  given 
by  Equations  (1.19),  (1.67),  (1.79),  (1.80),  (1.83),  and 

(1.84)  respectively.  The  results  of  the  Chi-square  test 
are  tabulated  in  Table  VI.  Also  shown  are  the  means 
and  standard  deviations  as  predicted  by  the  theory  and 
observed  in  the  simulation.  With  one  exception,  the 
computer  simulation  agrees  well  with  the  theoretical 
predictions.  One  of  the  important  results  of  the 
analysis  of  the  equal  intensity,  equal  v2  case 
(Section  1.2. 2.1)  is  that  P^(J)  is  independent  of  the 
number  of  sources.  Figures  24,  43,  44,  and  40  are 
examples  of  P^(f)  for  1/  3,  4,  and  20  equal  intensity 
(a2  =  2.5),  equal  v2  (v  =  .007  rad/sec)  sources. 

Figures  43  and  44  passed  the  Chi-square  test,  and  the 
histograms  clearly  support  this  prediction  as  well. 

As  a  final  example  of  multiple  source  simulations, 

I  treated  the  N  =  15  group  problem  of  Reference  [4].  As 
was  noted  in  Section  1.2* 1.4,  and  illustrated  in  Table  I, 
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TABLE  'VI 


Results  of  the  chi-square  test  for  RANDPHASE 
simulation  of  20  equal  intensity  (o?  =  2.5), 
equal  v2  sources  with  M  =  5  for  each.  The  input 
v  was  .00700.  Actual  v  was  .00699  rad/sec. 

The  units  are  the  same  as  in  the  figures. 


Fig. ' 

Variable 

x2 

X2:  .05 

Mean 

Std. 

Dev. 

Theory 

Data 

Theory 

Data 

37 

X 

9.26 

26.30 

100.00 

99.36 

22.36 

20.81 

33 

• 

X 

17.01 

23.63 

0.00 

-0.01 

0.22 

0.20 

39 

0 

11.08 

26.30 

9.94 

9.94 

1.11 

1.  40 

40 

• 

p 

19.62 

13.31 

0.00 

0.00 

0.01 

0.01 

41 

A 

15.26 

25.00 

19.89 

19.90 

0.98 

0.93 

42 

• 

A  • 

13.53 

26.30 

0.00 

0.00 

0.01 

0.01 

Computer  generated  histogram  of  p  for  L  s  20  sources 
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Fig.  43  Computer  generated  histogram  of  p  for  L  =3  sources 
and  equ.  (1.60) 


Computer  generated  histogram  of  p  for  L  s  A  sources 


the  first  eight  groups  are  the  most  energetic.  I 
applied  BURMRAN  therefore,  with  123  sources,  the 
appropriate  r  array,  5  paths,  I  »  600  samples,  and  an 
input  v  of  .007  rad/sec.  The  results  are  shown  in 
Figures  45,  46,  47,  and  47a  for  x»  X/  A,  and  $ 
respectively.  In  Table  VII  I  have  listed  the  Chi-square 
results  and  other  pertinent  data.  The  theoretical 
densities  are  given  by  the  Edgeworth's  series  approximations 
Equations  (1.45),  (1.108),  and  (1.53)  for  x,  and  A 
respectively.  For  5  I  have  plotted  the  Gaussian  using 
the  variance  obtained  from  the  data  (the  curve  was  not  fit) . 
The  excellent  performance  confirms  the  assertion  made 
in  Section  1.2  that  even  for  the  most  complicated 
realizations  of  Case  (c)  the  pdf  for'  p  will  be  Gaussian 
in  which  the  variance  is  in  fact  *  E[p2].  The  results 
for  x  and  A  are  not  very  good  (once  again  the  rate 
variables  out-perform  the  amplitude  variables) .  The  error 
appears  to  be  in  underpredicting  the  mean  as  indicated 
in  Table  vil.  However,  these  results  should  be  inter¬ 
preted  in  light  of  the  very  large  number  of  samples 
required  for  9r  and  §n  (369,000  for  each)  certainly 
introduced  some  correlation  error  due  to  the  repeat 
characteristics  of  the  random  number  generator  (see 
Table  IV) .  Both  0  and  9  barely  passed  the  Chi-square 
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TABLE  VII 


Results  of  the  chi-square  test  for  3URMRAN 
simulation  of  123  different  source  problem. 
The  units  are  the  same  as  in  the  figures. 


Fig. 

Variable 

x2 

Xn:  *05 

Mean 

Theory  Data 

Std. 

Theory 

Dev. 

Data 

45 

X 

36.46 

23.68 

5.72 

5.80 

1.65 

1.43 

46 

• 

X 

18.72 

23.68 

0.00 

0.00 

0.02 

0.01 

47 

A 

36.34 

27.59 

7.36* 

7.51 

1.46* 

1.07 

47a 

• 

P 

13.05 

30.14 

0.00 

0.00 

— 

o.oo: 

♦From  equ.  (1.60a)  which  is  the  Gaussian  assumption.  Note 
that  these  numbers  should  not  be  compared  with  those  given 
in  Section  1.2.1  because  this  simulation  included  only  the 
first  8  groups  so  there  will  be  a  slight  difference  in  the 
mean  values  (see  Table  I) . 
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Computer  generated  histogram  of  X  for  the  123  different 
source  problem  modelling  distant  shipping  noise  off  Bermuda, 


Fig.  46  Computer  generated  histogram  of  X  for  the  L  =  123  different 
source  problem  modelling  distant  shipping  noise  off  Bermuda, 


Fig.  47  Computer  generated  histogram  of  the  level  for  the  LSI23 

different  source  problem  modelling  distant  shipping  noise  off 
Bermuda,  and  equ.  (1.53). 


i&r 
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p  volts  -  rad /sec 
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Fig.  47a  Computer  generated  histogram  of  p  for  the 
LSI23  different  source  problem  modelling 
distant  shipping  noise  off  Bermuda,  and  the 
Gaussian  using  a  variance  given  by  E[/52| 
obtained  from  the  computer  generated 
samples. 
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test  with  results  of  x2  =  994.16  and  981.53  respectively, 
and  .03  *  1019.94.  In  any  case,  the  non-Gaussian 
skew  is  certainly  evident  in  the  histogram  in  Figure  45, 
once  again  dissuading  the  Gaussian  assumption  in  favor 
of  the  Edgeworth's  series. 

2 . 4  Simulation  Results  for  a  Finite  Bandwidth 

and/or  Modulated  Source 

Figures  43-51  are  examples  of  the  RANDPHASE 
simulation  of  when  the  source  is  undergoing  varying 
degrees  of  sinusoidal  phase  modulation.  The  histograms 
are  plotted  with  Equation  (1.169)  (Note  that  v  and  3a 
are  both  given  in  rad/sec) .  Figures  52-55  are  examples  of 
varying  degrees  of  uniform  frequency  modulation,  or 
alternatively  varying  amounts  of  bandwidth  on  ®._,.  The 
histograms  are  plotted  with  Equation  (1.176)  (Note  that 
v  and  A  are  both  given  in  rad/sec) .  In  Table  VIII  I 
have  compiled  the  results  of  the  Chi-square  test,  and  the 
values  of  8a,  A,  or  B,  as  appropriate.  All  these 
simulations  were  run  with  an  input  value  of  v  =  .007  rad/sec, 
5  paths,  and  I  *  700  samples.  With  the  use  of  Table  VIII 
compare  Figure  5  2  (for  B  =  v)  and  Figure  53  (for  B=2v), 
with  Figure  28,  the  non-mcdulated  pdf  for  ? .  This  shows 
quite  clearly  again  that  only  for  B  <<  2v  can  bandwidth 
effects  be  neglected.  As  with  the  other  simulation 
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TABLE  VIII 


Results  of  chi-square  test  and  other  data  for  the 
RANDPHASE  simulation  of  for  varying  degrees  of 
modulation'  and  bandwidth  with  W  »  5  paths. 


Fig. 

x2 

X2:  .05 

V 

rad/sec 

A 

rad/sec 

B 

rad/sec 

8  a 

rad/sec 

48 

19.30 

23.68 

.0070 

- 

- 

.0011 

49 

31.55 

23.68 

.00  70 

- 

- 

.0070 

50 

12.13 

23.68 

.0070 

- 

- 

.0210 

51 

22.39 

28.87 

.0070 

- 

- 

.0840 

52 

25.32 

26.30 

.0070 

.0035 

.0070 

- 

53 

12.72 

23.68 

.0070 

.0070 

.0140 

- 

54 

17.62 

23.68 

.0070 

.0210 

.0420 

- 

55 

22.76 

28.  87 

.0070 

.0700 

.1400 

— 

Histogram  of  computer  generated  samples  of  <f>  with  sinusoidal  phase 
modulation  of  the  source,  v*.007  ,  /3a=  .0011  and  equ.  (1.169) 


MULTIPATH  PHASE  RATE  RAD/SEC 


MULTIPATH  PHASE  RATE 


Histogram  of  computer  generated  samples  of 
sinusoidal  phase  modulation  of  the  source 


MULTIPATH  PHASE  RATE  RAD/SEC 


Histogram  of  computer  generated  samples  of  with  uniform  frequency 
modulation  of  the  source  ,  v  -  .007  .  A=  .021  and  equ.  ( 1.176) 


Histogram  of  computer  generated  samples  of  with 
uniform  frequency  modulation  of  the  source  , 
i/  =  .007,  A=.07and  equation  (1.176). 


results,  the  fit  of  the  theory  to  the  computer  data  is 
quite  good  with  only  one  of  the  eight  runs  failing  the 
Chi-square  test. 


CHAPTER  3 


DATA  ANALYSIS 

I  have  analyzed  two  sets  of  data  from  acoustic 
experiments  in  the  ocean.  These  data  are  compared  with 
the  theory  presented  in  Chapter  1.  In  Section  3.1,  data 
from  an  experiment  performed  by  R.  Porter  and  R.  Spindel, 
near  Eleuthera  [12]  are  shown  among  other  successful 
comparisons  to  support  the  theoretical  pdf's  for 
amplitude  rate. and  level  rate  derived  for  the  first  time 
in  this  thesis.  Also,  the  analytical  results  of  the 
modulation  theory  derived  in  Section  1.3  explain  heretofore 
unobserved  phenomena  of  the  Eleuthera  data  accurately. 

In  Section  3.2,  data  from  the  CASE  experiment  [13] 
conducted  in  the  Pacific  are  investigated  in  general, 
with  special  emphasis  in  light  of  the  modulation  theory 
of  Chapter  1.  Furthermore,  these  data  provide  more 
insight  into  the  parameter  v2  and  reveal,  as  well, 
shortcomings  of  our  current  understanding  of  this  vital 
ocean  acoustic  parameter  and  its  driving  mechanisms. 

3 . 1  The  Eleuthera  Experiment 

Data  made  available  by  the  Woods  Hole  Oceanographic 
Institution  were  acquired  during  a  long-range  acoustic 
propagation  experiment  conducted  near  Eleuthera  [12] . 
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The  data  analyzed  consisted  of  four  records  (see  Table  IX) 
in  which  two  cw  signals,  one  at  220  Hz,  and  one  at  406  Hz, 
for  each  record  were  transmitted  from  Eleuthera  to 
drifting  sonobuoys  approximately  300  km  northeast  towards 
Bermuda.  A  Doppler  position-tracking  system  [30]  was 
used  to  remove  mean  multipath  phase-rates  due  to 
sonobuoy  motion.  These  data  were  also  analyzed  by  Dyer 
and  Shepard  [15]  and  Hamblen  [9] . 


TABLE  IX 

Log  of  the  W.H.O.I.  Fluctuation  Data 


Record 

Date 

Time  (GMT) 

Record  Length  (h) 

447 

13 

Sept.  74 

1900  -  0310 

8.16 

448 

12 

Sept.  74 

0402  -  1204 

8.05 

449 

11 

Sept.  74 

1800  -  0145 

7.75 

424 

11 

Sept.  74 

0644  -  1527 

•  8.72 
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3.1.1  Analysis  of  the  Single  Source 
From  the  original  digitized  time  series  of  pcos<J>  and 
psin<J>,  time  series  of  A,  x>  and  <p  for  the  two  frequencies 
are  produced  for  each  record.  The  phase  data  were 
analyzed  previously  by  Dyer  and  Shepard  [15]  and 
Hamblen  [9] .  Details  of  the  phase  unwrapping  routine 
can  be  found  in  Reference  [15] .  Histograms  are  then 

generated  from  these  time  series.  For  the  histograms  of 

•  • 

A  and  x»  Equations  (1.7)  and  (1.4),  respectively,  are  fit 

such  that  the  value  of  v2  and  c£v,  respectively,  minimize 

the  Chi-square  statistic.  The  value  of  v2  for  each  run 

is  also  obtained  by  fitting  Equation  (A8)  to  the 

histograms  of  <p.  The  values  of  v2  obtained  from  <+> 

agree  closely,  as  expected,  with  those  obtained  previously 

by  Hamblen  [9] .  However,  the  value  of  v2  obtained  from 
• 

A  by  fitting  Equation  (1.7)  is  smaller  by  a  factor  of 

almost  2  in  every  case.  This  can  be  understood  as 

follows.  As  the  analysis  of  Section  1.3  shows,  any 

finite  bandwidth  effects  or  carrier  instabilities,  or  in 

this  particular  case,  any  errors  in  the  Doppler 

position-tracking  system,  would  manifest  themselves  in 

"spreading”  or  increasing  the  variance  of  p  while  leaving 
• 

A  unmolested. 


To  test  this  explanation  further,  v2  is  obtained  for 
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each  record  by  taking  the  ratio  of  c£v2  to  a|,  parameters 
which  were  found  previously  [9]  by  fitting  Equation  (Al) 
and  (A5)  to  histograms  of  p  and  p,  respectively.  Also, 
v2  is  obtained  by  taking  the  value  of  c£v  from  the  fit 
of  Equation  (1.4)  to  the  histograms  of  x  and  dividing 
by  the  previously  found  values  of  c2  [9] .  If  the 
reasoning  is  correct,  then  these  values  of  v2  should  agree 
with  the  values  of  v2  obtained  from  A  for  each  record  as 
they  are  all  derived  from  densities  which  are  insensitive 
to  angle  modulation  or  finite  bandwidth  effects.  The 
results  are  tabulated  in  Table  X,  which  shows  that 
without  exception  these  values  of  v2  agree  to  within  33% 
or  less  and,  furthermore,  they  are  all  as  hypothesized 
less  than  the  values  of  v2  obtained  from 

Except  for  Record  424,  the  values  for  v  obtained 
from  A,  noting  the  square  root  range  dependence  predicted 
by  Dyson,  Munk,  and  Zetler  [11] ,  are  reasonably  consistent 
with  the  values  for  v  obtained  previously  by  them  for 
the  Eleuthera  to  Mid-station  and  the  Eleuthera  to  Bermuda 
transmissions  for  the  406  Hz  source,  while  the  values 
of  v  obtained  from  <p  are  not  as  consistent  (see  Table  XI)  . 
Why  run  424  exhibited  the  larger  values  for  v  in  all 
cases  (see  Table  X  again)  is  unknown. 

In  order  to  obtain  a  measure  of  the  bandwidth  of  the 


*  W?  ~*5*V»**£**j*K^  * 


-217- 


TABLE  X 


The  values  of  v2  obtained  from  the  time  series  and 
histograms  of  the  amplitude  densities  are  compared 
with  the  values  of  v2  obtained  in  fitting 
Equation  (A8)  to  the  histograms  of  $.  The  values 
of  A,  the  half-bandwidth  measure  of  the 
"modulation”  is  tabulated  for  each  run  in  the  last 
column. 


Freq. 

(Hz) 

V2 

rad2 /sec2  (x 

IQ-4) 

Record 

a2v2/c2 

(a2\>/<J2) 2 

• 

A 

0 

A,  mHz 

a_,  mHz 

G 

220 

1.8 

1.5 

2.0 

4.7 

3.3 

2.0 

447 

406 

3.5 

3.6 

4.1 

7.9 

4.3 

2.6 

220 

6.7 

8.3 

6.1 

8.7 

2.2 

1.3 

448 

406 

2.9 

2.4 

3.0 

7.1 

4.4 

2.7 

220 

2.4 

2.4 

2.7 

4.4 

2.6 

1.6 

449 

406 

3.1 

2.4 

3.6 

7.6 

4.4 

2.7 

220 

10 

12 

11 

15 

2.8 

1.7 

424 

406 

11 

11 

15 

24 

5.8 

3.6 
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TABLE  XI 


The  average  values  of  v  in  mHz  for  the  406Hz 
source  for  Records  447,  448,  and  449 
obtained  from  4)  and  A  are  compared  to  the 
values  of  v  obtained  by  Dyson,  Munk,  and 
Zetler  (DMZ) . 


i 

9 

A 

DMZ 

Mid-Station 

DMZ 

Bermuda 

Range 

(km) 

300 

300 

550 

1250 

Measured  Equ.  Ray  Mix 

Measured 

v  in 

mHz 

4.4 

3 .0 

2.8  ■  (18)  2.6 

4.0 

2.2 

modulation  in  these  data,  I  fit  Equation  (1.176)  to  the 
histograms  of  4>  using  the  values  of  v2  obtained  from  the 
histograms  of  A.  I  determine  the  value  of  A,  the  half¬ 
bandwidth  of  the  modulation,  that  minimizes  the  Chi-square 
statistic.  I  likewise  fit  Equation  (1.178).  which  assumes 
a  Gaussian  rather  than  a  -uniform  modulation  function 
and  find  a_,  the  standard  deviation  of  the  carrier 

V3 

fluctuations.  Note  that  the  half-bandwidth  of  the 
modulation  as  given  by  1/e  of  the  best  fit  of 
Equation  (1.178)  is  approximately  equal  to  the  value  of 

A  obtained  from  fitting  Equation  (1.176) .  The  values  of 

*  •  • 
vz  from  the  A  histograms  are  used  because  A,  as  revealed 

by  Equation  (1.7),  is  independent  of  cr£  and  is  less 

sensitive  to  its  non- stationary  behavior.  In  Table  X, 

the  least  Chi-square  value  of  A  and  <?G  are  given  for  each 

run.  Because  the  frequency  stability  of  both  sources 

is  many  orders  of  magnitude  less  than  v,  the  modulation 

as  measured  by  A  or  can  be  attributed  to  the  error  in 

the  Doppler  position-tracking  system.  An  error  of  ±.017 

m/sec  in  measuring  the  velocity  of  the  receiving 

hydrophone  itself  results  in  a  half-bandwidth  value  for 

the  modulation  of  2.5  mHz  at  220  Hz  and  4.5  mHz  at  406  Hz, 

consistent  with  the  values  tabulated  in  Table  X.  The 

velocity  error  is  consistent  with  the  experimental  setup 
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and  the  performance  characteristics  of  the  Doppler 
position-tracking  system  [31] . 

•  • 

In  Figures  56-59,  the  histograms  of  x,  A,  and  $  and 
the  least  squares  fit  of  Equations  (1.4),  (1.7),  and 
(1.176),  respectively,  are  shown  for  all  the  records. 

Table  XII  lists  the  Chi-square  statistics  for  A  and  x 
and  Table  XIII  lists  the  Chi-square  statistics  for  $  when 
fit  by  Equation  (A8)  (which  assumes  no  modulation) , 
Equation  (1.176),  and  Equation  (1.178).  With  a  few 
exceptions,  the  fit  of  the  theory  to  the  data  is 
excellent.  Of  the  three  runs  that  failed  the  Chi-square 
test,  two  were  the  result  of  fitting  the  no-modulation 
density.  Equation  (A8) ,  to  $  which,  in  light  of  the 
foregoing  analysis,  is  suspect  from  the  start.  As 
Table  XIII  shows,  while  Equation  (1.176)  performed  better 
than  Equation  (A8)  in  five  out  of  the  eight  runs,. 

Equation  (1.178)  performed  better  than  both  Equations  (A8) 
and  (1.176)  for  all  cases  except  Run  424  where  performance 
was  better  than  Equation  (1.176)  but  slightly  worse  than 
Equation  (A8) .  In  fact,  a  Gaussian-like  error  in 
extracting  the  mean  phase  rate  is  characteristic  of  the 
Doppler  position-tracking  system  [31] . 
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fig.  38  Histograms  from  fttcord  449  for  x,  A,  and  ond  ins  isost  chi-souars  fit  of  sous  (1.4).  (1.7). and 
U.I76 1,  rtsptctivsly.  Rtsuits  for  220  Hz  art  on  ths  lift,  ond  406  Hz  on  tnt  right. 
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TABLE  XII 

Results  of  the  Chi-square  goodness-of-f it  test 
of  x  and  A  to  Equations  (1.4)  and  (1.7), 
respectively. 


Record 

Freq. 

Hz 

•  2 
X*X 

X^:.0S 

A:  X2 

X“:.05 

220 

1.81 

16.92 

3.45 

21.03 

447 

406 

3.52 

23.68 

1.86 

23.68 

220 

1.45 

31.41 

11.61 

28.87 

448 

406 

2.61 

19.68 

4.51 

23.68 

220 

1.65 

21.03 

1.20 

23.68 

449 

406 

2.96 

21.03 

6.44 

23.68 

424 

220 

2.30 

31.41 

32.31 

28.87 

406 

6.04 

31.41 

16.38 

28.87 
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TABLE  XIII 

Results  of  the  Chi-square  goodness-of-f it  test 
of  |  to  Equations  (AS ) ,  (1.176),  and  (1.178). 


Record 

Freq. 

Hz 

Eq.  (A8) : 

x2 

Eq.  (1.176): 
X2 

Eq.  (1.178): 

x2 

X2:.°5 

220 

16.24 

3.19 

2.78 

15.51 

447 

406 

16.54 

4.13 

1.19 

21.03 

448 

220 

406 

10.57 

22.95 

21.10 

6.34 

7.06 

3.68 

23.68 

21.03 

449 

220 

5.43 

3.26 

2.67 

21.03 

406 

14.87 

7.02 

3.67 

21.03 

220 

10.45 

23.84 

12.07 

26. 

.30 

424 

406 

9.78 

14.61 

12.71 

26. 

.30 
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3.1.2  Multiple  Sources 

Because  the  220  Hz  and  406  Hz  sources  were  transmitted 
and  recorded  simultaneously  for  each  record  and  were  each 
quadrature  demodulated  about  their  respective  center 
frequencies  and  summed  incoherently,  this  affords  an 
opportunity  to  check  the  equations  derived  in  Section  1.2. 2. 2 
for  two  different  sources.  For  each  record,  the  220  Hz 
source  was  designated  Source  one,  and  the  406  Hz, 

Source  two.  The  values  for  and  were  obtained 
from  Hamblen's  [91  previous  analysis.  The  values  for 
and  rad  /sec  are  obtained  from  Table  X.  In 

Figures  60-63,  I  have  plotted  representative  results. 

In  Figure  60  the  histograms  for  x  ana  p  are  plotted  with 
the  following  equations: 


P  (X) 
X2D 


2aii  -  2a‘2 


exp(-  j£r~)  '  ex?(~ 


11 


2j12 


(3.1) 


P. 


'2D 


(P) 


r*  —  nr* 
11  °12 


exp  (- 


)  -  exp(- 


p2 

5®n 


(3.2) 


Equation 

Equation 

Equation 


(3.1)  was  first  derived  by  Dyer  [4],  and 

(3.2)  is  simply  the  transformation  p*x  of 
(3.1).  In  Table  XIV,  I  have  listed  the  results  of 


r 


PX(X)  Pplp) 


im  of  X  for  two  different  sources  from  Record  449,  plotted 
95)  with  L  =  2. 
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Fig.  63  The  histogram  of  A  for  two  different  sources 
from  Record  447  plotted  with  equ.  (1.123). 


TAUIJ3  XIV 
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the  Chi-square  test  and  listed  the  values  of  and  v: 
used.  Note  that  the  theoretical  pdf's  for  these  plots 
were  not  fit  to  the  histograms.  Also  in  determining  the 
Chi-square  statistic,  I  have  approximated  the  area  under 
the  theoretical  pdf  for  each  class  interval  by  taking  the 
value  of  the  pdf  at  the  center  of  the  interval  times  the 
width  of  the  interval.  Figures  61-63  are  the  histograms 
for  x,  p,  and  A,  plotted  with  Equations  (1.95),  (1.114), 
and  (1.123),  respectively. 

In  general,  the  amplitude  densities  perform  rather 
poorly  as  expected,  while  the  rate  densities  perform 
exceptionally  well  supporting  the  theory  in 
Section  1.2. 2. 2.  As  noted  earlier,  the  non-stationarity 
in  a*  discovered  by  Hamblen  [9]  in  these  data,  as  well  as 
the  apparent  statistical  stability  of  the  rate  variables 
accounts  for  the  difference  in  performance  of  the 
amplitude  and  amplitude  rate  variables. 

3.1.2  Crossing  Rate  Statistics 

Following  Dyer  and  Shepard  [15] ,  I  have  analyzed 
the  crossing  rate  statistics  for  each  of  the  frequencies 
of  the  four  records  of  the  Eleuthera  experiment.  Even 
though  the  "modulation"  in  these  data  are  very  small 
(see  Table  X) ,  the  effect  on  the  crossing  rate  statistics 


1 


is  dramatic.  From  Equation  (1.18),  I  have  [15], 


G(P0) 


c^/27 


exp  [ 


(3.3) 


and  [15], 


G<V 


v 

2r 


(3.4) 


Equation  (3.4)  assumes  no  modulation, 
of  Equation  (3.3)  to  Equation  (3.4)  is 
i.e.  [15], 


G^) 

G<V 


(8rr) 


P  (pj 
,0  o 


When  the  ratio 
taken,  v  cancels, 


(3.5) 


However,  if  the  source  is  modulated  or  has  a  finite 
bandwidth  (i.e.,  B  *  2v) ,  then  in  fact  v  does  not  cancel 
when  the  ratio  of  the  two  is  taken: 


G(Pp> 

G(*0) 


(87r)*  71  Vpo)al 

9 


(3.6) 


where  v A  signifies  the  value  of  v  obtained  from  the 
amplitude  densities,  and  v*  the  value  of  v  obtained  from 
fitting  Equation  (A3)  (which  assumes  no  modulation)  to  the 
histogram  of  $  of  a  modulated  source.  This  is  tantamount 
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to  assuming  that  the  pdf  for  M  is  such  that  Equation  (1.165) 

yields  the  Longuet-Higgins  density  with  a  larger 

v  =  v*  >  v_.  If  there  is  no  modulation  of  the  source 
0  A 

and  B  <<  2v,  then  v.  * v*  (as  will  be  shown  in  Section  3.2). 

A  $ 

For  the  data  of  the  Eleuthera  experiment,  however, 

vat<v*  because  of  the  error  in  the  Doppler  position/tracking 

system.  The  ratio  vA/v*  can  be  considered  a  correction 

factor  to  account  for  modulation  effects  in  crossing  rate 

statistics.  However,  because  the  error  in  extracting 

the  mean  phase  rate  is  more  Gaussian  than  Longuet-Higgins, 

as  demonstrated  quite  convincingly  by  the  performance  of 

Equation  (1.178)  on  the  histograms  of  4>,  I  have  also 

computed  the  correction  factor,  using  the  analysis  of 

Section  1.3.3,  given  by  Equation  (1.209).  In  Table  XV, 

I  have  compiled  the  correction  factors  given  by  \>A/v*, 

and  Equation  (1.209).  I  obtain  vA/v*  by  taking  the  ratio 
•  • 

of  v  from  A  to  v  from  $  given  in  Table  X.  For 
Equation  (1.209)  cG,  and  v  (from  A,  converted  to  mHz)  are 
likewise  obtained  from  Table  X.  In  Figures  64-67 
reproduced  in  part  from  Reference  [15]  I  have  plotted 
the  new  curves  applying  the  vA/v^,  and  Equation  (1.209) 
correction  factors.  Note  that  because  the  correction 
factors  are  a  function  of  v  and  parameters  of  the 
modulation  there  is  one  curve  for  each  frequency.  The 


i 


236- 


TABLE  XV 


The  correction  factors  to  the  unmodulated 
value  of  G(p  )/G(4>  )  for  the  Eleuthera 
experiment.  Based  on  the  ratio  of  v. 

(from  A,  Table  X),  to  v*  (from 
Table  X) ,  and  a  Gaussian  frequency 
modulation  assumption. 


Record 

Frequency 

V»i 

Eq.  (1.209) 

220 

.65 

.41 

447 

406 

.72 

.44 

220 

.84 

.56 

448 

406 

.65 

.40 

220 

.78 

.47 

449 

406 

.69 

.42 

424 


220 


86 


58 


Equ.  (3.S) 


RECORD  447 

o  220  Hz 
x  406  Hz 


406  Hz  (  i^/fy) 
220  Hz 


406  Hz  Gaussian 
~  220  Hz  Gaussian 


Fig. 64  Ratio  of  phase  period  to  amplitude  period  ,  versus  amplitude , 
for  record  447.  Equ.  (3.5)  is  the  theoretical  ratio  fora 
non- modulated  narrowband  source.  The  correction  factors 
(Table 3Y),  and  Equ.  (1.209)  (Goussian  modulation) have 
been  applied  to  obtain  the  corrected  curves  for  220  Hz  and 
406Hz  accounting  for  the  modulation. 
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Fig.  65  Ratio  of  phase  period  to  amplitude  period  ,  versus  amplitude, 
for  record  448.  Equ.  (3.5)  is  the  theoretical  ratio  fora 
non-modulated  narrowband  source.  The  correction  factors 
(Table  and  Eau*  (1209)  (Gaussian  modulation)  have 
been  applied  to  obtain  the  corrected  curves  for  220Hz  and 
406Hz  accounting  for  the  modulation. 
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Fig.  66  Ratio  of  phase  period  to  amplitude  period,  versus  amplitude, 
for  record  449.  Equ.  (3.5)  is  the  theoretical  ratio  fora 
nan-modulated  narrowband  source.  The  correction  factors 
(Table  22),  and  Equ.  (1.209)  (Gaussian  modulation)  have 
been  applied  to  obtain  the  corrected  curves  for  220 Hz  and 
406  Hz  accounting  for  the  modulation. 
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Ratio  of  phase  period  to  amplitude  period  ,  versus  amplitude  , 
for  record  424.  Equ.  (3.5)  is  the  theoretical  ratio  for  a 
non-modulated  narrowband  source.  The  correction  factors 
(Table  EC),  and  Equ.  (l.209)(Gaussian  modulation)  have 
been  applied  to  obtain  the  corrected  curves  for  220 Hz  and406Hz 
accounting  for  the  modulation. 
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bracketed  < G  ( i>0 )  >  indicates  that  an  average  value  of 
G(4>0)  has  been  used.  The  range  of  values  results  from 
restrictions  placed  on  the  maximum  value  of  5  as  a 
consequence  of  phase  unwrapping  routines.  As  indicated 
in  the  figures,  pQ  has  been  normalized  such  that  the 
E[Pq]  ■  2.  Details  of  the  above  and  other  aspects  of  the 
data  reduction  can  be  found  in  Reference  [15] . 

The  corrected  curves  are  a  dramatic  improvement 
over  the  Dyer,  Shepard  theory.  Except  for  run  424  where 
the  Gaussian  performs  much  better,  the  Longuet-Higgir.s 
and  Gaussian  curves  perform  about  the  same  with  the 
former  tending  to  be  slightly  high  and  the  latter 
slightly  low.  The  variance  in  these  results  only 
indicates  our  uncertainty  in  the  exact  nature  of  the 
modulation  in  these  data.  The  important  results  are  as 
follows:  (1)  The  modulation  theory  discovered  the 

heretofore  unnoticed  modulation  in  these  data, 

(2)  accounted  for  its  effects  on  the  histograms  of  5, 
and  (3)  correctly  predicted  the  crossing  rate  statistics 
completely  explaining  the  consistent  over  prediction  of 
the  Dyer,  Shepard  theory. 


3 . 2  The  CASE  Experiment 

Data  from  the  CASE  experiment  [13]  were  analyzed  in 
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order  to  compare  the  theory  to  data  in  which  the 
modulation  is  known  and  controlled.  Three  configurations 
were  used:  (1)  fixed  source  on  a  seamount,  (2)  source 
towed  by  a  surface  ship,  and  (3)  source  mounted  on  a 
submersible.  The  signals  were  monitored  at  four  widely 
separated  fixed  deep  water  receivers  at  ranges  varying 
from  200  to  400  kn.  Runs  from  one  of  the  receivers  had  to 
be  subsequently  rejected  due  to  an  extremely  low  signal  to 
noise  ratio.  Two  carrier  frequencies  were  employed,  one 
at  15  and  the  other  at  33  Hz.  The  signals  were  frequency 
modulated  by  a  pseudo-random  function  generator  with  an 
average  period  of  107  sec  and  an  overall  pattern  that 
repeats  itself  every  640  sec  (see  Figure  68) .  The 
bandwidth  of  the  modulation  is  pre-selected,  and  for  the 
runs  analyzed  is  either  0  (no  modulation),  1/8,  1/4, 
or  1/2  Hz. 

A  total  of  88  runs  of  15  min  duration  each  was 

analyzed.  From  the  filtered  digitized  complex  data,  time 

•  • 

series  and  histograms  were  generated  for  A  and  0 .  As 
before.  Equations  (1.7)  and  (AS)  are  fit,  respectively, 
such  that  the  value  of  v2  obtained  minimizes  the  Chi- 
square  statistic. 

Of  the  88  runs,  12  were  not  modulated  and  the  values 
•  • 

of  v  obtained  from  A  and  }  for  those  runs  agree  (run  by  run) 


1 


to  within  approximately  10%.  In  Table  XVI,  I  have  listed 

these  values  as  well  as  the  range  frequency  and  receiver 

number  for  each  run.  For  all  of  these  runs,  the  source 

was  fixed  on  the  seamount.  For  9  of  these  12  runs,  the 

quadrature  components  exhibit  significant  nonzero  means. 

During  the  data  reduction  these  means  for  all  the  runs 

analyzed  regrettably  were  removed.  However,  runs  76-77 

had  no  significant  means  and  the  values  of  v  are 

consistent  with  those  from  runs  66-68  which  are  at  the 

same  frequency.  It  appears  that  extracting  the  means 

even  when  significant  (and  this  was  true  only  for  runs 

66-68,  71-73,  and  113-115)  is  a  negligible  factor. 

However,  runs  113-115  do  in  fact  contain  the  largest 

three  values  of  v  obtained  for  all  the  15  Hz  runs  and  are 

perhaps  suspect,  though  other  runs  with  insignificant 

means  also  exhibit  values  of  v  approaching  and  exceeding  1 

For  the  remaining  76  runs,  24  were  modulated  with  a 

1/8  Hz  bandwidth,  35  with  a  1/4  Hz  bandwidth,  and  17  with 

a  1/2  Hz  bandwidth.  As  expected  with  modulation,  the 

value  of  v2  obtained  from  fitting  Equation  (A8)  to  the 

histograms  of  <p  were  larger  than  those  obtained  from 
• 

A  by  large  factors  depending  upon  the  bandwidth  of  the 
modulation.  Because  the  modulation  pattern  (Figure  68) 
is  approximately  a  saw-tooth  and  the  spectrum  of  the 
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TABLE  XVI 


The  values  of  v  obtained  from  the  12  non- 
modulated  runs  analysed  in  the  CASE  experiment. 
In  addition  to  being  non-nodulated  the  source 
was  fixed  on  a  seamount. 


Case 

Record 

* 

Frea. 

(Hz)' 

V 

(rad/sec) 

Range  (km) 

RCVR 

• 

* 

V 

A 

66* 

33 

.094 

.  104 

250 

1 

67* 

33 

.353 

.354 

320 

2 

63* 

33 

.191 

.194 

450 

3 

71* 

15 

.407 

.421 

250 

1 

72* 

15 

1.136 

.  982 

320 

2 

73* 

15 

.741 

-727 

450 

3 

76 

33 

.042 

.044 

250 

1 

77 

33 

.  154 

•  147 

320 

2 

73 

33 

.106 

.096 

450 

3 

113* 

15 

1.349 

1.20.8 

250 

1 

114* 

15 

1-495 

1.272 

320 

2 

115* 

15 

1.279 

1.084 

450 

3 

*These  runs  showed  significant  means  in  the  quadrature 
components . 
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modulation  (Figure  69)  is  approximately  uniform,  the 
analysis  of  Section  I.B  is  applied. 

Equation  (1.176)  is  fit  to  the  histograms  of  $ 
using  the  value  of  v2  obtained  from  fitting  A  and  varying  A 
such  that  the  Chi-square  statistic  is  minimized.  The  value 
of  A  is  then  compared  to  the  actual  signal  bandwidth  set 
by  the  experimenter.  In  Figure  70  the  error  in  percent 
is  plotted  against  the  number  of  runs  that  exhibit  a  given 
error.  For  the  76  runs  with  modulation,  the  average 
error  in  determining  the  bandwidth  of  the  modulation  by 
this  method  is  8%.  Predictive  ability  appears  to  be  degraded 
somewhat  when  the  oceanic  fluctuations  are  of  the  same 
order  as  the  bandwidth  of  the  modulation,  i.e.,  when 
v  %  A  the  average  error  was  17%.  It  should  be  remembered, 
however,  that  the  modulation  is  not  exactly  uniform  and, 
therefore,  using  Equation  (1.176)  is  an  approximation 
to  begin  with,  and  this  approximation  is  worst  when  v  %  A. 
Finally,  I  note  that  the  error  is  consistently  negative, 
i.e.,  the  foregoing  method  undarpredicts .  The  reason 
for  this  is  unknown. 

Performance  on  the  Chi-square  test  was  also  quite 
good.  All  88  runs  passed  at  the  a-  .10  level  of 
significance  when  Equation  (1.7)  is  fitted  to  the 
histograms  of  A.  In  fitting  Equation  (1.176)  to  the 
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Fig. 70  The  error  in  percent  versus  number  of  runs  in  determining 
the  bandwidth  of  the  modulation  of  the  CASE  data. 

The  shaded  portions  represent  those  runs  in  which  i/»A. 
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histograms  of  $  for  76  runs  with  modulation,  only  2  runs 
fail.  For  the  12  runs  without  modulation,  only  1  run 
fails  when  Equation  (A8)  is  fit  to  the  histograms  of  ?. 

The  fit  of  Equation  (A8)  to  the  76  modulated  runs  is  as 
anticipated  extremely  poor;  29  of  the  76  runs  failed 
and  those  that  passed  did  so  by  a  very  small  margin. 

In  Table  XVII,  I  have  compiled  the  values  of  v 
(from  A) ,  A,  actual  A,  range,  range  rate,  source 
configuration,  frequency,  and  receiver  number  for  each  of 
the  76  modulated  runs.  Under  source  configuration  C 
indicates  fixed  on  the  seamount,  T  indicated  towed  by  a 
surface  ship,  and  S  indicates  mounted  on  a  submersible. 

Figures  71  and  72  are  representative  of  the 

modulated  runs  when  A  >>  v.  Figure  71  shows  the  time 

« 

series  and  histogram  of  A  and  the  best  fit  of 
Equation  (1.7) .  Figure  72  includes  the  time  series  and 
histogram  of  <j>  and  the  best  fit  of  Equation  (A8)  and 
Equation  (1.176).  Reference  (33]  contains  the  plots  for 
all  the  runs  as  well  as  the  FORTRAN  listing  of  the  data 
analysis  program  for  the  CASE  experiment. 

Comparison  of  v2  from  the  CASE  experiment  with  the 
Dyson,  Munk,  and  Zetler  model  [11]  is  unfavorable.  The 
values  of  v2  obtained  at  15  Hz  are,  in  fact,  greater  than 
those  obtained  at  33  Hz,  contrary  to  the  frequency  scaling 
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<£  RAD/  SEC 


Fig.  72  The  time  series  and  histogram  of  for  the  same 
run  shown  in  Fig.  71  is  plotted  along  with  the  best 
fit  of  equation  (A8)  which  assumes  no  modulation 
(a)  and  the  best  fit  of  equation  (1.176)  which 
includes  the  effect  of  source  induced  modulation  (b). 
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proposed  in  their  model,  and  these  values  are  in  turn 
several  orders  of  magnitude  greater  than  the  values 
predicted  by  the  model  (see  Table  XVIII  and  Figure  73) . 
However,  the  rms  single  path  phase  fluctuation  for  the 
frequencies  and  ranges  of  the  CASE  experiment  as  predicted 
by  an  internal  wave  model  only,  appear  to  be  <<  2t,  and  if 
so,  the  model  would  not  be  applicable  to  the  CASE  data. 

In  order  to  obtain  a  very  rough  estimate  of  applicability. 
Equation  (118)  (for  small  3),  Reference  [34],  has  been 
plotted  in  Figure  7  4  for  <<+>2>  =  (2tt)2,  for  range  in  km 
vs.  frequency  in  Hz.  That  portion  of  the  figure  above 
and  to  the  right  of  the  line  indicates  that  phase 
fluctuations  due  to  internal  waves  are  >  2v.  The  dotted 
lines  represent  the  CASE  experiment,  the  circles 
represent  the  Eleuthera  data.  Despite  the  assumptions 
involved,  the  good  results  of  DMZ  with  the  Eleuthera 
data  (Table  XI) ,  and  the  poorer  results  with  CASS,  could 
be  explained  by  Figure  74.  It  is  apparent  that  some  other 
mechanism  must  account  for  the  fully  saturated  phase 
random  process  evident  in  the  CASE  experiment.  Other 
possible  mechanisms  could  include  range  rate,  other  ocean 
dynamic  phenomena  such  as  Rossbv  waves  or  meso-scale 
eddies,  tidal  currents  magnified  by  bottom  interaction 
at  the  receiver  location,  or  rough  scattering  effects. 
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TA3LE  XVIII 


Comparison  of  the  average  values  of  v  obtained  from 
the  CASE  experiment  to  the  predictions  of  the  DMZ 
internal  wave  model.  Also  shown  are  the  variances  for 
the  values  of  v  obtained  from  CASE . 


v (mHz) 

Variance  of 

■Mm 

Case  v 

ESI 

CASE 

(mHz) 

Free. 

15 

o 

CO 

38* 

53 

35*  55 

(Hz) 

33 

t-4 

CO 

19+ 

20 

14+  15 

*These  values  were  computed  leaving  out  runs  71-73  and  113-115 
that  exhibited  significant  ncn-zero  means  in  the  quadrature 
components . 

+ These  values  were  computed  leaving  out  runs  6  6-6  3  that 
exhibited  significant  non-zero  means  in  the  quadrature 
components . 
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f  Hz 


.  2  1/2 
Fig. 74  Range  vs.  Frequency  for  <  <p  >  =  2ir.  The  portion 

of  the  graph  above  and  to  the  right  of  the  line  is  the 

fully  saturated  region  based  on  internal  wave  models. 

The  CASE  experiment  is  indicated  by  the  dotted 

lines  and  the  Eleuthera  data  by  the  circles. 
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Although  no  strong  correlation  was  found  between  v 
and  range  rate,  it  is  apparent  from.  Figure  75  that  the 
larger  values  of  v  were  obtained  at  0  range  rate.  I  also 
investigated  other  correlations.  I  found  none  between 
v  and  range,  or  time  of  day  (see  Figures  76  and  77 
respectively) .  Because  the  CASE  experiment  was  run  over 
approximately  one  and  one-half  months  with  experimental 
runs  occurring  on  many  different  days,  tidal  effects 
cannot  be  inferred  from  the  times.  At  the  present  time, 
the  dates  of  each  of  the  runs  are  unavailable,  •  making 
tidal  checks  impossible. 

In  Figures  78  and  79  I  have  plotted  the  values  of 
v  vs.  receiver  for  15  Hz  and  33  Hz  respectively.  Because 
each  receiver  was  monitoring  a  given  experimental  run 
simultaneously  I  have  connected  the  values  of  v  at  each 
receiver  by  lines  for  each  run.  The  time  of  day  is  noted 
to  the  left  of  the  receiver  1  values.  An  "N"  next  to 
the  time  indicates  that  there  was  no  range  rate,  either 
the  source  was  fixed  on  the  seamount  or  motionless  for 
that  run.  For  15  Hz  (Figure  78) ,  receiver  2  consistently 
sees  a  higher  value  for  v  than  receiver  1,  and  receiver  3 
consistently  sees  a  larger  value  than  receiver  1.  Between 
receivers  2  and  3  the  results  are  mixed.  For  33  Hz 
(Figure  79) ,  receiver  2  consistently  sees  a  v  greater  than 


rad/sec 


Time  (ZULU) 

Fig.  77  A  plot  of  v  vs.  time  of  day.  The 
circles  represent  15  Hz,  the  x’s 
Hz. 
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Fig.  79  Plot  of  v  vs.  receiver  for  33  Hz.  The  lines 

connect  values  of  v  obtained  simultaneously  at 
each  receiver.  The  start  times  for  each  run  (ZULU) 
are  indicated  to  the  left  of  the  receiver  I  values. 
The  N  indicates  a  zero  range  rate  run. 
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receiver  3.  For  8  of  the  12  runs,  smaller  values  of  v 
are  seen  at  receiver  3  than  are  seen  at  receiver  1. 

Between  1  and  2  the  results  are  mixed.  These  results 
suggest  the  possibility  of  rough  scattering  or  a  tidal 
current  phenomena  which  would  be  receiver  dependent. 

Note,  however,  that  the  no  range  rate  runs  ("N")  for 
both  frequencies  all  exhibit  the  same  pattern,  with 
low  values  at  receiver  1,  the  highest  values  at  receiver  2, 
and  low  values  again  at  receiver  3.  For  all  receivers 
range  rate  groups  v  at  lower  values  than  the  no  range 
rate  runs,  the  latter  having  a  consistently  higher  mean. 
Thus,  as  observed  earlier  (Figure  75) ,  there  does  appear 
to  be  some  small  correlation  with  range  rate  to  the  extent 
that  either  there  is  range  rate  in  which  case  v  is 
independent  of  the  amount,  or  there  isn't,  in  which  case 
v  exhibits  a  higher  mean  and  variance.  This  is  more 
graphically  illustrated  in  Figures  80-85  in  which  I  have 
plotted  v  vs.  range  rate  now  separating  the  values  by 
receiver  as  well  as  frequency. 
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Range  Rate  kts 

Fig.  82  A  plot  of  v  vs.  range  rate  for  the  33  Hz  signal 
measured  at  receiver  3. 
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Fig.  83  A  plot  of  v  vs.  range  rate  for  the  15  Hz 
signal  measured  at  receiver  I. 


Fig.  84  A  plot  of  v  vs.  range  rate  for  the  15  Hz 
signal  measured  at  receiver  2. 
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Fig.  85  A  plot  of  v  vs.  range  rate  for  the  15  Hz 
signal  measured  at  receiver  3. 


DISCUSSION  AND  CONCLUSIONS 


The  major  contributions  of  this  thesis  are  twofold. 
First,  on  the  theoretical  front  the  understanding  of 
the  statistics  of  acoustic  signals  propagated  to  long 
ranges  in  the  ocean  has  been  significantly  advanced  by 
the  derivation  of  many  new  pdf's,  particularly  for  the 
rate  variables  and  joint  pdf’s  of  amplitude  and  amplitude 
rate,  as  well  as  the  theoretical  development  of  the  effects 


of  modulation.  Second,  the  application  of  this  new 
theory  to  ocean  acoustic  data  has  revealed  new 
understanding  of  the  effects  of  finite  bar.dwidths  and/cr 
modulation  on  the  statistics,  as  well  as  a  clearer  picture 
of  the  limitations  of  current  models  relating  v2  to  ocean 
or  experimental  phenomena. 

To  be  more  specific,  the  completion  of  the  family  of 
pdf's  for  the  single  source  and  particularly  the  derivation 
of  ?^(A)  has  permitted  direct  measure  of  v2  from  amplitude 
quantities  which  are  independent  of  bandwidth  and/or 
modulation  effects.  However,  more  research  is  needed  to 
understand  the  effects  of  non- stationarities  in  on  the 
amplitude  and  amplitude  rate  pdf's.  Although  c£  does  not 
appear  explicitly  in  (A) ,  non-stationarities  in  c2  will 
possibly  have  some  small  effect.  The  effect,  however,  will 
be  related  to  the  rate  of  change  of  a ^  and  not  to  the 
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absolute  change  which  would  affect  those  variables  in 
which  a2  appears  explicitly,  so  that  measurement  of  v2  via 
A  is  the  preferred  method. 

The  Eleuthera  data  analysis  reveals  that  for 
measuring  v2  use  of  the  time  series  and  histograms  of  $ 
and  Equation  (A8)  apply  only  when  B<<2v,  A  < <  v ,  and  the 
Doppler  shift  due  to  relative  source/receiver  motion  is 
<<  2v.  It  is  true  that  relative  scurce/receiver  motion 
induces  Doppler  modulation,  which  is  path  dependent  as 
each  path  has  a  different  arrival  angle  at  the  receiver. 
However,  in  forming  A,  or,  in  fact,  any  of  the  amplitude 
or  amplitude  rate  variables,  the  mean  Doppler  is  removed 
and  it  is  only  the  path  to  path  differential  Doppler  that 
remains.  For  $,  however,  the  mean  Doppler  as  well  as  the 
differential  Doppler  contribute,  and  naive  use  of 
Equation  (A8)  on  <p  will  result  in  error.  Even  when  the 
mean  phase-rate  is  removed,  as  was  done  with  the  Eleuthera 
data,  we  discover  that  a  new  criterion  must  be  met,  namely 
that  the  error  in  removing  the  mean  phase  rate  must  be 
<<  2v . 

Recognition  of  the  A  method  of  measuring  v2  is  alone 
one  of  the  major  contributions  of  this  research.  However, 
a  more  thorough  analysis  is  in  order  to  determine  the 
effect  of  the  differential  Doppler  on  the  pdf's  of 
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amplitude,  amplitude  rate,  and  phase  rate  when  there  is 

relative  source/receiver  motion.  Because  arrival  angles 

are  small  for  long-range  propagation,  it  is  anticipated 

that  this  effect  will  be  small  if  not  negligible  in  most 

cases.  For  the  Eleuthera  data,  the  mean  drift  rate  along 

the  transmission  path  was  on  the  average  100  m/hr  [12] . 

Assuming  a  maximum  arrival  angle  of  14°  the  maximum 

-4 

differential  Doppler  shift  is  1  x  10  Hz  at  220Hz  and 
-4 

2  x  10  Hz  at  406Hz,  an  order  of  magnitude  less  than  v 
for  these  data.  For  the  CASE  data,  the  differential 
Doppler  was  also  negligible  for  all  the  runs  analyzed. 

In  addition  to  range  rate,  other  mechanisms  which 
affect  the  value  of  vz  must  be  researched.  Analysis  of 
the  CASE  data  revealed  that  some  mechanism  other  than 
internal  waves  must  account  for  the  fully  saturated 
phase  random  nature  of  the  data.  My  preliminary 
correlations  indicate  that  a  receiver  related  phenomenon 
such  as  slope- inf luenced  tidal  currents,  or  rough 
scattering,  would  be  a  good  place  to  start.  However,  other 
mechanisms  including  Rossby  waves,  meso-scale  eddies,  etc. 
cannot  be  counted  out.  I  must  confess  to  some  feeling  of 
uneasiness  surrounding  the  great  run  to  run  variance  of 
the  measured  values  of  v2  in  the  CASE  experiment.  Clearly 
a  mechanism  of  some  temporally  varying  nature  must  be 
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accountable.  This  again  points  toward  a  tidal  phenomena. 
The  Eleuthera  data  reveal  a  temporally  stable  v2,  even 
more  so  than  Hamblen  [9]  discovered,  because  his 
analysis  includes  what  is  now  recognized  as  a  Doppler 
error. 

The  analysis  of  effects  of  finite  bandwidth  and/or 
modulated  signals  coupled  with  the  new  confidence  in 
obtaining  v2  minus  these  effects  has  revealed  a 
potentially  powerful  technique  for  separating  source  and 
ocean  effects  in  the  received  signal.  This  technique 
uncovered  the  effect  of  the  Doppler  error  in  the  Eleuthera 
data  with  great  improvement  in  crossing  rate  predictions, 
and  predicted  the  bandwidth  of  the  modulation  of  the  CASE 
data  with  an  average  error  of  8%.  This  technique  could 
be  refined  by  use  of  more  sensitive  statistical  tests 
such  as  the  Kolmogorov-Smirnov  test.  Also  the  pdf's  for 
many  different  kinds  of  modulation  could  be  tried,  not 
only  to  determine  such  parameters  of  the  modulation  as 
its  bandwidth,  but  also  the  nature  of  the  modulation 
itself  when  it  is  unknown,  by  comparing  the  performance 
of  the  selected  densities  on  the  histograms. 

Aside  from  the  applications  of  the  amplitude 
statistics  to  distant  shipping  noise  problems  [4,3]  the 
analysis  of  the  statistics  of  multiple  source  cases  is  a 


relatively  untapped  reservoir  especially  for  the  rate 
statistics.  In  order  to  improve  predictions  in  the  tails 
of  the  densities  whose  exact  solutions  remain  unknown, 
the  Chernov  bound  or  "tilted"  density  [25]  approach 
mentioned  earlier  should  be  investigated.  I  believe, 
however,  that  the  analysis  presented  is  complete  enough, 
and  the  problems  remaining  few  enough  to  warrant 
investigation  into  the  applications  of  this  analysis.  The 
coherent  source  analysis  may  be  a  reasonable  model  for 
determining  the  statistics  for  the  signal  plus  ocean 
noise.  The  multiple  source  cases  analyzed  may  well  be 
applicable  to  other  noise  problems.  Table  III  provides 
a  comprehensive  summary  of  the  state  of  this  analysis 
at  the  present  time. 

The  analysis  contained  in  Appendix  C,  relying  on 

some  results  of  earlier  investigators  in  fields  other 

than  acoustics,  completes  the  solution  of  the  statistics 

for  the  amplitude  for  N  <_  3  paths.  For  N  _>  4  I  have 

shown  with  the  computer  simulation  that  limiting  pdf's 

suffice.  More  research  is  needed  into  the  exact  nature  of 

Pi  (8  ),  the  pdf  for  the  single  path  phase.  As 
9n  n 

previously  noted  by  Hamblen  [9]  and  as  supported  by  the 
analysis  in  Appendix  C,  the  rate  densities  of  amplitude  and 
multipath  phase  for  M  small  cannot  be  found  unless  (8  ) 


-279 


is  known.  Mo  investigators  to  date  have  related,  to  my 

knowledge,  any  knowledge,  experimental  or  theoretical, 

about  PI  (9  )  aside  from  the  usual  Gaussian  or  uniform 
9n  n 

assumptions  as  I  have  done.  In  Appendix  C  I  have 
solved  for  P.(x)  when  W  *  2  paths  making  these 

A 

assumptions . 

In  summary,  the  theoretical  analysis  of  long  range 
acoustic  propagation  presented  in  this  thesis  has  been 

ivi4v\  cJclHc  -$rotn  Oco oi-Sc  &x pe ri <v\em ts 

supported  by  both  computer  simulation  and  comparison ^with 
extremely  satisfactory  results.  I  hope  that  the  analysis 
presented  and  the  conclusions  reached  will  be  of  use  to 
other  researchers  in  the  future. 
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APPENDIX  A 

FIRST  ORDER  DENSITIES  FOR  THE 
SINGLE  NARROW3AND  SOURCE 


Assembled  below  are  the  pdf's  associated  with  a  single 
narrowband  source  and  their  characteristic  functions  D-41  . 

Notation: 


P  (x)  is  the  pdf  of  random  variable  x 

M  (u)  is  the  characteristic  function  of  x 
X 

U  is  the  expected  value  or  mean  of  x 
X 

a2  is  the  variance  of  x 


10  log1Qe  *  4.34... 


A.  Amplitude  Densities 


a)  v»>  ■  ^;exp(- 


]  ,  D  >  0  ;  Rayleigh 


(Al) 


up  =  a^/1  ,  a2  =  a2{2  -  ir/2) 


oj2a? 

Mp(u)  =  exp (-  . 4  -)  D_2(-iuja1) 


where  Dn(z)  is  the  parabolic  cylinder  function  [2  4]. 


(2) 


px(x)  = 


-  exp  [- 

2al 


-3  ,  x  >  0 


Exponential 


y 


x 


4a 


1 


M  (w)  ■  - = - 

X  1  -  i2a^ 


(3)  P .  ( A)  »  — —  exp  [—  -  -i-  exp  (^)  ]  ,  -®<A<®; 
A  2ea£  £  2a£  E 


Lag-Rayleigh 


UA  =  etln  2a^  -  y] 


Ma(u)  =  (2a|)ieur(l  +  ieu) 


where  y  ■  Euler's  constant  =  .5772... 
T(z)  is  the  gamma  function  £241. 
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1/2it  0  <  3<2:t 


otherwise 


;  Uniform 


( A4 ) 


2  IT 

u4>  "  77  '  a4>  “  3 


M^(aj)  *  —  exp  (-iuj-rr)  sinuTT 


B .  Amplitude  Rate  Densities 


?.(p)  a  - - -  exp  [-  — - - ]  /  -®<p<“  ?  Gaussian 

p  /2io*v2“  ^~2"2 


2  a£v‘ 


(A5) 


ui  -  0  '  as  ■ 


(ai)  *  exp(-  ~  ui2a2v2) 
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(6)  P 


M 


(7)  P 


M 


where 


1  I  v  1 

.  (x)  s  -  exp[-  — -H  ,  -oo<x<«  ;  Laplace  (A6) 

x  4va2  2 


y .  51  0 ,  a?  *  8ciJ\)2 
X  X  1 


U)  =  - i - 

*  1  +  4w2a£v2 


(A) 


l/2v 


2e  [1  + 


A2 


4e2v2 


-1 


3/2 


-co<,\<cn  .  Longuet-Higgins 


(A7) 


y  *  =  0  , 


•^(aj)  =  2ev  |  aj  |K1  (2sv  |  w  j  } 


K^(z)  is  the  modified  Bessel  function  of  order  one. 
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(8)  For  no  modulation  only: 


*:(♦> 


i£> _ _ 

:  *  T72 


2U  +  $rl 


»<$<®  ;  Longuet-Higgins  (A8) 


0  '  «J 


M*  (uj)  =  v|uj|K^(v|<d|  ) 
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APPENDIX  B 

THE  ERLANG  AND  RELATED  PDF's 

The  Erlang  pdf  (Equation  1.19,  Section  1.2. 1.1)  was 
derived  by  A.  K.  Erlang  in  the  early  1900's  in  connection 
with  waiting  times  in  telephone  operations.  It  is  the 
probability  of  the  time  until  the  n  arrival  in  a 
Poisson  process,  or  the  density  of  the  nth  order  inter¬ 
arrival  time  [35].  Note  that  the  first  order 
interarrival  time  (L*l)  is  an  exponential  pdf.  L  can  be 
extended  to  include  noninteger  values  by  merely  replacing 
the  factorial  of  the  Erlang  with  the  T  function  and  hence 
the  name  gamma  pdf.  Although  mathematically  valid,  the 
simple  physical  interpretation  of  the  Er Icing  breaks  down 
when  noninteger  values  of  L  are  introduced.  In  the  case 
of  distant-shipping  noise,  L  is  an  integer. 

In  the  language  of  queing  theory,  it  is  interesting 
to  note  that  determination  of  the  pdf  of  x  for  Case  (c) 
of  Section  1.2.1  in  all  its  generality  is  the  same  as  the 
solution  for  the  pdf  of  the  nth  order  interarrival  time 
in  a  renewal  process  in  which  the  first  order  interarrival 
times  are  independent  random  variables  distributed 
according  to  Equation  (1.19)  with  different  and  L^ . 

Though  originally  derived  independently,  the 
slightly  less  general  chi-square  pdf  with  n  degrees  of 
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freedom  can  be  derived  from  Equation  (1.19)  by  merely 

letting  a  -  1/2  and  L  «  n/2  CL  4  . 

If  we  let  a  ■  m/ux  and  L  »  m,  and  make  the 
1/2 

transformation  R  *  x  in  Equation  (1.19),  we  have  the 
density  of  the  short  time  rms  pressure: 


Pr(R) 


2mmR2m-1  e-(m/ux)R2 

r<a)u“ 


(Bl) 


In  terms  of  the  statistics  for  x»  »  »  (u  /a  )2,  which  is 
the  inverse  of  the  normalized  variance  of  x* 

Equation  (Al)  was  first  proposed  by  M.  Nakagami  in 
1943  [ 7 ]  to  describe  the  envelope  of  long  range  h.f. 
radio  wave  propagation  undergoing  rapid  fading  and  is 
known  as  the  "m"  distribution.  When  m»l,  Equation  (Bl) 
is  the  Rayleigh  distribution  as  expected.  One  can  note 
that  the  phase-random  model  of  long  range  acoustic 
propagation  is  analogous  to  the  rapid'  fading  of  long  range 
h.f.  propagation.  Thus,  the  connection  between  the  ”m” 
and  Erlang,  or  gamma  pdf,  is  more  than  just  a  functional 
similarity;  they  describe  in  alternative  language  the  same 
process . 
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APPENDIX  C 

THE  STATISTICS  FOR  M  <  3  PATHS 


The  solution  for  the  pdf  of  the  amplitude  of  a  vector 
which  is  the  sum  of  many  vectors  added  with  random  phases 
has  been  a  problem  of  long  standing  interest.  Lord 
Rayleigh  was  probably  the  first  to  investigate  this 
problem  in  two  papers,  the  first  published  in  1880  [1] , 
and  the  second  in  1899  [36].  He  was,  however,  concerned 
with  the  limit  when-  the  number  of  summed  vectors  is  large, 
and  he  derived  the  density  of  the  amplitude  for  this 
limiting  case  which  bears  his  name.  The  first  investigators 
to  tackle  the  problem  of  small  N  (note  that  the 
characterization  of  summed  paths  as  random  vectors  is 
mathematically  identical)  was  Kluyver  (1905)  [37],  and 
Pearson  (1906)  [38] ,  the  latter  in  connection  with  mosquito 
migration! 

For  arbitrary  A!  and  unequal  amplitudes,  A^,  Kluyver 
[37]  obtained  the  general  solution  in  integral  form  as 
follows : 


A/Pp<P) 


N 

J  (up)  II  J  (uA.  )du 

W  4—1  °  A 


(Cl) 


For  the  problem  of  sums  of  independent  multipaths,  we  have 
that  A^  ■  r  for  all  i  and  thus: 
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w Pp(p)  ■  P  |  u  JQ(UP)  [JQ(ur)  ]N  du  (C2) 

o 

Pearson  [38]  obtains  solutions  of  this  equation  for  N*2,3 
and  for  N^4  he  obtains  solutions  in  terms  of  series  of 
Bessel  functions. 

C.l  Solution  for  M  »  2 

Though  I  derived  this  solution  independently,  it  was 
Pearson  [38]  who  obtained  it  first.  For  two 

paths  I  have  for  the  quadrature  components: 

X  »  r(cos0.  +  cos0,) 

1  *  (C3 ) 

Y  »  r(sin©1  +  sin©2) 

where  0^-and  ©2  are  independent  random  variables  distributed 
uniformly  between  0  and  27r.  Forming  the  sta  mean  square 
pressure  x  from  Equations  (C3), 

2  2  2.2 
X  *  r  (cos  0^  +  2cos0^cos02  +  cos  ©2  +  sin  0^ 

+  2sin0^sin02  +  sin202)  (C4) 

Using  some  trigonometric  identities  I  have 

X  -  2r2  +  2r2cos(91  -  02)  (C5) 

Taking  advantage  of  the  symmetry  of  the  cosine  the  random 
variable  defined  by  9^-62  behaves  as  it  were  uniformly 
distributed  between  0  and  tr.  Thus,  forming  the  cumulative 
distribution  function: 


-  *  ^ 
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2"  <  X 


<  *<x)  *? 


1  f 


___-l  /X-2rv 
cos  (—  «— ) 

2rz 


<29,  0  <  9  <  it 


and 


2V*>  "  -Sx  2P<X(X) 


I  obtain  the  final  solution: 


2PX<X)  "  w 


0  <  X  <  4r 


where , 


ux  ■  2r2,  °x  ■  2r4 


Making  the  transformation  to  p  and  A: 


,P  (p)  -  r  .  1-  r  ' 


77  /4r^-p2 


0  <  p  <  2r 


Up  -  4r/ir,  ap2  -  r2(2  -  16/ir2) 


and,  A/2e 

2PAU>  ■  FT  /4rAeA 7? 


,  -®  <  A  <  2eJln2r 


UA  -  2einr,  c2  ■  e2[ir2/6  -  2(inr)2] 


From  Equation  (CS)  I  have 


X  ■  2r^(92  -  8^  sinOj^  -  92) 


(C6) 


(C7) 


(C8) 


(C9) 


(CIO) 


(Cll) 


•  I 


/ 
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It  is  clear  from  Equation  (Cll)  that  for  small  N  exact  knowledge 
of  P§  (9  )  is  required  before  rate  densities,  or  joint 
densities  of  amplitude  and  amplitude  rate  can  be  found. 

I  will  now  perform  the  calculations  required  to 
find  2Pv(x)  as  stating  Pq  (6^)  is  Gaussian  and  then  assuming 

A  •• 


Pa  ( 0_)  is  uniform.  For  the  former: 
n  .  - 

9n2 


<9J 

9n  n 


2v‘ 


/  2irv2 


l9n 


<  » 


(C12) 


The  pdf  for  t  •  sin  (9^  -  92)  is  given  by 


Pt(t) 


1 1 1  <  1 


The  pdf  for  z 
PZU) 


92  "  9i  is 


inspection 

f  |  z|  <  • 


(C13) 


(C14) 


Because  z  and  t  are  independent  random  variables  the  joint 
density  of  z  and  t  is  the  product  of  Eqrts.  (C13)  and  (C14)  . 
I  seek  the  pdf  for  the  product  u  ■  zt,  and  explorting 
symmetry, 


Pu(u) 


u/t 

o 


P 


2ft 


( Z ,  t) 


dzdt 


(C15) 
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Differentiating  under  the  integral  sign  I  obtain 

top  (-  -£-> 

■A-  - d« 

v/if  Jo  td-tV^ 


Pu(u) 


(C16) 


irvi 


I  now  make  the  change  of  variable  a  u  and  then  integrate 

.  2 

and  make  the  final  change  of  variables  x  a  2r  to  obtain: 


f  *2 

f  *2 

x 

Ko 

X 

— 2 - -  exp 

4r  ttv/tt 

.  32r4v2j 

t32r4v* 

<  « 


(C17) 


where  KQ(z)  is  the  Modified  Bessel  function  of  order  0. 


For  9n  uniformly  distributed 


P5  (6_) 
9n  n 


2/3v 

0 


elsewhere 


(CIS) 


r  *  2  1  2 

Note  the  factors  in  Eqn.  (C18)  insure  E[9n  J  “  v  as  required. 
Now  the  density  for  z-^-S^is 

,  |  z  {  <2  /3  v 


P,(«) 


1 

2 /Tv 

1° 


i-iaH 

2/3v 


(C19 ) 


elsewhere 


The  pdf  for  t  »  sinO^  -  92)  is  given  by  Eqn.  (C13)  . 
Defining  u  *  zt  I  integrate  over  the  joint  density  as 
before.  The  integration  is  non  trivial  and  requires 


|  - 


/ 
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careful  concern  for  the  limits : 

fi  rU/t 

Pu(u)  *  h  |  PZ/t(2,t)dzdt 

5  Jo 

u 

2/3v 

■u/2/Tv  f2/lv 

*  pz,t(2't)  «*at  (c20> 

O  Jo 

In  each  case  I  perform  the  integration  over  z  first  and  then 

differentiate  under  the  remaining  integral,  and  integrate. 

.  2 

After  making  the  final  change  of  variables  x  a  2r  u  I  have: 


C.2 


Solution  for  N  •  3 

For  this  case  following  Pearson  [3?], 
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and  lastly. 


3P.(A) 


I  ,3A, 

eXP(4?> 


- 5—377  yR  (j»  Y) /  -»  <  A  <  2s£nr 

2eir  r  7 

,3A, 

axP(4F>  „  , 

- y  3  yy  K  (7,  7)  ,  2eJ.nr  <  A  <  2eZn3r 

2eir  rJ/^  *  Y 


,  elsewhere 


(C24) 


where 


16r3  exp  (jj) 


[exp(^0  +  r]3[3r  -  exp(~)] 


•  •  •  9  • 

The  rate  density  for  x>  P,  or  A  assuming  a  Pg  (0  ) , 

n  n 

are  sufficiently  more  complicated  than  *V  =  2  that  I  have 
not  solved  for  them. 

C.3  Computer  Simulation 

Figures  C1-C9  are  histograms  of  the  computer 
generated  samples  of  the  amplitude  variables  for  H  «  2,3, 
and  4  paths.  The  exact  pdf  for  U  *  2,  and  3  is  shown  with 
the  limiting  pdf.  Note  that  for  .V  =  4  the  performance  of 
the  limiting  pdf  is  good  enough  (see  Table  V  for  Chi-s!quare 
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test  results)  to  warrant  use  of  the  limiting  forms  for 

M  >  4 .  In  Figure  CIO,  I  have  plotted  Equation  (C21)  with 

the  histogram  of  x  for  N  =  2  (recall  that  the  computer 

simulation  assumes  Pi  (0  )  is  uniform) ,  as  well  as  the 

9n  n 

limiting  pdf.  Equation  (1.4) . 


computer  generated  histogram  for 


The  computer  generated  histogram  for  N  s  2  paths,  plotted 
with  the  exact  pdf,  equ.  (C9),  and  the  limiting  pdf  (N-*»a>), 
equ.  (Al). 


8  0 


computer  generated  for  histogram  for  N  *  3  paths, 
ed  with  the  exact  pdf,  equ.  (C22),  and  the  limiting 
(N-^co),  equ.  < Al). 


computer  generated  histogram  for  N  =4  paths,  plotted 
the  limiting  pdf,  equ.  (Al). 


ompu ter  generated  histogram  for  N  =  2  paths,  plotted 
the  exact  pdf,  equ.  (CIO),  and  the  limiting  pdf 
to),  equ.  (A3). 


The  computer  generated  histogram  for  N  *  3  paths,  plotted 
with  the  exact  pdf,  equ.  (C24),  and  the  limiting  pdf  (N-^a 
equ.  ( A3). 


computer  generated  histogram  for  N  *  4  paths,  plotted 
limiting  pdf,  equ.  (A3). 


Computer  generated  histogram  for  N  =  2  paths,  plotted  with 
the  exact  pdf,  equ.  (C2I),  and  the  limiting  pdf  (N-^©), 
equ.  (1.4). 
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APPENDIX  D 

AMPLITUDE  PARAMETER  VARIATION 

In  this  appendix  I  will  briefly  introduce  the  analysis 
involved  in  amplitude  parameter  variation.  This  analysis 
is  the  first  step  towards  understanding  the  statistics 
of  received  signals  that  are  either  nonstationary  in  r 
or  a£,  or  when  r  or  is  purposely  modulated  at  the 
source. 

For  the  first  case,  I  consider  that  r  is  a  random 
variable  distributed  uniformly  between  a  lower  limit  of 
b  _>  0,  and  an  upper  limit  of  a  >  b.  It  is  evident  from 
Equation  (2.1c)  with  L*1  that 

PAV  »  rp  (Dl) 

where  now  p  can  be  considered  a  random  variable,  distributed 
Rayleigh  according  to  Equation  (Al)  but  in  which 
»  N/2.  Note,  however,  that  there  is  an  implied 
assumption  that  N  ^  4.  From  Equation  .(Dl)  and  the  fact 
that  r  is  a  measure  of  the  gain  of  the  signal  which  can  be 
controlled  at  the  source,  or  of  the  change  in  pathwise 
signal  strength,  it  is  true  that  r  and  p  are  independent 
random  variables.  This  assumes  of  course  that  during  the 
experiment  the  number  of  paths  remains  constant.  Thus, 
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forming  the  joint  density  of  r  and  p  I  obtain  the 
cumulative  distribution  function. 


PAV  PAV 


P<PAV(pAV> 


(a 


-i—  f  f  °  P  exp(-  ^Idrdp 

-b)0l  p...  v  1 


b 
a 


AV 


r  a 


(a  -  b) 


a 

|  p  exp  (-  -^yjdrdp 


(D2) 


Performing  the  integration  and 
I  obtain  the  final  result: 


taking 


7p 


AV 


■<‘>av<!’av> 


“V 


1 

(a  -  b) 


ir 

7 


Erf  (■ 


AV 


bc^/7 


Erf  ( 


PAV 

aa^/2 


(D3) 


(a  +  b) 


a+b  /"S  ir  8ab 
"T"  aiy  3  ”7"  3  (a+b)  2 


Note  that  when  a  ■ 
variation  values 


b  *  r,  u  and  a 
PAV 

(see  Equation  A1 


assume  their  no 

PAV 

,  Appendix  A) ,  and 
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likewise. 


lim  Po  (pAV)"Po(p) 
a-**r  PAV  AV  p 

b-*r 


With  the  appropriate  transformations,  I  get: 


pxav(xav> 


(xAV) 

2a^  (a  -  b) 


7T 

2 


^XAV  *  XAV 

Erf  ( — ~)  -  Erf  ( — — ) 


bo1/2 


acr^/2 


(D4) 


PA  (AAV} 
lAV  AV 


AV 

i  exp(Tr}  ft 
e  2a1(a-b)  "J 


Eri 


AAV  - 
exp(-^) 


ubaL/I 


-  Erf 


uav  - 

exp(^r] 


i  ac^/2  ^ 


(D5) 


Note:  »  M/2  in  Equations  (D2-D5) . 

In  order  to  compare  the  effect  of  certain  levels  of 
uncertainty  in  r  on  the  pdf'  for  p,  the  following 
procedure  must  be  employed.  Given  that  r  is  uniformly 
distributed  between  b  and  a,  its  average  value  is  a+b/2. 
Plot  Equation  (Al)  with  cr£  »  ^  (-y^)  and  then  plot 
Equation  (D3) .  If  this  is  not  done,  then  one  may 


l 


r 
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erroneously  conclude  that  the  effect  of  uncertainty  in 
r  is  quite  dramatic  when  in  fact  it  is  not! 

I  have  taken  various  values  of  a  and  b  such  that 
*  1  and  plotted  Equation  (Al)  with  Equation  (D3)  in 
Figure  Dl.  Note  that  with  50%  (3dB)  uncertainty  in  r, 
the  variation  from  Equation  (Al)  is  not  terribly 
significant;  however,  for  uncertainty  >  50%,  the  effect 
becomes  very  noticeable  in  the  pdf.  In  Figure  D2  I  have 
plotted  Equation  (A3)  with  Equation  (D5) .  As  before 
with  3dB  of  uncertainty,  the  effect  is  small.  However, 
for  ^  7dB  or  greater  uncertainty,  the  effect  is  large  and 
in  fact  is  very  close  to  the  limiting  form  (100%  or 
°°dB  uncertainty) .  Nakagami  [7]  obtains  a  similar  result 
in  relation  to  uncertainty  of  the  mean  intensity  in 
dB  of  multipath  BF  propagation.  In  obtaining  the 
limiting  pdf  (i.e.,  100%  uncertainty),  it  is  necessary  to 
let  b=0  and  a  **  2r  in  Equations  (D2-D5)  (note,  Erf(°°)  =  1) 
and  rewrite  the  equations  analytically  first,  otherwise 
considerable  computer  time  and/or  overflows  will  result. 

Next,  I  assume  r  varies  sinusoidally  about  a  mean 
given  by  a  and  an  amplitude  of  b.  Thus, 


Pr(r) 


1  _ 1 _ 

*  /S2  -  (r  -  a)  2~ 


|  r— a |  <  b 


i 


p,<?> 


0  I  2  3  4  5  6  7 

P 


Fig.  Dl  Comparison  plot  of  the  effect  of  amplitude 

variation  on  the  pdf  for  p.  Equ  (Al)  is  plotted 
with  equ  (03)  with  various  levels  of  uncertainty. 


with  various  values  of  uncertainly. 
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The  cumulative  distribution  function  is  given  by 


(PAV) 


/bz  -  (r-a) 


0  * 

exp(-  ^j)drdp 


_ p _ 

>/b2  -  [r  -  a T2“ 


exp( 


•)drdp 


.  .  . (D6) 

Integrating/  and  differentiating  with  respect  to  pav  and 
simplifying  as  much  as  possible,  I  obtain 


PAV 


AV 


(a  -  b) 2aJ 


exp 


AV 


2(a-b)2a2J 


, P  AV  v  2 

(S=E) 


1  [  exD  (—  y  -  -  -  — -  —  —  — ■  ■  ■  dy 

55f  J  .  ,  ^  /yTE^M  V 


i  PAV» 

(i+H) 


.  .  . (D7) 
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Unfortunately,  the  analytical  solution  to 
Equation  (D7)  has  not  been  found.  For  the  special  case 
when  a»b,  however,  the  analysis  can  be  concluded: 


PAV  AV 


AV 


.1 


4a2^?  o2  2aal 


)  exp(- 


AV 


16a2a| 


PAV  PAV 

"l  D  3(5ir>  +  "av0  i<k£> 
*7  1  “5  1 


(D8) 


where  D^(z)  is  the  parabolic  cylinder  function.  As  a 
check,  the  integral  of  Equation  (D8)  over  pav  from  0  to  » 
is,  in  fact,  1. 

As  a  final  example,  I  consider  the  case  in  which  we 
allow  to  be  a  random  variable  and  Equation  (Al)  to  be 
the  conditional  pdf  of  p  given  a£.  Allowing  c£  to  be 
uniformly  distributed  between  b  and  a,  and  applying  Bayes 
Theorem,  I  get: 


(PAV) 


(D9) 


I 


-317- 


then: 


PPav(pAV}  "  a^E  [vT?*  “  VtIt*]  (D10) 

where  E^(z)  is  the  exponential  integral.  Applying  the 
identity  [18]: 


E1(z)  ■  -y  -  lnz  -  ,  ( |  argz  |  <  tt ) 


where  y  is  Euler's  constant  »  .5772...  I  obtain 


Pp  ^PAV 
PAV  AV 


PAV 

a-b 


CO 
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n«l 


(-1) 


n  ,pAV.n  1 

(TF}  HHT 


00 

-  I 
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(DU) 


Although  I  did  not  plot  Equation  (Dll) ,  the  numbers  I 
obtained  for  relative  levels  of  uncertainty  in  agree 


x'.'.auakui'! 
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closely  with  the  results  of  uncertainty  in  r,  and  the 
effect  of  randomization  (i.e.,  <_  3dB  uncertainty)  on 
the  amplitude  densities  appears  also  to  be  small  while 
limiting  forms  are  obtained  for  uncertainty  _>  7dB. 

A  more  thorough  investigation  of  this  phenomena 
could  uncover  the  statistical  dependencies  on 
identifiable  non-stationarities,  such  as  propagation 
loss  due  to  changing  ranges  or  other  oceanic  or 
experimental  factors  which  result  in  the  temporal 
dependence  of  r  or  c£.  The  effects  on  the  statistics 
can  then  be  analyzed  using  the  procedures  employed  above. 


/ 
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